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I.   INTRODUCTION 

A.   OVERVIEW 

3riefly  stated,  direct  time-domain  techniques  for 
transient  radiation  and  scattering  are  broadly  concerned 
with  the  response  of  some  object  of  interest,  either 
equipment  or  a  target,  to  non-monochromatic  electromangetic 
sources.  The  analysis  is  done  in  the  time-domain,  with  time 
treated  as  an  explicit  independant  variable  in_  most  methods. 

In  the  past,  frequency  domain  analysis  has  been  the 
primary  means  of  analyzing  non-monochromatic  EM  phenomena. 
Time-domain  solutions  in  closed-form  were  almost  impossible 
to  obtain,  except  in  certain  very  simple  cases.  The  amount 
of  computational  work  required  was  correspondingly 
significant.  No  practical  aides  were  available  to  assist  in 
performing  the  operations  prior  to  the  development  of  fast 
computer  systems. 

Another  limitation  was  the  requirement  that  the 
time-domain  response  be  as  nearly  an  impulse  as  possible  and 
thus  possess  a  broad  band  of  frequencies.  Hardware  was 
required  to  produce  very  short  pulse  waveforms. 
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The  situation  changed  radically  with  the  widescale 
introduction  of  computers  capable  of  performing  rapid 
computations  involving  large  numbers  of  iterations  in 
relatively  short  periods.  This,  coupled  with  the  rapid 
maturing  of  short-pulse  hardware  technology,  soon  made 
possible  the  application  of  time-domain  techniques  to  a  wide 
variety  of  eletromagnetic  problems. 

Siller  and  Landt  [Ref.  1]  list  some  advantages  of 
direct  time-domain  solutions  over  frequency  domain 
treatments  of  transient  problems.   They  are: 

1.  Greater  solution  efficiency  for  many  types  of  EM 
problems. 

2.  More  convenient  handling  of  non-linearities. 

3.  Improved  physical  insight  and  interoretablility  of  EH 
phenomenon. 

4.  Availability  of  wide-band  information  from  a  single 
calculation. 

5.  Opportunity  to  isolate  interactions,  using  time  range 
gating  (e.g.,  pulse  reflections  from  distant  objects 
not  of  interest,  etc.,). 

The  major  disadvantage  is  the  increased  complexity  of 
the  computer  implementation  of  a  time-domain  program  as 
opposed  to  one  involving  frequency  domain  techniques.   There 
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is  also  a  corresponding  difficulty  in  the  software 
development  and  its  practical  use.  Finally,  the 
computational  time  may  be  significantly  greater  than  that 
required  of  other  techniques. 

The  applications  for  time-domain  techniques  are  diverse. 
Transient  analysis  can  be  used  to  directly  evaluate  the 
transient  characteristics  of  new  antenna  designs, 
determining  such  parameters  as  driving  point  currents, 
radiated  far  fields  and  near  fields,  and  assisting  in 
determining  antenna  gain.  These  methods  may  also  be  used  in 
determining  the  degree  of  coupling  due  to  electromagnetic 
pulses  (EMP)  generated  because  of  lightening  during  a 
nuclear  detonation.  Finally,  time-domain  techniques  may  be 
directly  applied  in  the  solution  of  the  inverse  scattering 
problem  involving  the  receiving  and  scattering  properties  of 
a  large  category  of  target  types. 

As  time-domain  transient  techniques  are  further  refined, 
it  can  be  expected  that  a  greater  number  of  uses  will  be 
identified. 

B.   3EVIER  OF  TRANSIENT  EH  ANALYSIS  TECHNIQUES 

There  are  three  broad  categories  of  techniques 
identifiable  in  the  literature  involving  transient  analysis. 
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They   are    direct   (time- domain)     techniques,    transform 
(frequency      domain)       techniques,       and     hybrid 
techniques.  [Ref.  2]    This   report  is   oniy   interested   in 
direct  techniques  and  will  not  review  the  others. 
1 .   Physical  Optics  Inverse  Scattering 

The  earliest  techniques  developed  involved  physical 
optics  inverse  scattering.  Kennaugh  and  Cosgriff  [Ref.  3] 
in  1958  proposed  that  the  impulse  response  of  a  scattering, 
conducting  body  is  simply  the  second  derivative  of  the 
projected  area  function  of  the  body  if  physical  optics 
currents  are  postulated  on  the  surface.  If  a  target  is 
illuminated  with  an  incident  ramp  plane  wave,  then  the 
projected  area  function  of  the  target  is  directly 
proportional  to  the  far-scattered  ramp  response.  This  is 
only  an  approximate  relationship,  however,  since  it  is  based 
on  the  assumption  of  physical  optics  currents. 

A  more  direct  and  precise  ramp  response  technique 
was  suggested  by  Kennaugh  and  Moffat  [Ref.  4]  in  1965. 
They  predicted  satisfactory  approximate  ramp  response 
signatures  could  be  obtained  with  narrower  bandwidth 
interrogating  signals  than  that  required  for  the  impulse 
response. 
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The  above  techniques  were  generally  applicable  for 
targets  which  were  ro rationally  symmetric  and  for  whom  the 
direction  of  the  electromagnetic  field  incidence  was  axial. 
In  1976,  Young  [ Hef .  5]  demonstrated  that  a  volume  estimate 
can  be  obtained  from  the  ramp  response  waveform.  He  further 
extended  the  ramp  response  technique  to  arbitrarily  shaped 
targets  with  the  incident  field  being  other  than  axial. 
2.   "Exact"  Inverse  Scattering  Method 

Another  method  which  improves  upon  the  physical 
optics  inverse  scattering  theory  is  the  "exact"  inverse 
scattering  method  first  suggested  by  Bennett  et.  al., 
[Ref.  6]  in  1974,  and  further  refined  in  1977.  [Ref.  7] 
The  inverse  scattering  problem  is  solved  using  an  inversion 
of  a  space-time  integral.  The  method  takes  into  account  the 
backscatter  response  of  the  target,  the  direct  effect  of  the 
incident  electromagnetic  field,  and  the  correction  currents 
flowing  on  the  target  surface.  Theoretically,  by  completely 
describing  the  interactions  of  the  two  components  of  the 
surface  currents,  a  close  approximation  of  the  target  shape 
is  possible.  The  actual  implem  en  -cation  is  limited  to 
"perfectly"  conducting  scatterers.  It  is  this  method  which 
has  been  selected  for  use  in  the  Naval  Postgraduate  School 
Time  Domain  Radar  Laboratory  (TDRL) . 
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3.   Other  Methods 

Other  methods  that  have  been  developed,  although  to 
a  less  satisfactory  degree  of  accuracy  than  the  inverse 
scattering  methods,  include  radar  target  identification 
methods.  Only  a  finite  number  of  solutions  (targets  shapes) 
are  allowed  as  opposed  to  the  nearly  infinite  number 
possible  using  inverse  scattering.  The  models  of  targets 
are  stored  for  parameter  comparison  with  actual  returns. 
Certain  parametric  aspects  of  a  wideband  illuminating 
signal,  such  as  natural  resonance  freguencies  and 
polarization  of  the  return,  are  used.  Other  similar  methods 
use  amplitude  and  phase  differences  or  special  properties  of 
Rayleigh  region  scattering.  All  these  methods  are  limited 
in  practical  application. 

Still  another  important  and  powerful  method 
considers  multifrequency  target  illumination  using  existing 
radars  to  simultaneously  illuminate  targets  with  VHF,  OHF,  L 
and  S  band  radar  frequencies  with  circular  polarization. 
[  Ref .  8]  Circular  polarization  eliminates  any  polarization 
angle  dependance  the  signal  might  have.  This  approach 
allows  for  the  inclusion  of  all  radar  system  parameters 
(such  as  antenna  gain,  noise  level,  range,  etc.)  in  the 
identification     algorithm.      The    technique    further 


21 


demonstrates  that  the  target  shape  information  may  be 
obtained  with  a  high  degree  of  accuracy  from  radar  systems 
correctly  tailored  to  the  problem  of  general  target 
identification. 

C.   REVIEW  OF  TDHL  PROGRESS  AT  NPS 

Three  requirements  were  initially  identified.  First,  a 
method  had  to  be  chosen  from  the  techniques  listed  above. 
Secondly,  the  construction  of  a  physical  laboratory  to 
perform  the  required  measurements  in  support  of  the  chosen 
technigue  had  to  be  designed  and  built.  Thirdly,  an 
alogorithm  had  to  be  developed  employing  the  selected 
technique  for  laboratory  use. 

As  stated  above,  the  "exact"  inverse  scattering  method 
suggested  by  Bennett  was  the  scheme  chosen.  To  support 
this,  work  began  in  January  1980  where  evaluation  of  various 
designs  to  be  used  in  constructing  a  rime  Domain  Radar 
Laboratory  at  the  Naval  Postgraduate  School  was  carried  out 
by  Capt.  C.  W.  Hammond,  USMC,  as  a  thesis  project.  rforfc  was 
essentially  completed  with  the  erection  of  the  laboratory 
and  interfacing  of  hardware  in  September,  1980.  [Ref.  9] 

In  support  of  the  TDRL,  LCDR  :ieir  Morag,  Isreali  Navy, 
began,   as  a  thesis  project,    the  development  of  an  "exact" 
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inverse  scattering  algorithm  in  September,  1980.  The 
algorithm  was  completed  in  Sarch  1931  [ttef.  10]  However, 
the  program  was  developed  in  FORTRAN,  a  language  not  used  by 
the  TDRL  computer  system.  A.  translation  to  3A5IC  and 
adaptation  of  the  coda  was  necessary  to  reduce  memory 
storage  requirements  and  to  ansure  computational  times  were 
of  reasonable  lengths. 

The  final  step  in  the  establishment  of  the  TDRL  was  the 
validation  of  the  laboratory  in  confirming  the  work  of 
Hammond  and  Morag.  This  is  one  of  the  principal  aims  of 
this  project. 

D.   THESIS  OBJECTIVES 

The  objectives  of  this  thasis  are  to: 

1.  Develop  a  working  software  program  which  will  acquire 
and  prepare  transient  response  data  from  various 
targets  for  input  to  the  time-domain  integral 
equation. 

2.  To  show  that  through  signal  processing,  the  effects 
of  the  receiving  antanna  can  be  neglacted,  and  the 
response  to  arbitrary  specified  exciting  waveforms 
can  be  constructed,  or  "synthesized",  using  FFT 
techniques. 
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3.   Validate  the  TDRL. 

The  main  effort  has  been  directed  toward  the  development 
of  a  user  oriented  computer  program  which  will  allow  the 
determination  of  the  second  objective  and  prepare  the 
transient  domain  data  for  input  into  the  numerical 
time-domain  integral  equation.  The  shape  of  the  object 
being  observed  is  displayed  both  before  and  after  correction 
factors  are  applied.  The  shape  of  the  object  was  developed 
from  the  ramp  response  of  the  impulse  function  applicable  to 
the  target.  The  program  was  also  made  compatible  with  the 
input  of  data  to  an  alternate  natural  resonance  frequency 
target  identification  technique,  using  Prony's  method,  which 
is  being  developed  by  LT  Demetrius  Papaspiridakos,  Hellinic 
Havy,  as  a  thesis  project. 

This  report   first  describes   the  physical   structure  of 

the  TDRL.   It  then  derives  the  relevant  time-domain  integral 

equations   for  both   the  radiation   and  inverse   scattering 

solution  for   simple  axisymmetrical,    perfectly  conducting, 

closed  surfaces   in  Chapter  III.     tfext,   it   discusses  and 

develops  the  numerical  solutions   to  the  integral  equations. 

In  Chapter  V,    the  implementation  of  the   complete  computer 

code   is  presented.     And  finally,    measurements  that   are 

relevant  to   the  validation  and  confirmation   of  previously 

developed  concepts  are  reported. 
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This  thasis  is  the  third  to  be  produce!  for  this  major 
project  at  the  Naval  Postgraduate  School.  It  marks  a  point 
of  departure  in  that  essentially  all  the  preliminary  work 
which  will  allow  for  the  effective  use  of  the  TDRL  has  been 
performed  with  validation.  The  laboratory  is,  for  the  most 
part,  ready  for  applied  research  involving  transient 
time-domain  analysis. 
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II.   PHYSICAL  DESCRIPTION  07  TDRL 

A.   INTRODUCTION 

As  stated  in  Chapter  I.,  the  primary  thrust  of  this  work 
is'  to  provide  a  practical  working  alogorithm  to  serve  as  the 
principal  link  between  the  acquisition  of  transient  data  for 
targets  situated  on  the  TDRL  scattering  range,  and  the  use 
of  that  data  in  the  inverse  scattering  algorithm. 

The  work  performed  on  this  project  is  a  direct  extension 
of  the  preliminary  work  performed  by  Hammond  in  the 
construction  of  the  Naval  Postgraduate  School  Time  Domain 
Radar  Laboratory  (TDRL) ,  and  by  Morag  in  the  development  of 
the  numerical  integral  equation  algorithm  for  the  solution 
of  the  inverse  scattering  problem  involving  axisymmetric 
targets. 

The  results  obtained  by  Hammond  and  Morag  are 
essentially  complementary.  However,  the  independant 
approaches  each  took  in  their  work,  particularly  in  the 
selection  of  unrelated  processing  machines  for  data 
processing,  made  necessary  certain  practical  modifications 
to  the  TDRL.  It  is  germane  to  review  the  basic  work  of 
Hammond  and  discuss  those  modifications  mai=.   Such  a  review 
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will  also  allow  for  the  setting  of  a  frame  of  referance  to 
aid  in  the  explanation  of  the  linking  algorithm  developed 
and  results  obtained  in  this  project. 

3.   PHYSICAL  DESCRIPTION  OP  TDRL 

The  physical  set-up  of  the  TDRL  is  shown  in  Figure  2. 1. 
As  can  be  seen,  the  TDRL  can  be  broken  functionally  into 
three  groups: 

1.  The  Impulse  Generator  Source  Group 

2.  The  Imaging  Plane  Group 

3.  The  Sampling  Receiver  and  Signal  Processing  Group 
Basically,   the  TDRL  is  an  open  system  designed  to  determine 
the  transfer   function  of   a  target   from  its   backscattered 
radiation.     Scattering   measurements  are   bi-static.     A 
description  of  each  functional  group  follows. 

1 .   Source  Group 

This  group  consists  of  a  single  itea,  a  Video 
Communications  UHF  Impulse  Generator,  Model  1000.  This 
generator  is  the  source  of  the  transmitted  pulse  and  also 
provides  the  triggering  signal  which  synchronizes  the  Signal 
Processing  Group  to  the  transmitted  signal. 

The  Impulse  Generator  is  specified  to  provide  a 
fixed  amplitude  45-   volt  pulse  with  a  maximum   rise  time  of 
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Figure  2.3.   Source  Pulse  froa  Generator 

370  picoseconds.  Figure  2.3  is  a  typical  waveform  as 
measured  on  a  Tektronix  S-6  Sampler  Head  by  the 
manufacturer . 

The  output  of  the  Impulse  Generator  as  measured  at 
the  TDRL  as  measured  by  the  Operating  System  is  given  in 
Figure  2.4,  and  Figure  2.5  as  photographed  from  the  DPO 
directly.  To  set  full  waveform  presentation  on  the  measuring 
oscilloscope,  30dB  attenuation  was  applied  to  the  generator 
output.  Taking  the  attenuation  into  consideration,  the 
measured  peak  value  of  voltage  (1.4604  volts)  gives  the 
actual  output  of  the  generator  to  be  46.2  volts.  This  is 
about  3%   better  than  rated.    The  rise  time  of  the  output  is 
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TARGET:  RUN  DATE     DIST    ANT  TGT   REMARKS 

1   6-26-81   N/A    H9   ???   GENERATOR  WAUEFORM 


HAXIMUM  PEAK  UALUE  (UOLTS) 
MINIMUM  PEAK  UALUE  (UOLTS) 
RMS  UALUE  (UOLTS) 
MEAN  UALUE  (UOLTS) 
NUMBER  OF  WAUEFORMS  AUERAGEO  ■  I 


1.1284 

-9.3328 

9.4214 

8.9888 

OPTIMIZATION  UALUE  «  8 


1.86 


8.38 


8.88 


DIRECT   UAUEFORM 


E-9  SEC 


9.39 


1.88 


1.38 


2.98 


Figure  2.4.      Source   Waveform   Measured   at   TDRL 

about    250      picoseconds,      also    better      than   specified      by   the 
manufacturer,    by   about    48%. 

Some  notable  features  of  the  generator  used  in  TDRL 
measurements  are  a  variable  trigger  repetition  rate  from  10 
Hz  to  1  MHz  and  an  adjustable  delay  from  50  to  150 
nanoseconds.  Ml  outputs  are  inatched  into  a  50-ohm  load. 
The  generator  allows  rapid  acquisition  of  highly  coherent, 
repetitive   waveforms.  This    is      an    important      feature   when 

attempting   to   maximize    signal-to-noise      ratio    through   signal 
averaging. 
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Figure  2.5.   DPO  Waveform — OHF  Iapulse  Generator 


The  Video  Communications  Impulse  Generator  replaces 
a  Tektronix  Type  109  impulse  generator  used  by  Hammond.  The 
completely  solid-state  electronics  of  the  Mark  1000  has 
proven  more  reliable  in  producing  a  relatively  noiseless, 
pulse-to-puise  coherent  signal  than  did  the  reed-switched 
Type  109.  The  Type  109  had  a  fixed  pulse  repetition  rate  at 
about  720  pulses  per  second.  This  is  a  very  limiting  factor 
when  attempting  to  noise  average  a  signal. 

The  impulse  source  is  connected  to  the  transmitting 
antenna  via  a  3.2  meter  RG-SA/U  coaxial  cable.  This  cable 
serves  both  as  a  conductor  of  the  source  pulse  and  a  short 
delay  element  to  the  transmitted  signal. 
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2.      Imaging    Plane    Group 

The   Imaging   Plane    Group    consists   of    four   elements 
the      transmitting      antenna,      the      receiving      antenna,         the 
target,    and  the    imaging  plane. 

a.  Transmitting  \ntsnna 

The  transmitting  antenna  is  a  6.4-meter  endfed 
monooole.  It  extends  vertically  from  the  ground  plane  and 
produces  a  vertically  polarized  EM  field. 

b.  Receiving  Antenna  —   Transverse  Electromagnetic 
Horn  (TEH) 

Hammond  explored  a  variety  of  receiving  antenna 
types  to  be  used  in  the  TDRL.  The  most  useful  is  the  TEM 
horn  antenna  in  that  it  has  the  following  advantages  over 
the  other  antennas  considered: 

1.  Better  directivity. 

2.  Less  Attenuation. 

3.  Mobility  on  the  image  plane. 

4.  A  readily  variable  characteristic  impedance. 

5.'   Relatively  frequency  indspendant  over  a  wideband. 

The  TEM  Horn  usad  in  the  TDRL  is  physically  an 
open  sided  pyrimidal  structure  when  deployed.  It  is  a 
single  piece  of  flat  copper,  about  1mm  thick  and  triangular 
in  shape.    Figure  2.5  gives  the   relevant  dimensions.    The 
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Figure  2.6.   TEH  Horn  Dimensions 

hole  at  the  apex  of  the  antenna  allows  for  connecting  to  the 
receiving  cable. 

The  receiving  horn  antenna  is  mounted  in  a  three 
dimensional  trapezoilal  support  mount  with  quarter-inch 
plexiglass  walls  on  the  top,  bottom,  and  sides.  The  front 
and  rear  are  open  to  allow  easy  insertion  of  the  copper 
plate.  The  dimensions  and  shape  of  the  mount  are  given  in 
Figure  2.7   All  dimensions  are  outside  wall  measurements. 

The  flat  copper  antenna  is  slipped  into  the 
front  of  the  plexiglass  mount  and  protrudes  through  the 
back.    k   coaxial  cable  is  connected   to  the  back  with  an  rf 
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c.       FRONT 
Figure  2.7.      rEM  Horn   Support    Mount   Dimensions 

end-connector  and  hand-tightaned  screw.  See  Figure  2.13  for 
a  diagram  of  the  assembly  of  the  connection  point.  The 
antenna    is   physically    supported    within      the    mount   by   grooves 
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cut  into  and  running  along  the  length  of  each  side.  The 
grooves  serve  to  fix  the  flare  angle  of  the  antenna  at 
approximately  24.10. 

The  dimensions  of  the  mount  were  chosen  for 
reasons  of  practicality,  being  the  most  appropriate  for  the 
material  available  and  providing  a  reasonable  flare  angle 
and  antenna  characteristic  impedance  for  signal  reception. 
The  mount  also  provides  stabilty  to  the  Horn  under  various 
weather  conditions.  Thus,  by  placing  the  horn  in  the 
plexiglass  support  mount,  physical  and  electrical 
characteristics  are  fixed,  assisting  in  correlating  run  to 
run  measurements.  The  T2M  Horn  as  it  would  be  mounted  on 
the  imaging  plane  is  drawn  in  Figure  2.8.  The  support  mount 
is  not  shown  for  clarity. 

The  addition  of  the  plexiglass  mount  should  have 
negligible  effect  on  the  received  signal.  Figure  2.9  shows 
the  response  of  the  supported  horn  antenna  to  a  pulse 
generated  by  the  source.  This  compares  favorably  with  the 
pulses  obtained  by  Hammond  without  the  plexiglass 
mount.  [Hef .  1 1 ] 

As  noted  in  Figure  2.8,  the  electric  field  at 
the  aperture  of  the  horn,  assuming  a  uniform  TEM  wave  being 
excited  between  the  horn  and  the  ground  plane,  is  given  by: 
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Figure  2,8.   TEH  Horn  Receiving  Antenna 
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V  -jkr 

E    =   IT     e  (2.1) 


where:  V0  =  amplitude  of  the  incident  voltage  at 
the  feed    point. 

h      -   height   of   aperture 

k      =  wave   number    =   2 7T/  \ 

r0    =  distance        of      observer        from        the 
referance   point. 

It    can   than      be   easily    shown   that      the   receiving 

antenna   impulse   response   in   the   far    field   is:    [Hef.    12] 

Sr(t)    =    -    c    £      WZT    -5     (t    -    r   /c)  (2.2) 

OLo 

where:         *     =   impulse 

Z^  =  load  impedance 

and  the  impulse  response  of  the  transmitting  antenna  is: 

t  60e  W 

S*(t)  =  -  -^  £  5  (t-ro/c) 

where:  r  =  distance  from  the  aperture  to  an 
observer  on  the  image  plane  in  the 
boresight  direction. 

Thus  the  transmitting  impulse  response  is  simply 

the  time  derivative  of  the  receiving  impulse  response.    The 

impulse  response  is   a  delayed  impulse,   delayed   in  time  by 

r0/c,    and  multiplied  by  a   negative  constant  which   is 

dependant  on  the  width  of  the  antenna  aperture,   w,   and  the 

load  impedance  Zu . 
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TARGET:  RUN  DATE     D1ST   ANT  TCT  REMARKS 

1   6-26-81   N/A    Hfl   ???   GENERATOR  WAVEFORM 


MAXIMUM  PEAK  VALUE  CUOLTS) 
MINIMUM  PEAK  UALUE  <UOLTS> 
RMS  VALUE  <UOLTS) 
MEAN  VALUE  <VOLTS> 
NUMBER  OF  WAVEFORMS  AVERAGED  -  1 


9.2328 

-8.0932 

0.8796 

8.0888 

OPTIMIZATION  VALUE  -  8 


E-1U0LT 


DIRECT  WAVEFORM 


E-9  SEC 


8.59  1.80  i.38  2.08  2. 50  3.00  3. 30  4.00  4.30  3.00 

Pigure  2.9.   Typical  TES  Horn  Received  Waveform 

Generally,  the  following  limitations  are 
applicable  to  the  TEM  Horn  during  most  TDRL 
measurements:  [Ref.  13] 

1.  The  antenna  is  not  in  the  far  field  of  the 
transmitting  antenna. 

2.  The  flare  angla  of  ths  TEM  Horn  is  not  wide  enough  to 
receive  energy  from  all  positions  along  the 
transmitting  antenna. 

3.  Wavelengths  at  high  frequencies  may  not  propagate  in 
a  true   TEM   mode,    resulting   in   attenuation   and 


distortion. 


39 


The  net  result  is  a  deviation  from  the  ideal 
results  of  the  impulse  responses  as  stated  above.  When 
Figure  2.9  is  compared  with  Figure  2.4,  distortion  of  the 
pulse  as  received  by  the  antenna  is  noted.  However,  for 
TDRL  purposes,  the  reproduced  waveform  is  still  reasonably 
faithful  to  the  original  source  pulse. 

c.   Receiving  Antenna  --  Monopoles 

Monopole  antennas  are  also  used  on  the  TDRL 
scattering  range.  One  of  eleven  monopoles,  ranging  in 
length  from  85  millimeters  to  285  millimeters,  in  20 
millimeter  lengths,  can  be  chosen.  Each  monopole  is 
constructed  of  brass  rod  and  is  threaded  on  one  end  to  allow 
for  ease  of  connection  to  the  receiving  elements. 

The  monopoles  may  be  considered  center-fed 
dipoles  with  their  images  in  the  ground  plane.  In  effect, 
they  are  oscillating  electrical  dipoles  in  free  space. 
Oscillation  at  frequencies  related  to  the 
quarter-wavelengths  of  the  monopoles  is  to  be  expected. 

Figures  2.10,   2.11  and   2.12  confirm  the  above. 

The  effects  of  radiation  damping,   dispersion  and  reflection 

due  to   the   impedance  discontinuity   of   the   antenna   are 

clearly  demonstrated.    From   Figure  2.4  the  pulse   width  of 

the  source   pulse  is  seen  to  be  about  1.2  aanoseconds.   This 

corresponds  to  a  wavelength  of  about  360  millimeters. 
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target:  run  date    dist   ant  tgt  remarks 

4  s-29-81  1.2711  2  ???  acquire  85mm  antenna  response 


MAXIMUM  PEAK  UALUE  <UOLTS> 
MINIMUM  PEAK  UALUE  <UOLTS> 
RMS  UALUE  <UOLTS) 
MEAN  UALUE  <UOLTS> 
NUMBER  OF  WAUEFORMS  AUERAGED  ■  1 


8.9577 

-8.8782 

8.8237 

8.8888 

OPTIMIZATION  UALUE  •  9 


DIRECT  HAUEFORM 


5.38 

8.88 

-5.881. 

E-9  SEC 

8. 58  1.88  1.58  2.88  2.58  3.89  3.59  4.99  4.59  5.99 

Figure  2«10.   85aa  Monopole  Pulse  Response 

Figure  2.10  is  the  short-pulse  response  of  the 
35mm  monopole.  This  antenna  is  almost  exactly  one 
quarter-wavelength  in  height  in  relation  to  the  source 
pulsewidth.  It  correspondingly  produces  a  nearly  pure 
sinusoidal  variation,  damped  in  time,  as  predicted  by 
classical  antenna  theory. 

Figures  2.11  and  2.12  are  the  observed  responses 
for  the  185  millimeter  and  285  millimeter  monopoles.  As  can 
be  seen,  the  antennas  tend  to  resonate  at  their 
guarter-wavelength  periods  of   2.5    and   3.8   nanoseconds 


41 


TARGET:  RUN  DATE     DI3T   ANT  TCT  REMARKS 

3  6-29-81   1.27M  7   ???  ACQUIRE  183NM  ANTENNA  RESPONSE 


MAXIMUM  PEAK  UALUE  (UOLTS) 
MINIMUM  PEAK  UALUE  (UOLTS) 
RMS  UALUE  (UOLTS) 
MEAN  UALUE  (UOLTS) 
NUMBER  OF  WAUEFORMS  AUERAGED  «  1 


8.0859 

-8.8987 

8.8328 

8.8888 

OPTIMIZATION  UALUE  ■  8 


DIRECT  WAUEFORM 


E-2U0LT 


3.88 


8.88,, 


-5.88 


E-8  SEC 

8.89     8.28     8.48     8. 69     9.88     1.89 
Figure  2.11.   185nm  Sonopole  Pulse  Response 

respectively.     Although   the   amplitudes   of   these   two 
aonopoles   do  not  differ  greatly   (less  than  .08%  between 
their  maximum   peak  values) ,    the  longer   antenna  tends   to 
distort  the  incident  wave  to  a  greater  extent.   This  is  most 
readily   observed  when   comparing   the   root-mean-square 
voltages  of  the  two  waveforms.     The  185  millimeter  antenna 
provides  more  effective  power  to  the  receiving  circuits  than 
does  either  of  the  other  two  moaopoles.    The  285  millimeter 
antenna  dissipates  a  considerable  amount   of  energy  in  the 
long  period  oscillations;   the  85   millimetar  antenna  is  too 
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small  to  acquire  the  full  power  of  the  transmitting  antenna. 
Thera  may  then  be  a  "best"  oionopole  length  for  a  particular 
source   pulse.  What    this   length      might  be   and      whether   the 

monopoles  are  selective  of  targets  according  to  their 
dimensions    will   be  discussed    more   fully  in      Chapter  7. 

TARGET:    RUN  DATE  DIST        AHT  TGT     REMARKS 

2      6-29-81      1.27H      12      ???      ACQUIRE   293MM   ANTENNA  RESPONSE 


MAXIMUM   PEAK   VALUE   (UOLTS) 
MINIMUM  PEAK   UALUE    (UOLTS) 
RMS   UALUE    (UOLTS) 
MEAN   UALUE    (UOLTS) 
NUMBER   OF   HAUEFORMS   AUERAGED   ■    1 


9.8861 
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9.9289 
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DIRECT   WAUEFORM 
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9.88 


8.88 


-3.88 


8.99 


8.38 


1.89 


1.39 


2.99 


E-8  SEC 


Figure    2.12.       285a«  Honopole   Pulse   Response 

All  receiving  antennas  are  coupled  to  the  Signal 
Processing  Group  through  a  coaxial  end-coupling  fitted 
through     the      imaging      plane.  The      connection      point      is 

diagrammed      in    cut-away      in   Figure      2.13.  The   antenna      is 
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Figure  2.13.   Receiving  Antenna  Connection  Point 
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coupled  as  close  to  the  ground  plane  as  possible.  This  is 
particularly  important  foe  the  TEM  Horn  to  prevent 
undesirable  oscillations.  As  drawn,  the  threaded  portion  of 
the  monopole  (or  hand  tightened  screw  of  the  TEH  Horn)  is 
joined  to  a  connector  pin  by  butting  the  surfaces  of  the 
antenna  and  pin  together.  A  coaxial  cable  is  then  connected 
to  the  end  connector  and  the  rf  signal  is  transferred  to  the 
receiving  circuits. 

This  arrangement  has  worked  well.  No 
reflections  or  excessive  attenuation  of  the  incoming  signal 
due  to  the  joint  were  observed.  Oscillations  on  the  TEM 
Horn  were  minimal. 

Table  I   lists  the  antennas  currently  available 
to  the  TDRL  user.    The  significant  dimensions  and  pertinent 
characteristics  are  reviewed, 
d.   Targets 

The  targets  used  on  the  TDRL  scattering  range 
are  of  various  dimensions.  Currently,  eleven  are  used  in 
imaging  experiments.  They  are  all  axisyaietric,  closed 
surfaces  with  edges  and  vertices.  The  targets  may  be  broken 
into  three  catergories  according  to  general  shape: 

1.  Half-cylinders 

2.  Half-spheres 
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TABLE    I 


TDRL   Antennas 


ANTENNA 
NUMBER 

DESCRIPTION 

DIMENSIONS 

COMMENTS 

1 

TEM  Horn 

Figure  2.6 

Broadband 

2 

Dipole 

85mm 

Oscillating  Dipole 

3 

Dipole 

105mm 

4 

Dipole 

125mm 

5 

Dipole 

145mm 

6 

Dipole 

165mm 

1 

7 

Dipole 

185mm 

8 

Dipole 

205mm 

9 

Dipole 

225mm 

10 

Dipole 

245am 

11 

i 

Dipole 

265mm 

12 

Dipole 

285mm 

1 

3.       Half- cones 
Table    II   lists   the  targets      and   provides   general   information 
about    each. 

There  are  six  half-cylinders  used  as  targets. 
They  are  all  machined  from  aluminum  stock.  The  largest  is 
12-inches  long  with  a  3-inch  radius.  The  smallest  is 
3.06-inches   long   and    1.5-inches    in   radius. 

A  single  half-cone  is  also  used  as  a  target. 
The  cone  is  machined  from  the  same  stock  as  the 
half-cylinders. 
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Four  spheres,  varying  in  radius  from  5-inches  to 
1. 5-inches  are  the  four  remaining  targets.  The  spheres  are 
basically  styrofoam  balls  cut  in  half.  Sach  half  is  covered 
with  heavy-duty  aluminum  foil  to  provide  the  proper 
reflection   of   transient   EM. 

Sach  of  the  targets,  when  placad  on  the  ground 
plane,  appears  to  the  incident  pulse  to  be  mirrored. 
Therefore,  the  receiving  antenna  sees  not  the  half-shapes, 
but  the   full  cylinder,    cone   cr   sphere. 

No  particular  problems  were  noted  when  using  the 
targets  except  for  the  need  to  anchor  the  spheres  with 
weights   during   inclement   weather.  Otherwise,      they    have  a 

tendency   to    move   on   the   plane   due   to    the   wind, 
e.      Image  Plane 

The  final  element  of  the  Image  Plane  Group  is 
the  image  plane  itself.  It  is  a  square  structure  of  welded 
Reynolds  Type   AN-190    aluminum   sheet.  It    measures   about    11 

meters   on    a   side   for   an  area   of    about    121    square   meters.      It 
is    four  millimeters   thick. 

Based  on    the  above,    the   range   clear   time,      which 

establishes    the   low   frequency   cut-off   of  the    measurement,    is 

greater   than      36   nanoseconds.  The   clear      time   establishes 

■che   first   arrival     at    the    measurement    point      of   a   reflection 

from  objects   surrounding  the   plane   and   the    end    of   the    plane. 
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TABLE  II 


TDRL  Targets 


TARGET 

i i 

TYPE 

i — 

RADIQS 

LENGTH 

MATERIAL 

1 

Cylinder 

3" 

12" 

I 

Aluminum  Stock 

j 

2 

Cylinder 

1.5" 

12" 

Aluminum  Stock 

j 

3 

Cylinder 

3" 

6" 

Aluminum  Stock 

4 

Cylinder 

1.5" 

6.0  6" 

1 

Aluminum  Stock 

L 1 

5 

Cylinder 

3" 

3.06" 

Aluminum  Stock 

6 

Cylinder 

1.5" 

3.06" 

Aluminum  Stock 

7 

Sphere 

r~ 

5" 

Al  and  Styrofoam 

8 

Sphere 

tl« 

Al  and  Styrofoam 

9 

Sphere 

3"      ; 

Al  and  Styrofoam 

._ .... 

10 

Sphere 

1.5" 

Al  and  Styrofoam 
i                   j 

11 

Cone 

3.18" 

..,„ 

5.94"     Aluminum  Stock 

L  _         . 1 

The  main  purpose  of  the  image  plane  is  to 
isolate  the  targets  and  tha  antennas  from  the  underlying 
structure  of  the  TDRL.  Instrumentation  cables  are 
introduced  from  under  the  ground  plane  and  are  thus  kept 
from  interfering  with  the  measurement.  The  plane  also 
serves  as  the  primary  support  structure  for  both  the 
transmitting  and  receiving  antsnna. 

No  significant  problems  were  encountered  with 
the  image  plane.  However,  soma  disadvantages  of  this  plane 
arrangement  are  that  only  objects   having  a   symmetry  plane 
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can  be  measured,  and  polarization  and  incident  angles  are 
limited.  But,  for  the  work  intended,  the  image  plane  serves 
its  purpose  well. 

3.   Sampling  Receiver  and  Sicjnal  Processing  Group 

The  Signal  Processing  Group  consists  of  three 
elements: 

1.  Digital  Processing  Oscilloscope  (DPO) 

2.  Mini-Computer 

3.  Hard  Copier 

The  Signal  Processing  circuits  are  connected  from 
the  output  of  the  receiving  antenna  to  the  input  of  the  DPO 
via  2.58  meters  of  RG-8A/U  coaxial  cable.  The  cable 
provides  some  attenuation  and  delay  of  the  source  signal 
provided  by  the  impulse  generator.  Additionally,  a  trigger 
signal  is  sent  to  the  DPO  from  the  pulse  generator  via  .54 
meters  of  RG-53C/U  cable. 

The  lengths  of  the  signal  and  trigger  cable  are 
chosen  for  convenience  and  for  minimum  attenuation.  They 
may  be  varied  for  special  purpose  applications  such  as  to 
accommodate  a  different  trigger  source  by  providing  a 
necessary  delay.  Care  must  be  taken  not  to  introduce  too 
much  delay  or  the  received  signal  will  not  be  windowed  on 
the  CRT  display  of  the  DPO.   The  maximum  practical  length  of 
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cabling,  including  that  to  the  transmitting  antenna  from  the 
pulse   generator,   should   bs  less   than  20   meters  for   the 
current  equipment  configuration  of  the  TDRL. 
a.   Digital  Processing  Oscilloscope 

The  Tektronix  Digital  Processing  Oscilloscope  is 
a  computer  compatible,  analog  to  digital  system.  It  is 
comprised  of  two  elemental  units;  a  Tektronix  7704A  General 
Purpose  Oscilloscope  System,  and  a  P7001  Processor. 
Modification  of  the  7704A  Oscilloscope  by  the  addition  of 
various  plug-in  modules  allows  for  tailoring  of  the  whole 
DPO  system  to  the  specific  needs  of  the  IDRL.  A  cursory 
description  of  the  various  elements  of  the  DPO,  as 
applicable  to  the  TDP.L,  follows. 

The  7704A  Oscilloscope  system  is  composed  of  two 
parts;  the  D7704  Display  Jnit  and  the  A770U  acquisition 
unit.  The  Display  unit  provides  results  of  signal 
acquisition,  both  real  time  and  as  stored  in  the  P7001.  It 
does  this  by  a  visual  interface  with  the  DPO  in  the  form  of 
an  Oscilloscope.  Additionally,  the  CRT  readout  is  carried 
over  into  the  DPO  to  provide  appropriate  scaling  on  the 
displays. 

The  A770£*  acquisition  unit,  through  its  modular 
plug-ins,   provides  the  means  to   acquire  the  analog  signal. 
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The  plug-in  units  used  in  the  TDRL  are  the  s-6  Sampling  Head 
and  the  S-53  Trigger  Recognizer,  inserted  into  a  7S12 
TDR/Sampler. 

The  7S12  provides  for  measurement  of  recurring 
fast-rise  time  signals;  i.e.,  the  source  pulse  and  the 
target  back-scatter.  It  features  the  means  of  determining 
vertical  deflection  factors  in  millivolts;  horizontal  sweep 
factors  in  seconds;  a  HIGH  RESOLUTION  switch  which  reduces 
the  waveform  noise  and  jitter  by  signal  averaging  (10 
averages  per  sample  point) ;  a  TIMS-DISTANCE  scale  which 
allows  one-way  distance  measurements  in  air  dielectrics;  and 
a  locate  switch  which  increases  the  time/division  and 
intensifies  a  portion  of  the  display  to  locate  the  time 
window  relative  to  the  total  waveform.  Finally,  the  7S12 
provides  the  means  for  the  inserting  and  interfacing  of  the 
plug-in  heads  into  the  DPO. 

The  S-6  Sampling   Head  is  a  50-ohm   loop  through 

input  sampling  unit.   It  provides  the  means  of  inputting  the 

signal  of  interest   to  the  DPD .    This  unit   has  a  specified 

rise   time   of   30   picoseconds   or   less   and   a  bandwidth 

equivalent  to  dc  to  11.5  GHz  at  3  dB  down.    It  has  a  signal 

range   of   + 1volt  to   -1volt  (dc   to   ac   peak)    or   1volt 

peak-to-peak.    The   S-6  establishes  the   baseline  measuring 

capabilities  of  the  TDRL. 
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The  S-53  Trigger  Recognizer  permits  use  of  the 
7512  as  a  general  purpose  sampler.  It  produces  a  stable 
trigger  signal  from  input  signals  from  DC  to  1  GHz.  It  is 
capable  of  recognizing  signals  having  amplitudes  ranging 
from  10  millivolts  to  2  volts  peak-to-peak  into  50  ohms.  A 
trigger  pulse  from  the  pulse  generator  is  fed  into  this  unit 
at  the  same  rate  as  the  source  signal.  Trigger  to  signal 
delay  is  specified  as  15  nanoseconds  or  less. 

The  second  element  of  the  DPO,  the  P7001 
Processor,  provides  the  link  from  the  DPO  to  the 
mini-computer.  This  unit  receives  the  analog  signal 
acguired  from  the  Acquisition  unit  and  digitizes  and  stores 
or  discards  the  information,  along  with  appropriate  scale 
factors,  as  selected  by  the  operator.  The  method  used  by 
the  digitizer  is  pseudo- randoi  sampling  with  one  sample 
taken  every  6.5  microseconds.  This  is  the  maximum  rate  of 
sampling.  Any  transient  longer  than  5  milliseconds  or 
COHERENTLY  REPETITIVE  signal  up  to  the  frequency  limit 
specified  by  the  Sampling  Head  that  can  be  displayed  on  the 
D7704  Display  Unit  can  be  stored  in  internal  memory  along 
with  its  scale  factors  .  The  information  stored  can  then  be 
redisplayed  on  the  CRT  and/or  sent  to  the  minicomputer  for 
processing.     There  are   a  maximum   of  512   samples   taken 
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horzontally  (9-bit  word)  per  waveform.  Although  either  X- 
or  T-  (TIMS)  may  be  specified  for  this  axis,  the  TDRL 
primarily  employs  the  time  axis.  Vertical  resolution  of  the 
Processor  is  1024  levels  (10-bit  word).  This  axis  is  used 
to    measure   the   voltage    of    the   back-scattered   signal. 

b.  Tektronix  4052  Mini-Computer 

The  minicomputer  used  in  the  TDRL  is  a  Tektronix 
4052.  It  is  a  high  performance,  integrated  graphics  system 
employing  LSI,  bi-polar,  16-bit  technology.  The  unit  has  a 
64k-byte  memory  with  a  300!c-byte  magnetic  cartridge  built 
in.  Extended  BASIC  is  the  high  level  programming  language 
employed  to  perform  processing  functions.  Although  this 
unit  has  many  features  worthy  of  note,  the  above  are  the 
salient   ones   for   the   TDRL. 

c.  Hard   Copier 

The  Hard  Copier  used  is  the  Tektronix  4631  Hard 
Copy  Unit.  It  provides  permanent,  dry  copies  of  the  graphic 
and  alphanumeric  information  displayed  on  the  CRT  storage 
screen  of  the  4052  mini-computer. 

All  units   described  for   the  Signal   Processing 

Group  were  especially   selected  to  perform  as   an  integrated 

system  requiring  minimum  interfacing.    The   GPIB  bus  is  the 

means  of  transferring  information  from  one  unit  to  the  other 

in  the  TDRL. 
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Ill,   THEORY 

This  chapter  is  devoted  to  the  general  derivation  of  the 
equations  which  are  applicable  to  the  time-domain  analysis 
of  the  impulse  and  camp  responses  of  axisymmetric,  metallic 
bodies  which  are  illuminated  by  an  axially  directed  incident 
electromagnetic  field.  Where  sources  providing  detailed 
derivations  are  numerous  and/or  adequate,  the  derivations  in 
this  chapter  will  be  general,  providing  only  the  major 
results  of  interest  with  specific  constraints  noted. 

The  time-domain  integral  equations  which  provide 
solutions  for  the  transient  electromagnetic  problem  will  be 
considered  first.  They  will  be  followed  by  a  discussion  of 
the  measurement  equations  and  procedures  employed  in  the 
TDRL.  Finally,  various  methods  used  to  reduce  the  noise 
error  in  acquired  signals  will  be  discussed. 

A.   TIHE-DOHAIN  INTEGRAL  EQUATIONS 

One  primary  problem  in  transient  analysis  of  target 
impulse  responses  is  to  determine  the  shape  of  an  unknown 
object  when  the  incident  field  and  the  response  of  the 
scatterer  to  that  field  are  a.  priori  information.  Numerous 
time-domain   techniques  have   been   developed  which   perform 
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this  function.  As  noted  in  Chapter  I,  tha  earliest  method 
used  involved  the  physical  optics  response  of  a  body  to  an 
incident  electromagnetic  field.  The  surface  currents 
produced  an  approximate  impulse  response  that  was 
essentially  the  second  derivative  of  the  projected  area 
function  of  the  scatterer.  [Bef.  5]  However,  this  approach 
was  limited  in  its  accuracy  in  that  the  relationship  between 
the  impulse  response  and  the  derivatives  is  exact  only  at 
the  leading  edge  of  the  scattered  field  response,  a  single 
point  in  time.  Interactions  between  currents  on  the  target 
body  which  continue  to  radiate  for  a  significant  time  after 
the  incident  field  leading  edge  has  propagated  further  in 
space  alters  the  backscattsred  response  of  the  target. 
Therefore,  for  an  accurate  solution  to  the  problem,  these 
additional  currents  must  be  of  consideration. 

Bennett  [Ref.  14]  first  proposed  a  direct  method  which 
"corrected"  the  response  by  computing  the  currents  flowing 
on  the  scatterer  surface.  From  these  currents,  the 
scattered  field  could  be  calculated  "exactly". 

In  general,  the  remainder  of  this  section  follows  the 
format  of  the  method  used  by  Bennett  in  the  derivation  of 
integral  equations  for  exact  solutions. 
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1 •   Exact  Integral  Solution 

The  expression  for  the  lagnetic  field  H  at  an 
arbitrary  point  in  space  (not  on  the  scatterer  surface)  is 
given  by: 


H(r 


?.t)   -H^r-.t)   *^    /s    i[-\+  jfc£]  JCr'.t)  x  ^     ds '        (3.1) 


5         R 

where:  H(r,t)   =  total  magnetic  field  at  (r,t) 

V   (?, t)  =  incident  magnetic  field  at  (r,t) 
J(r»,t)  =  surface  current  at  (?',t) 

r    =  position     vector    to     the 

observation  point 
"r1    =  position     vector    to     the 

integration  point 
R    =  |c-r'  | 


A         r   -  r' 

**  "ft— 


—      t  =    time  in        light-meters  (one 

light-meter   is    the      time  it    takes 

light   to   travel   one    meter) 

T      -    t-R   =    retarded  time 

3y   specializing  the     arbitrary   space   point   to     "r,      a 

point    on   the   surface   of  the   scatterer,    and   then    applying  the 

boundary  conditions    (i.e.,    by   causuality,      incident   field   is 
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zero  prior  to  arrival  of  the  incident  pulse  at  a  sample 
point) ,  an  integral  equation  for  the  curreat  density  on  the 
surface  of   the   scatterer   can   be   found: 

J(?,t)    =    zln  x  H^r.t)    ♦  £    rs  -an  x    {[i  +  |  ±]    J(?,x) 

K. 

x   aR)    ds'"  (3.2) 

— * 
where:   r  =  position  vector  to  a  point  located  on 

the  surface  of  the  scatterer 

an  =  unit  vector  normal  to  the  surface 

The  first  term  on  the  right-hand  side  of  equation 
(3.2)  is  the  source  term.  It  represents  the  direct 
influence  Df  the  incident  field  on  the  current  at  the 
observation  point  (r,t).  It  is  the  physical  optics 
approximation  for  the  surface  current  when  applied  to  the 
illuminated  side  of  the  scatterer.  The  integral  on  the 
right  side  of  equation  (3.2)  represents  the  influence  of 
currents  at  other  surface  points  on  the  current  at  (?,t). 
Note  that  the  influence  of  other  currents  or  the  current  at 
(r,t)  is  delayed  by  R.  This  allows  the  surface  current 
density  equation  to  be  solved  by  an  iterative  numerical 
procedure  in  the  time-domain,  rather  than  the  familiar 
frequency-domain  matrix  inversion  process. 

Once   the  surface   currents   have   been  found,    the 
far-scattered  field  can  be  calculated  by  using: 
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roHSc?'T)  =TF     aT    ;     $      C?',t)    x  a   }'ds«  (3.3) 


where:    r0  =    distance   to   the   far-field    observer 

£r   =    unit   vector   from   the   integration    point 

to  the   far-field   observer 

With     the        substitution      of      the        surface     current 

expression   into   the   above    equation,      and  the   assumption   that 

the  incident    field   is    a  ramp,    the   result   for    the   baclcscatter 

direction   is: 


s    tJCRU'TJ    x   ar}    ds'  t3-4) 


where:    r   HA      =  backscatter      ramp   response     of      the 

target 

s(t  )  =  physical  optics  silhouette  area  of 

the  target 

r      ■   distance   of   the      far-field  observer 
o 


from   the  origin 


t     =  t  +r 
s      © 


H1 
^     =  IH1 


Jc    a  Jc ,   correction  currents  resulting 
from  the  incident  ramp  waveform. 
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This  result  provides  the  exact  relationship  between 
the  target  response  and  target  geometry.  A  particular 
unknown  sample  value  at  a  given  point  in  space  and  time  is 
determined  by  the  exciting  field  at  that  same  space-time 
point  and  by  the  scattered  fields  from  earlier,  more  distant 
locations.  Note  that  the  interactions  between  unknown 
samples  are  displaced  in  time  by  an  amount  equal  to  the  time 
required  for  a  field  to  propagate  between  the  samples  at  the 
speed  of  light.  The  unknown  samples  can  be  solved  at  any 
time  step,  provided  all  sample  values  at  earlier  times  are 
already  known.  Thus,  by  determining  the  unknown  samples, 
i.e.,  the  correction  currents  Jl,  and  adding  the 
contribution  of  the  target  area  function,  s(t  ),  the  target 
ramp  response  can  be  utilized  to  recover  information 
involving  the  unknown  dimensions. 
2-   Inverse  Scattering  Solution 

The  targets  employed  in  the  TDRL  are  all 
rotationally  symmetric  scatterers,  similar  to  the  one 
diagramed  in  Figure  3.1.  The  general  constraints  applicable 
to  the  problem  are  that  the  scatterer  is  symmetric  about  the 
z-axis,  the  incident  field  is  axially  incident,  and  the 
far-field  is  comouted  in  the  backscatter  direction. 
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Figure  3.1. 


Geometry  of  aotationally  Symmetric  Scattering 
Problem 
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The  contour  of  a  rotationally  symmetric  target  is 
completely  described  by  the  radius  vector  p  (z)  which 
varies  as  a  function  of  z.  The  projected  area  function  in 
the  plane  orthogonal  to  the  direction  of  the  propagation  of 
the  field  (z-axis)  can  be  expressed  simply  as: 

S  =  irp2(z) 

Substituting  the  projected  area  function  into  the 
equation  for  the  backscatter  ramp  response  of  the  target  and 
solving        for  (z)  ,  an     inversion        equation     for        the 

rotationallly  symmetric   case   is   obtained: 

PU)    -      UV^.t)    -   ATF     4      Ue»CUt)    x   ar)    ds']-;H}1/2    (3.3) 

Thus,  p  (z)  is  given  in  terms  of  the  measured  ramp  response 
(a  field  measurement) ,  and  in  terms  of  correction  currents 
at  earlier  times,  which  have  been  previously  computed  or  are 
known  to  be  zero. 

In  the  TDRL,  an  iterative  approach  is  employed  using 
estimates  of  the  entire  target  geometry  to  solve  for  the 
contour  function.   The  procedure  involves  five  steps: 

1.   Neglecting   the  surface    integral   providing   the 

correction  currents,    obtain   the   physical   optics 

estimate  of  O   (z)  . 
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2.  Determine  the  correction  currents  using  the  numerical 
solution  of  equation  (3.2)  for  surface  current 
density. 

3.  Apply  the  correction  currents  to  equation  (3.5)  to 
obtain  a    next   estimate   of     O  (z)  . 

pxCz)    =    [2rQ  H*   (?,t)]1/2 

4.  Compare  the  new  value  Q  (z)  with  p  (z)  to  see  if 
the  change  is  less  than  some  small  number,  an  error 
factor.  The  normalized  mean  square  error  is  defined 
as: 


,   CH^Ct)  -  h|   (t))2 
2  _  __2 Rk-1 

K  "       ruSr^.2 


(Hjct)] 


where:     £         =  the  value  to  be  minimized 

IS 

in      order   to      obtain      the 
most  accurate  shape 

description 

c 

H     (t)  =  the  kth  estimation  of  the 

backscattered  field. 
Thus,  the  difference  between  the  estimated  and  actual 
backscattered   fields   are  directly   linked   to   the 
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difference  between  successive  contour  functions,  and, 
therefore,  successive  shapes  of  the  body. 
5.   Repeat  the  iteration  process  until  a  minimum  error  is 
obtained. 
The  above  procedure  was  used  successfully  by  Morag 
[Ref.  15]    in  the  development   of   an  algorithm  involving 
theoretical  analysis   of  targets.    This  algorithm   has  been 
integrated   with   the   overall  TDRL   algorithm   for  inverse 
scattering. 

B.   BACKSCATTERED  WA7EF0RB  HEASURESEHTS 

The  TDRL  algorithm  is  designed  to  minimize  the  effects 
of  the  impulse  responses  of  the  transmitting  and  receiving 
antennas.  Figure  3.2  diagrams  the  parameters  involved. 

As  can  be  seen,  there  are  two  measurement  paths  of 
interest.  The  first  is  the  Direct  Path  between  the 
transmitting  and  receiving  antenna.  The  only  factors 
affecting  the  transmitted  waveform  are  the  intervening 
distance  and  the  impulse  responses  of  the  receiving  and 
transmitting  antennas  to  the  source  pulse.  The  second  path 
is  the  Augmented  Path,  from  the  transmitting  antenna  to  the 
target  and  back  to  receiving  antenna.  Factors  effecting  the 
transmitted  pulse  in   this  path  include  the   antenna  impulse 
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Figure  3.2.   Signal  Paraaeters  for  TDRL  aeasurenents 
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responses,  the  impulse  response  of  the  target  and  the 
corresponding  intervening  distances. 

Three  waveforms  are  measured:  the  Direct  Waveform,  the 
Incident  Waveform  and  the   Augmented  Waveform.    During  the 

I  2 

measurements,  h  (t)  and  h  (t)  ,  responses  due  to  the  delays 
of  the  intervening  distances,  are  windowed  out  by  using  the 
time-range  gating  feature  of  the  DPO. 

The  Direct  Waveform  is  measured  first.  The  transmitting 
and  receiving  antennas  are  boresighted.  No  target  is  on  the 
imaging  plane.  The  measured  wave  is  then  the  convolved 
responses  of  the  source,  and  transmitting  and  receiving 
antenna  impulse  responses.  3y  using  the  properties  of 
convolution,  the  frequency  domain  representation  converts 
the  convolution  process  into  a  product.  The  frequency 
domain  representation  cf  the  Direct  Waveform  is: 

A(ejw)  =  S(eJw)Ht(eJw)Hr(eJ^ 
where:  A(eiu>)   =  Direct  Waveform 
S  (e*u  )   =  source  pulse 


n%*») 


=  transmitting  antenna  response 


H(eicJ)   =  receiving  antenna  response 
The  effects  of  the  intervening  distance  have   been  time 
gated  out.    Figure  3. 3  is   a  representative  Direct  Waveform 
measured  in  the  TDRL. 
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TARGET:    RUN   DATE 

i      8-11-8! 


OIST        ANT   TGT 
1.2711      1         1 
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CYLINDER  ACQUISITION  AVERAGE 
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Figure  3.3.   Typical  Direct  Wavefora — Cylinder 
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The   Incident      Waveform   is      measured  next.  The   imaging 

plane  is  as  set-up  for  the  Direct  Waveform,  but  the 
time-gating  is  further  delayed  by  JAt  to  eliminate  the 
effects  of  the  distance  intervening  between  the  receiving 
antenna     and      the     target      site.  The     resulting      measured 

waveform   is    the   Direct    Waveform    delayed  by   <j£t: 

B(eja))    =   A(ejuJ)    e^^dLt) 

where:  3(e3u>)  =  Incident  Waveform 
Figure  3.4   provides  a   representative  time-domain  Incident 
Waveform  measured  in  the  TDRL. 

The  final  waveform  measured  is  the  Augmented  Waveform. 
It  is  a  convolution  of  the  Direct  Waveform  with  the  impulse 
response  of  the  target  and  summed  with  the  Incident 
Waveform.  For  this  measurement,  the  transmitting  and 
receiving  antenna  remain  boresighted.  A  target  is  placed  on 
the  image  plane.  The  measured  frequency  domain 
representation  is: 

where:  C(ew)   =  Augmented  Waveform 
H  (eiw)  =  Target  Response 
Figure  3.5   provides  a  representative   time-domain  Augmented 
Waveform   measured   in   the   TDRL    for   a    axisymmetric 
half-cylinder. 
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Figure  3-4.      Typical  Incident    SJavefora — Cylinder 
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TARGET:  RUN  DATE 
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Figure  3-5.   Typical  Augmented  Waveform — Cylinder 
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As  can  be  seen  from  equation  (3.6) ,  the  target  response 
is  not  independant  of  the  effects  of  the  impulse  responses 
of  the  antennas.  To  eliminate  this  dependance,  the 
backscattered  waveform  is  determined  and  then  the  remaining 
antenna  terms  are  theoretically  eliminated  by  a  division 
process. 

The  Backscattered  Waveform  is  the  difference  of  the 
Augmented  and  the  Incident  Waveforms: 

D(:>)  =  C(e^)  -  B(eJu) 
=  A(ej(J)  HS  (eja)) 

Figure  3.6  is  the  time-domain  3ackscattered  waveform 
resulting  from  the  subtraction  of  the  waveforms  in  Figures 
3.5  and  3.4  respectively. 

3y  next  dividing  the  Backscattered  Wave  by  the  Direct 
Waver  all  antenna  influences  are  eliminated: 

hVu)  =  d^ 

A(e^) 

o 

The  target  impulse  response,  h  (t)  ,  may  now  be  determined  by 

deconvolution  of  the  target  impulse  frequency  domain  values. 
Figure  3.7  is  a  representative  time-domain  target  impulse 
response  for  a  cylinder.  Noise  errors  have  been  minimized 
by  an  adaptive  filtering  method  to  be  described  shortly. 
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TARGET:  RUN  DATE     DIST   ANT  TGT 
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Figure   3.7. 
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C.   NOISE  CONSIDERATONS 

As  noted  in  the  previous  section,  deconvolution  of  the 
target  backscattered  impulse  response,  less  the  impulse 
responses  of  the  antennas,  is  a  central  operation  to  the 
effective  use  of  the  IDRL.  The  reliability  of  the  acquired 
signal  to  be  a  true  response  of  a  target  to  an  incident 
field  bears  directly  on  the  accuracy  of  the  final  results  as 
determined  by  solution  of  the  inverse  integral  equations  for 
time-domain  analysis.  Any  srrors  in  the  data  representing 
the  Incident  Waveform  and  Augmented  Waveform  degrade  the 
cancellation  of  the  division  process  and  introduce  errors  in 
the  backscattered  target  response.  Also  present  at  the 
input  to  the  DPO  is  noise  due  to  thermal  electron  excitation 
in  the  transmission  lines  and  antenna  surfaces  as  well  a 
that  due  to  external  atmospharic,  cosmic  and  man-made  noise 
sources  being  received  by  the  horn  antenna.  When  the 
backscattered  impulse  response  is  transformed  to  the 
time-domain,  high  frequency  noise  in  most  cases  forms  the 
dominant  feature  of  the  waveform.  An  example  of  a  typical 
impulse  response  of  a  cylindar  is  given  in  Figure  3.8. 
Noise  errors  have  not  been  removed.  Figure  3.9  is  the  ramp 
response  for  the  same  target. 
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TARGET:  RUN  DATE     DIST   ANT  TGT  REMARKS 
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Figure  3.8.   Impulse  Response — Noisy  Signal 
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TARGET:  RUN  DATE     DIST   ANT  TGT   REMARKS 
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Figure   3.9.      Ramp  Response — Noisy   Signal 
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In  the  TDRL,  typical  sources  of  errors  in  the  Incident 
and  Augmented  Waveforms  are  sampling  errors  and  noise  in 
signal  acquisition,  leakage  and  aliasing  errors  in  the 
Fourier  Transf omations ,  and  rounding  errors  in  the  computer. 
The  error  of  primary  concern  is  that  involving  noise  in 
signal  acquisition  and  in  particular  the  enhanced 
high-frequency  noise  resulting  from  decon volution  of  the 
antenna  responses  to  obtain  the  impulse  and  ramp  responses 
of  the  target.  This  noise  has  been  directly  traced  to  the 
structure  of  the  data  within  the  input  array,  as  will  be 
noted  in  Chapter  IV. 

Two  methods  are  employed  to  reduce  noise  to  acceptable 
levels:  acquisition  averaging  and  an  optimal  compensation 
technique.  [ Hef .  16] 

1 .   Acquisition  Averaging 

Acquisition  averaging  is  aimed  at  reducing  signal 
acquisition  errors  and  noise  by  averaging  several 
acquisitions  taken  within  a  short  period  of  time.  During 
the  process,  a  more  serious  error  may  be  generated,  that  due 
to  signal  drift  with  time.  Signal  averaging  has  been  most 
useful  in  ensuring  a  good  target  ramp  response.  It  has  had 
noticeable  effect  on  making  the  information  content  of  the 
impulse  response  more  available.     Figure  3.10  demonstrates 
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an  impulse  waveform  with  210  acquisitions  for  the  same 
target  as  for  Figure  3.8.  Theoretically,  the  RMS  noise 
level  should  be  reduced  by  the  square  root  of  the  number  of 
acquisitions,  N.  N  =  10  in  Figure  3.3,  and  N  =  210  in 
Figure  3.10.  The  noise  reduction  should  be  about  4.5  times. 
Note  that  the  noise  level  in  Figure  3.10  is  indeed  about  one 
fourth  that  in  3.8. 

Figure   3.11   is  the  corresponding   ramp  response. 
Signal  acquisition  only  is  employed. 
2-   QE.liM.1  Compensation  Technique 

The  optimal  compensation  technique  involves  the 
design  of  a  compensator  ( deconvolution)  function  operating 
on  the  convolution  output  of  the  Incident  and  3ackscattered 
Waveforms,  to  produce  a  relatively  noise-free  estimated 
Impulse  Response  for  the  target.  The  compensator  is  applied 
to  the  convolution  product  of  the  target  impulse  response 
and  the  Direct  Waveform  to  produce  the  required 
deconvolution  result.  The  design  as  used  in  the  TDRL 
involves  an  iteration  process  on  a  single  variable  with  a 
"man-in-the-loop"  to  determine  optimal  values  of  the 
variable. 

Figure  3.12  diagrams  the  design  of  a  frequency 
domain  optimal  compensator.    The   compensation  principal  in 
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TARGET:  RUN  DATE     DIST   ANT  TGT   REMARKS 

1   8-11-81   1.27M   1    1    CYLINDER  ACQUISITION  AVERAGE 
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Figure  3.10. 


Iapulse    Response:    Mo  Optimization — Acquistion 
Averaged 
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Figure  3.11. 


Raap  Response:  Ho  Optimization — Acquistion  Av- 
eraged 
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Figure  3.12-      Frequency  Domain  Optimal  Compensator  Design 

deconvoiution   is   to   design   a   transfer    function    Q  (e     )       to    be 
applied     to   the      backscatter   response      Dfe'5'*')         to   yield      an 

estimate  of    the   target    impulse    response   H.  (e  ua)  .      That    is: 

c 


HS(e^)  =  D(eJw)Q(eJM) 


(3.7a) 


=  ACeJW)HS(eJu)Q(e>) 


(3.7b) 


The  main  design  criteria  is  to  minimize  the  energy 
in  the  noise  terms  that  introduce  error  in  the  estimate.  The 
error  energy,  E   is  given  by: 


Ee  "  4   I  He  CeJW)  -  H*  (e^) 


dw 


(3.8) 


where:  X\=    frequency  band  of  interest. 


S  .  *o 


When  the   error  energy  is  minimized,    H^(e  )   will 

equal  H  (eio>)   and  the  impulse   responses  will  be  identical, 

except  for  the  addition  of  highly  enhanced  random  noise.   In 

this  case,  the  compensator  transfer  function  is: 
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3ut,  this  will  result  in  an  infinite,  unbounded  value  for 
Q  (e  )  as  A  (e  )  goes  to  zero,  yielding  i  noisy  form  for 
H  (e^"9).  This  noise  can  be  limited  by  constraining  Q(e  ) 
to  be  a  bounded  function.  Noting  that  the  impulse  response 
of  the  target  is  bounded  by  physical  constraints,  it  is 
possible  to  keep  the  energy  due  to  the  convolution  of  the 
impulse  response  and  the  designed  transfer  function  bounded, 


i.e. : 


Ec  =  fo    I  H^ju)  Q(e>)  !2  d"  C3-10) 


The   problem   now     is  to  minimize   S        while   keeping   E 
finite.        The   total   energy   in      the   backscattered    wave   can   be 
defined  as: 


E   =   E      +XE      ,      X    >_  0  (3.11a) 


=   ?     [|    HS    (e>)    -    HS    (e^)     |2    ♦    X    |  H(eJu)Q(ej u)  |  2]  dw    (3.11b) 


where:    \  -   optimization    parameter. 
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Thus,  by  minimizing  the  total  error  energy  through 
the  selection  of  the  optimal  value  of  X  *  ^he  error  energy 
in  the  estimate  is  also  minimized.  The  result  yields  an 
optimum  compensator  Q  +  (eJU5)  for  the  unknown  target  impulse 
response  h  (t)  .  The  result  is  that  the  whole  system  output 
is  low-pass  filtered  with  minimal  effect  on  waveform  shape. 

Adhering  to  the  design  constraints  listed  above,  it 
can  be  shown  that  the  form  of  the  compensating  transfer 
function  Q  (e  )  for  a  givan  Direct  Waveform  A(eioJ)  is: 
[Ref.  17] 


Q0>)  =    A*{eJ") C3J2) 

[|A(eJw)|2*A] 


where:  A* (e   )  =  complex  conjugate 
From  equation  3.7,   the  compensation  deconvolution 
process  yields: 


S  .  Dfe>U*(eJU))  (3.13) 

e    [  A(e>)|Z  +  X] 


Figure  3.13  shows  the  affect  that  optimization  has 
on  the  impulse  response  with  no  sgnal  acquisition  averaging. 
Optimization  is   0.5.    The  target  impulse   response,   while 
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target:  run  date    dist   ant  tgt 

1   3-11-31   1.27N   1    1 


REMARKS 

TYPICAL  CYLINDER  RESPONSES 


MAXIMUN  PEAK  UALUE 9. 0011  nU 

MINIMUM  PEAK  UALUE -0.9010  nU 

RMS  UALUE 0.  0004  nU 

MEAN  UALUE 0.9090  fiU 

NUMBER  OF  WAUEFORMS  AUERACED  =  1      OPTIMIZATION  UALUE  »  9.5 
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8.59  1.00  1.58  2.88  2.58  3.98  3.58  4.88  4.58  5.80   nS 


Figure  3,13, 


Impulse  Response:  .5  Optimization- 
No  Acquisition  Averaging 
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still  noisy,  is  clearly  discernable.  Figure  3.13  should  be 
compared  with  Figure  3.11.  Figure  3.14  shows  the  effect  of 
optimization  increased  to  10.  The  noise  is  almost  totally 
eliminated.  The  trade-off  is  that  the  optimization  has 
reduced  the  amplitude  of  the  response  and  spread  it  in 
frequency.  This  is  a  general  observation.  However,  it 
appears  that  the  degree  of  effect  on  the  waveform  is  a 
function  of  the  ideal  optimization  factor.  Chapter  V  will 
explore  this  concept  more  fully.  However,  it  should  be 
noted  that  if  the  optimization  factor  is  large  enough  to 
reduce  the  large  scale  features  of  the  target  then  it  is  too 
large  and  should  be  reduced — this  is  the  "man  in  the  loop's" 
job. 

The  combination  of  waveform  averaging  and 
optimization  on  the  impulse  response  is  demonstrated  in 
Figures  3.15  and  3.15.  With  a  relatively  small  amount  of 
optimization,  the  target  impulse  response  is  readily 
apparent.  See  Figure  3.16.  Signal  acquisition  averaging 
also  yields  an  important  advantage.  Since  noise  is  reduced 
directly  by  the  square  root  of  the  number  of  acquisitions, 
and  because  no  filtering  of  any  sort  is  employed  for 
averaging  alone,  the  detailed  features  of  the  actual  signal 
are  not  eliminated  as  with  filtering.   The  smaller  amount  of 
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optimization  needed  results  in  minimum  effect  by  the 
compensator  on  the  shape  of  the  output  waveform.  The 
trade-off  is  the  increased  liklihood  of  error  introduced  due 
to  signal  time  drift  during  the  acquistion  averaging  period. 
Figures  3.17  to  3.21  show  the  effects  of  signal 
averaging  and  optimization  on  the  Ramp  Responses.  Each 
figure  is  presented  in  the  same  order  as  were  the  Impulse 
Responses.  Signal  averaging  is  very  beneficial  in  obtaining 
an  accurate  ramp  response.  Optimization  has  a  greater 
tendency  to  effect  the  wave  shape  of  the  ramp  response  than 
the  impulse  response  of  the  target.  But  less  optimization 
is  required  to  reduce  the  error  noise.  It  is  generally 
possible  to  use  the  ramp  response  without  optimization.  The 
main  effect  of  optimization  is  to  spread  the  frequency 
content  of  the  ramp  response.   Amplitude  effects  are  minor. 
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TARGET:  RUN  DATE 

1   8-11-81 


DIST   ANT  TGT  REMARKS 

1.27M   1   1    TYPICAL  CYLINDER  RESPONSES 


MAXIMUM  PEAK  UALUE 9.8094  nU 

MINIMUM  PEAK  UALUE -9.9994  nU 

RMS  UALUE 9.9892  nU 


MEAN  UALUE , 

NUMBER  OF  WAUEFORMS  AUERAGED  =  1 


8.8899  nU 

OPTIMIZATION  UALUE 


18 


IMPULSE  RESPONSE 


-8.0002 


8.30  1.80  1.50  2.00  2.50  3.00  3.50  4.00  4.50  5.00   nS 


Figure  3.1^. 


Impulse  Response:  10  Optimization — No  Acquisi- 
tion Averaging 
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TARGET:  RUN  DATE     DIST   ANT  TGT   REMARKS 

1   8-1 1-81   1.27M   I    1    CYLINDER  ACQUISITION  AVERAGE 

MAXIMUM  PEAK  VALUE 3.0994  nU 

MINIMUM  PEAK  UALUE -9.0891  nU 

RMS  UALUE 9.9933  nV 

MEAN  UALUE 9.9099  nU 

NUMBER  OF  WAUEFORMS  AUERAGED  =  21      OPTIMIZATION  UALUE  *  9 


«U 


9.99S9 


9.9989 


-8.8858 


IMPULSE  RESPONSE 


8.58  1.89  1.58  2.88  2.58  3.88  3.58  4.86  4.58  5.98   nS 


Figure  3.15. 


Impulse  Response:  No  Optimization — 210  Acqui- 
sitions Averaged 
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TARGET:  RUN  DATE 

1   8-11-81 


DIST   ANT  TCT 
1.27M   1    1 


REMARKS 

CYLINDER  ACQUISITION  AUERAGE 


MAXIMUM  PEAK  UALUE 0.8019  nU 

MINIMUM  PEAK  UALUE -9.0009  «U 

RMS  UALUE 9.0003  mU 


MEAN  UALUE 

NUMBER  OF  WAUEFORMS  AUERACED  =  21 


0.9090  nU 

OPTIMIZATION  UALUE  «  9.5 


flU 

9.9819 

9.9895 

9.9999 

-9. 9995 


IMPULSE  RESPONSE 


- 

\ 

j!\u 

\ 

\A    ' 

Va 

J1 

\ 

X 

\A 

At^> 

f\ 

\ 

*\ 

/•V1 

"  V 

A? 
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r 

V 

8.58    1.98    1.58    2.90    2.58    3.88    3.50    4.88    4.58    5.88       nS 


Figure  3. 16. 


Impulse    Response:    .5   Optimization — 210   Acquis- 
tions   Averaged 
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TARGET:  RUN  DATE     OIST   ANT  TCT  REMARKS 

t   8-11-91   1.2711   1  1     TYPICAL  CYLINDER  RESPONSES 

MAXIMUM  PEAK  UALUE +53.65  nU-M 

MINIMUM  PEAK  UALUE -700.00  mU-M 

RMS  UALUE +373.03  nU-M 

MEAN  UALUE -283.42  mU-M 

NUMBER  OF  WAUEFORMS  AUERAGED  »  1      OPTIMIZATION  UALUE  «  0 


nU-M 


RAMP  RESPONSE  ~  TIME 


-298 


-480 


-600 


- 

- 

/ 

Wvv 

<v/" 

/* 

- 

- 

J 

i 

- 

tym 

..     ... 

9.33  1.98  1.30  2.00  2.30  3.00  3.50  4.00  4.30  5.08   nS 


Pigure  3,17. 


Ramp  Response:  No  Optimization — No  Acquisition 
Averaging 
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TARGET:  RUN  DATE 

1   8-11-81 


DIST   ANT  TGT  REMARKS 

1.27M   1    1    TYPICAL  CYLINDER  RESPONSES 


MAXIMUM  PEAK  UALUE +47.77  nU-M 

MINIMUM  PEAK  UALUE -647.48  nU-M 

RMS  UALUE t-367.  1 1  nU-M 

MEAN  UALUE -279.  80  nU-M 

NUMBER  OF  WAUEFORMS  AUERAGED  =  1      OPTIMIZATION  UALUE  -  9.3 


nU-M 

RAMP   RESPONSE   - 

•   TIME 

8 

- 

-208 

-488 

- 

-680 

• 

8.38  1.00  1.58  2.00  2.30  3.80  3.30  4.00  4.30  5.00   nS 


Figure  3.18. 


Ranp  Response:  .5  Optimization — No  Acquisition 
Averaging 
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TARGET:  RUN  DATE     DIST   ANT  TGT 
1   8-11-81   1.27H   1    1 


REMARKS 

CYLINDER  ACQUISITION  AUERACE 


MAXIMUM  PEAK  UALUE +69.81  nU-M 

MINIMUM  PEAK  UALUE -481.43  mU-M 

RMS  UALUE +231 .  1 1  mU-M 

MEAN  UALUE -178.58  «U-M 

NUMBER  OF  UAUEFGRHS  AUERAGED  »  1      OPTIMIZATION  UALUE 


18 


nU-M 

58 

8 

-58 

-188 

-158 

-288 

-256 

-388 

-358 

-488 


RAMP   RESPONSE    —    TIME 


^ 

s 

y 

■ 

- 

- 

■ 

» 

- 

L  / 

- 

- 

8.58    1.88    1.38    2.88    2.58    3.88    3.38    4.88    4.58    5.88       nS 


Figure   3.19. 


Raap  Response:  10  Optimization — No  Acquisition 
Averaging 
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TARGET:  RUN  DATE     OIST   ANT  TGT  REMARKS 

I   8-11-81   1.27M   1  1    CYLINDER  ACQUISITION  AVERAGE 

MAXIMUM  PEAK  VALUE +73.58  nU-M 

MINIMUM  PEAK  UALUE -4S7.95  nU-M 

RMS  UALUE +239.  59  nU-M 

MEAN  UALUE -223.  75  nU-M 

NUMBER  OF  WAVEFORMS  AVERAGED  =  21      OPTIMIZATION  UALUE  =  3 


nV-M 


RAMP  RESPONSE  --  TIME 


-238 


-483 


- 

- 

■ 

- 

■1 

JK**^ 

3.58  1.80  1.58  2.00  2.50  3.80  3.50  4.00  4.50  5.00   nS 


Figure    3.20.      Ramp  R?sP2^?i 
y  tions   Averaged 


No  optimization— 210  Acquisi- 
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TARGET:  RUN  DATE 

1   3-11-81 


DIST   ANT  TCT  REMARKS 

1,2711   1   1    CYLINDER  ACQUISITION  AVERAGE 


MAXIMUM  PEAK  VALUE +72.25  nU-M 

MINIMUM  PEAK  VALUE. -465.94  nV-M 

RMS  UALUE +288.81  nV-M 

MEAN  UALUE -222.20  mV-M 

NUMBER  OF  WAVEFORMS  AUERAGED  *  21      OPTIMIZATION  UALUE  »  8.5 


wV-H 


RAMP  RESPONSE  —  TIME 


-288 


-408 


- 

- 

• 

^y 

0.53  1.88  1.58  2.00  2.50  3.80  3.50  4.08  4.50  5.00   nS 


Figure  3.21. 


Raap  Response:  .5  Optimization — 210  Acquisi- 
tions Averaged 
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II-   QEI1ATING  SYSTEM  COMPUTER  IMPLEMENTATION 

A.  INTBODOCTION 

k  set  of  software  programs  called  the  "Operating  System" 
controls  the  input,  processing  and  output  of  information 
obtained  in  the  TDRL.  It  is  the  quality  and  throughness  of 
the  Operating  System  that  directly  determines  the  quality  of 
the  results  obtained.  In  this  chapter,  there  will  be  a 
discussion  in  detail  of  the  implementation  of  the  Operating 
System.  It  will  begin  by  discussing  the  physical  structure 
of  the  the  host  processor.  A  brief  overview  of  the 
Operating  System  will  follow,  defining  the  basic 
sub-programs  that  make-up  the  total  algorithm.  Finally,  a 
detailed  description  of  eaoh  program  will  be  provided, 
emphasizing  the  method  of  implementation  of  the  numerical 
solutions,  where  applicable,  and  giving  details  of  a 
procedure  to  run  tha  individual  programs. 

B.  THE  TEKTRONIX  4052  GRAPHIC  COMPUTING  SYSTEM 

The  host   processor  used  in   the  TDRL  is   the  Telctronics 

4052  Graphic  Computing   System.    It  is  a   versatile,   stand 

alone   microprocessor.     As  such,    its   capabilities   and 

limitations  are   important  driving  functions  that   shape  the 

structure  of  the  Operating  System. 
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The  4052  has  three  memories.   They  are: 

1.  Read  Only  Hemory  (ROM) — a  permanent  memory  that 
contains  the  System's  intelligence. 

2.  Random  Access  Memory  (RAM) — working  memory  of  the 
processor's  CPU . 

3.  Line  Buffer — small  temporary  memory  that  allows  data 
to  be  written  and  edited  on  a  display  screen  prior  to 
releasing  it  to  the  RAM.  An  important  editing 
feature. 

Built-in  peripherals  for  the  Graphic  System  (GS)  are  a 
smart  keyboard  (primary  input) ,  a  CRT  visual  display  screen 
(primary  output)  and  a  magnetic  tape  unit  (mass  storage). 
Three  extenal  peripheral  devices  extend  the  versatility  of 
the  system.  They  are  a  Tektronix  4631  Hard  Copy  Unit  (makes 
paper  copies  of  display  information) ,  a  4924  Digital 
Cartridge  Tape  Drive  (additional  mass  storage)  and  a  Digital 
Processing  Oscilloscope  (data  acquisition).  A  block  diagram 
for  the  system  is  shown  in  Figure  4.1.  [Ref.  18] 

The  Processor  is  the  main  computing  device.  It 
corresponds  to  the  Central  Processing  Unit  (CPU)  found  in 
larger  systems.  It  maintains  the  "firmware"  that  allows  the 
GS  to  direct  system  operations,  decode  instructions  and 
perform  arithmetic  and  logic  operations. 
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pBMBW, 


«»<10  COPT 
UNIT 

I  ! 

Figure   U.  1 .      Block  Diagram   of  TDEL  Graphic   System 

The  General  Purpose  Interface  provides  the  means  for  the 
system  to  communicate  information  to  or  receive  data  from 
any  external  peripheral  that  is  compatible  with  IEEE 
Standard     U88-1975.  This      interface      is      a      bit-parallel, 

byte-serial  system  capable  of  transfering  ASCII  code  or 
machine  dependant  binary  cods  up  to  8-bits  wide  at  a  maximum 
rate  of   250- Ic  bytes/second. 


96 


A  magnetic  tape  data  cartridge  is  used  in  the  internal 
magnetic  and  peripheral  tape  drive  units  to  provide  up  to 
300-k  bytes  of  additional  mass  storage  per  tape.  Programs 
and  data  may  be  stared  on  tape  in  specifically  defined 
"files".  Piles  may  range  in  size  from  768-bytes  up  to  the 
limit  of  the  magnetic  tape  storage  capability.  Files  may  be 
either  ASCII  or  binary  in  format.  All  files  in  the 
Operating  System  are  stored  in  Binary  format  to  take 
advantage  of  a  peripheral  "firmware"  package  called  the 
binary  loader.  3inary  data  transfers  between  the  processor 
compiler  and  magnetic  tapes  are  faster  than  when  stored  in 
ASCII.  This  is  because  the  processor  uses  binary  format  for 
internal  processing  and  the  conversion  from  ASCII  to  binary 
format  is  eliminated  when  data  is  stored  binary. 
Additionally,  a  stored  binary  program  can  be  transfered  to 
the  processor  memory  from  a  specified  peripheral  device 
without  disturbing  variables  and  associated  values 
previously  defined.  This  feature  is  not  available  to  ASCII 
programs. 

The  GS  permanent  ROM  has  been  programmed  to  respond  to 
3ASIC  (Beginner's  All- Purpose  Symbolic  Instruction  Code),  a 
high-level  programming  language  first  developed  at  Dartmouth 
College.     The  Operating   System   has  been   correspondingly 
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written  in  this  language.  However,  BASIC  as  used  by  the  GS 
has  been  extended  to  fit  the  language  to  the  specialized 
capabilities  of  Graphic  System.  These  extensions  are  in  the 
areas  of  graphics,  file  system  access,  unified  handling  of 
input/output  operations,  matrices,  character  string 
manipulation,  high-level  language  interrupt  handling,  and 
operating  system  facilities.  A  key  difference  between 
extended  GS  BASIC  and  most  other  BASIC  languages  is  that 
most  keywords  and  their  parameters  can  be  evaluated 
independantly  of  program  control.  The  result  is  a  rich, 
versatile,  yet  simple  to  use  programming  tool. 

Three  peripheral  ROM  firmware  packages  have  been  added 
to  the  TDRL  system  to  tailor  it  to  specific  needs  of  the 
laboratory.  They  are  the  Signal  Processing  ROM  Pack  No.  1, 
Signal  Processing  ROM  Pack  No.  2  (FFT)  ,  and  the  EDITOR.  A 
description  of  each  is  provided  in  Appendix  A. 

C.   OVERVIEW  OP  THE  TDRL  OPERATING  SYSTEM 

The  TDRL  Operating  System  has  been  specifically  designed 
to  be  higly  operator-processor  interactive.  This  provides 
the  greatest  degree  of  flexibility  in  the  acquisition, 
processing  and  evaluation  of  electromagnetic  time  domain 
transients  for  target  imaging. 


98 


The  total  Operating  System  occupies  83818-bytes  of  mass 
storage  on  TDRL  Library  Tape  1  in  three  files.  It  requires 
an  additional  201628-bytes  of  storage  in  62-files  on  TDRL 
Library  Tape  2  (Data  Storage) .  The  Operating  System  is 
logically  divided  into  three-individual  programs: 

1.  INPUT  —  for  the   acquisition  and  averaging   of  signals 
of  interest. 

2.  MATH — performs  the  significant  ma-chmatical  processing 
of  the  acquired  signals  for  inverse  scattering. 

3.  GRAPH — provides  a   visual  display  to  the   operator  of 
the  results  of  the  processing. 

TABLE  III 
Mass  Storage  TDRL  Programs 


PROG  HAM 


FILE 


T~T 


BYTES 


INPUT 


24552 


12 


HATH 


14256 


GRAPH 


1 


426  24 


22 


Mass  storage  requirements  of  the  Operating  System 
programs  are  given  in  TABLE  III.  The  storage  requirements 
listed  in  TABLE  III  are  for  the  programs  only  and  do  not 
include  requirements  for  the  necessary  definition  of  arrays, 
strings  or  numeric  constants  and  variables. 
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All  programs  are  designed  to  be  as  nearly  independant  of 
each  other  as  possible,  after  an  initial  run,  to  allow  for 
maximum  flexibility  and  utility  in  data  usage  and 
evaluation.  The  most  dependant  program  is  MATH  which  cannot 
be  used  without  inputs  from  either  INPOT  or  GRAPH.  The  most 
independant  program  is  GRAPH  which  is  capable  of  being  run 
with  inputs  from  mass  storage  only. 

A  general  flow  diagram  of  the  operating  system  is  given 
in  Figure  4.2.  Specific  named  subroutines  have  not  been 
shown.  The  next  section  will  describe  in  detail  the 
individual  programs  of  the  Operating  System  and  their 
specific  requirements  and  processing  methods. 

D.   OPEBATING  SYSTEM  PROGRAM  DESCRIPTION 

This  section  will  describe  in  detail  the  individual 
programs  and  their  subroutines  that  make-up  the  Operating 
System  for  the  TDRL.  A  general  description  of  the  purpose 
of  each  program  will  be  provided  first.  Sources  of  inputs 
and  outputs  will  be  noted.  A  general  flow  diagram  for  the 
specific  program  will  be  provided  showing  how  all 
subroutines  interact.  This  general  overview  will  then  be 
followed  by  a  more  detailed  description  of  each  subroutine 
in   the  program.    The  purpose   of  the  subroutine  will   be 
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noted.  Sources  of  input  and  output  will  be  defined.  A 
discussion  of  applicable  algorithms  and  the  flow  of 
information  will  be  also  provided.  Examples  of  visual 
displays  applicable,  potential  problems  and  correction 
methods,  will  be  discussed  when  warranted. 

As  noted,  the  following  commentary  will  make  frequent 
referance  to  program  generated  messages.  Examples  of  these 
messages  will  not  be  given  in  the  body  of  this  text,  but  can 
be  found  in  Appendix  F,  a  ganeral  run  of  the  Operating 
System . 

1.   Program  INPUT 

INPOT  provides  the  means  to  initialize  and  drive  the 
total      Operating   System.  This      sub-program      is   stored      on 

Library  Tape  1,  File  2.  Memory  space  required  for  its 
implementation  in  the  processor  CPU  is  24552-bytes.  With 
all  arrays  and  constants  defined,  a  total  of  42469-bytes  of 
memory  space  may  be  required.  A  listing  of  INPUT  can  be 
reviewed   in    Appendix   C. 

INPUT  is  the  front-end  program  for  the  Operating 
System  and  as  such,  will  generally  be  the  first  program  to 
be  initiated  for  TDRL  measurements.  Program  start-up  is 
begun  by  placing  Tape  1  into  the  4052  built-in  magnetic  tape 
drive,      and    placing      Data    Storage  Tape   2   in      the   4924    Remote 
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Tape  Drive.  INPUT  is  then  normally  placed  into  the  CPU's 
working  memory  under  operator  control  by  the  issue  of  the 
following  commands  from  the  keyboard: 

FIND  2 
CALL  "BOLD" 
Program  execution  is  next  begun   by  typing  the  command  "RUN" 
from  the  keyboard. 

If  this  is  the  first  attempt  to  run  the  Operating 
System  in  conjunction  with  the  Digital  Processing 
Oscilloscope  (it  is  assumed  the  DPO  is  "ON"  whenever  INPUT 
is  employed) ,  the  following  message  may  appear  on  the  GS  CRT 
following  the  issuance  of  "FIND  2": 

"Ho  SRQ  on  unit  in  immediate  line  —  message  nuaber  43," 
If  the  program  is  being  run  from  any  line  prior  to  line  140 
(the  usual  case),  this  message  may  be  ignored  and  the 
program  will  automatically  correct  the  error  message  and 
continue  processing.  If  the  program  is  being  run  subsequent 
to  line  140  (the  case  when  returning  to  INPUT  from  some 
other  internal  file) ,  the  following  command  must  be  entered 
from  the  keyboard: 

POLL  N,fl}1 
Following  this,  the  regular  sequence  of  commands  may  then  be 
issued  to  start  the  system. 
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As  noted,  INPUT  may  alternately  be  loaded  and 
initiated  under  program  control  from  GRAPH  (File  4) . 
Execution  will  begin  following  line  180  in  this  case.  This 
return  to  INPUT  occurs  in  normal  Operating  System  (OS) 
operation  after  an  initial  run  has  been  completed  and 
information  is  then  desired  for  an  entirely  new  target. 

The  program  specifically  provides  for  the  means  to 
input  information  concerning  the  transient  electromagnetic 
response  of  a  target.  It  acquires,  averages  and  stores  all 
pertinent  waveforms.  INPUT  also  allows  for  the  easy  storage 
of  specific  antenna  and  target  dimension  parameters  for 
comparison  when  the  target  is  known  exactly.  This  data  is 
permanently  maintained  in  a  library  on  Tape  2  in  Files  5-49. 
INPUT  further  initializes  most  of  the  parameters  used  later 
in  Program  GRAPH  for  the  visual  display  of  the  processed 
results. 

Data  inputs  are  from  the  keyboard  or  internal 
memory.  Outputs  are  to  internal  memory,  mass  storage  and 
File  3  (MATH) .  INPUT  is  highly  operator-processor 
interactive.  Numerous  directives  and  visual  displays  guide 
the  user. 

The  structure  of  the  program  consists  of 
1 1-subroutines.    One  (DRIVER)   is  of  general  application  to 


104 


the  program.  Five  (ACAVG,  INWA7,  DIWA7,  AUGWA7,  SCATWA7) 
are  specifically  related  to  input,  acquisition,  averaging 
and  storage  of  relevant  information  concerning  the 
time- domain  transient  response  of  the  target  in  question. 
The  remaining  five  subroutines  (ANTLOC,  TGTLOC,  NEWPAR, 
ANTIN,  TGIIN)  provide  special  services  for  the  operating 
system.  A  general  flow  diagram  for  input  is  provided  in 
Figure  4.3  A  description  of  each  subroutine  follows, 
a.   DRI7ER 

This  is  the  main  driving  routine  for  INPUT.  It 
consists  of  129  steps  (40%  of  the  program)  from  line  100 
through  line  1260.  It  is  initiated  directly  under  operator 
control  at  line  100,  or  under  program  control  from  GRAPH  in 
line  180.  Data  is  input  from  the  keyboard  or  from  the  mass 
storage  files  located  on  Tape  2.  DRIVER  will  pass  program 
control  directly  to  MATH  when  directed. 

Lines  100-300  initialize  various  variables  and 
dimension  various  arrays.  Line  150  clears  the  General 
Purpose  Interface  3us  (GPI3)  of  the  DPO  service  request. 
Flag  ?7  (indicates  whether  the  acquired  data  is  destined  for 
use  in  Time  Domain  or  Prony  Methods)  is  set  to  0  (Time 
Domain)  in  line  160.  P8,  the  run  number,  is  initialized  to 
zero  in  line  170  and  will  be  incremented  by  one  in  line  210 
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for  every  run  (including  the  first).  Various  string 
constants  are  initialized  in  line  240-250  and  280-290.  Flag 
?9  (indicates  the  need  for  initializing  input  parameters)  is 
set  to  1  (input  parameters  required).  Lines  320-530  provide 
a  menu  of  subroutines  and  functions  which  are  available  for 
Operator  selection.  The  selection  of  a  process  from  the 
menu  is  input  from  the  keyboard  into  systei  variable  P0.  An 
example  of  the  menu  is  given  in  Appendix  F. 

Lines  550-570  direct  the  program  flow  to  the 
appropriate  portion  of  DRIVER.  If  the  operator  has  selected 
one  of  the  first  four  subroutines,  and  this  is  an  initial 
run  of  the  OS,  then  the  time  for  a  full  screen  sweep  of  the 
DPO  must  be  inputted  to  the  program  from  the  keyboard.  The 
sweep  rate  is  put  in  in  seconds  in  T1  in  line  640  and  must 
be  greater  than  or  equal  to  the  time  required  for  the  the 
DPO  beam  to  sweep  completely  across  the  DPO  screen. 

Program  control  is  now  passed  to  line  850-1070 
where  the  following  run  descriptive  parameters  are  input: 

1.  Date 

2.  The   distance  the   target   is   from  the   transmitting 
antenna 

3.  The  antenna  number  selected 

4.  The  target  number  selected 
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5.       Explanatory   remarks 
These    are   all  keyboard   entries.         An   exampls    of    the  messages 
that     query      entry     of      the   above      parameters      is      given      in 
Appendix   F. 

If  the  run  is  not  the  first  for  the  Operating 
System,  a  different  series  of  inputs  are  queried.  In  this 
case,  program  control  is  passed  to  lines  670-840  which 
generate  messages  on  the  GS  asking  the  operator  if  there  is 
a  change  in  the  original  location  of  the  target,  if  a  new 
antenna  is  to  be  used,  or  if  a  defined  target  is  to  be 
changed.  The  operator  inputs  the  alpha  character  "T"  or  "N" 
for  ••yes"  or  "no"  from  the  keyboard  to  reflect  the  current 
situation.       Appendix    F    displays   these    messages. 

Lines  1000-1040  generate  the  input  of  the 
antenna  parameters.  These  parameters  are  drawn  from  files 
stored  on  Data  Tape  2.  Each  antenna  file,  located  in  one  of 
the  files  numbered  5-24,  is  a  63-element  array,  returned  in 
A5.      In   elements    1    through   60   are   located: 

1.  the      run    number      for     this   particular      antenna-target 
combination 

2.  the      antenna   length      compensating      parameter  for      the 
antenna-target    combination 
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3.  the  antenna  amplitude  compensator  for  the 
antenna-target  combination 
ks  can  be  seen,  antenna  parameters  are  actually  three 
elements  within  the  larger  array,  each  three  elements 
identifying  a  unique  antenna-target  combination.  The 
parameters  in  elements  51-63  are  system  antenna  parameters 
used  for  unknown  targets.  They  are  simply  the  sum  of  all 
run,  antenna  length  compensator  updates  and  antenna 
amplitude  compensator  updates.  The  average  of  the  last  two 
elements  is  applied  to  signals  acquired  if  the  target  is  not 
defined  from  the  target  library.  This  is  shown  in  lines 
1024-1025. 

The  antenna  amplitude  compensator  is  a  very 
sensitive  function  of  the  distance  the  target  is  from  the 
transmitter  due  to  EM  propagation  attenuation.  In  line 
1025,  for  unknown  targets,  and  in  line  4  152  for  defined 
targets,  these  losses  are  taken  in  consideration  by  finding 
the  square  of  the  ratio  of  the  targets  true  location  on  the 
plane  with  referance  to  the  transmitting  antenna  and  the 
control  distance  (1.27m). 

After  the  definition  of  variables  is  complete, 
DRIVES  passes  control  to  the  appropriate  function  in  line 
1090.   Line  1100  returns  the  display  to  the  Menu  for  further 
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function  selection  or  to  line  1170  where  the  DRIVER  passes 
control  to  HAIH  and  terminates  any  further  INPUT  processing. 
In  summary,  INPUT  DRIVER  requires  input  of 
selected  program  routines  and  initial  parameters  by  the 
operator  from  the  keyboard,  and  requires  the  OS  to  determine 
whether  the  current  run  is  the  initial  for  the  day  or  a 
subsequent  one.  The  parameters  are  saved  for  later  graph 
requiremts.  Funtion  selection  drives  the  appropriate 
subroutines  and  routines.  A  flow  diagram  for  GRAPH  DRIVER 
is  provided  in  Figure  4.4. 

b.  Subroutine  ANTLOC/TGTL3C 

ANTLOC  and  TGTLOC  are  similiar  subroutines  in 
that  they  locate  the  appropriate  antenna  and  target  storage 
arrays    filed   on   Tape   2.  Antenna    parameters   are   located   in 

Files  5-24.  Target  parameters  are  located  in  files  26-45. 
Lines  4000-4060  contain  ANTLOC.  TGTLOC  is  found  in  lines 
4100-4160.  TGTLOC  has  the  additional  function  of  defining 
the  system  antenna  parameters  for  the  current  run,  as  shown 
in  lines  4  149-4153.  There  are  no  operator-procecessor 
'interactions   during   execution   of    these   subroutines. 

c.  Subroutine  DIWAv/INWAV/AUGWAV 

DIWAV,  INWAV  and  AUGWAV  are  similar  subroutines. 
They  drive  ACAVG  to  input  the  direct  waveform,   the  incident 
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waveform,  and  the  augmented  waveform.  They  allow  the 
selection  by  the  Operator  of  the  number  of  signals  to  be 
acquired  for  each  waveform  by  the  DPO  (total  composite 
waveforms  from  the  sampler)  and  averaged  by  the  OS.  Each  of 
these  subroutines  is  accompanied  by  a  message  on  CRT 
describing  the  procedure  to  be  employed  in  acquiring  the 
particular  signal.  These  messages  are  reproduced  in 
Appendix  F.  DIWAV  is  located  at  lines  3000-3170.  It  saves 
data  in  File  50  of  Tape  2.  INWAV  is  located  at  lines 
3200-3380.  It  saves  data  in  File  51.  AUGW&V  is  located  at 
lines  3400-3590.  It  saves  data  in  File  52. 
d.   Subroutine  ACAVG 

ACAVG  acquires  and  averages  the  appropriate 
waveform  interactively  requested  by  the  operator  from  the 
processor  keyboard.  It  also  removes  any  dc  signal  that 
might  be  included  in  the  acquisition.  This  feature  allows 
the  operator  to  position  the  desired  signal  anywhere  on  the 
DPO  screen  and  still  maintain  a  true  zero  axis  in  later 
graphic  displays.  Note,  the  signal  waveform  must  be 
completely  within  the  DPO  screen  area.  ACAVG  is  74-steps 
long  (23%  of  the  program)  and  is  located  in  statements 
2000-2730.  It  is  driven  by  one  of  the  appropriate  input 
waveform  subroutines   previously  described.     It   returns 
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results  to  that  subroutine  for  later  storage.  It  also 
displays  system  information  for  the  current  signal  being 
acquired  to    the   DPO   screen. 

Lines  2090-2130  initialize  &CAVG  by  requesting 
the  operator  to  input  the  number  of  waveform  signals  to  be 
averaged     in        N2.  This      value        is     than        returned      for 

verification.  In  lines  2140-2540,  the  OS  writes  the  current 
acquisiton  number  of  the  waveform  to  the  DPO  screen, 
acquires  and  stores  the  signal  into  DPO  memory  and  then 
transfers   the   signal   to   the   1052    for    further    processing. 

Lines  2170-2370  perform  the  count  function. 
Depending  on  the  number  of  waveforms  desired,  the  count  will 
be  for  every  signal  (N2  lass  than  21)  ,  for  every  five 
signals  (N2  less  than  10  1)  or  for  every  ten  signals  (N2 
greater   than    or   equal    to   101).  This   is   done    to   reduce  the 

time   required   for  acquisition   of   all   signals. 

Lines  24  10-243  0  instruct  the  DPO  to  store  and 
hold  the  waveform  in  DPO  memory.  Line  2420  stops  processing 
to  allow  for  the  DPO  to  fully  acquire  the  waveform  in  a 
single   sweep.      T1    is    the   length    of    the    wait. 

In  lines  2440-2450,  the  GS  instructs  the  DPO  to 
transfer  the  digitized  signal  from  DPO  memory  to  the 
microprocessor   temporary   storage   array   X0.         Lines   2470-2500 
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then  search  the  first  eight  values  of  the  array  for  any 
obviously  out  of  range  values  input  by  poor  DPO  blanking  on 
the  signal  front  end.  If  such  a  value  is  found,  the  array 
element  is  set  to  discrete  value  511,  corresponding  to  a 
scaled  value  of  0.  The  waveform  is  next  added  to  the  sum  of 
the  previous  waveforms  for  the  current  acquisition  in  line 
2530.  The  resultant  matrix  is  the  time  sampled  values  of 
the  signal  in  BO: 

N  " 

BO  =  Z      X0.  (4.1) 

i=l    x 

If  further  acguistions  are  necessary,  the  entire  process 
described  above  is  repeated.  If  no  further  acquisitons  are 
to  be  done  on  the  current  signal,  the  ensemble  average  time 
sampling  matrix  X0  is  determined  in  line  2550: 

(4.2) 
X0  =  B0/N2 

iny  dc  offset  is  eliminated  in  line  2560. 

It  is  interesting  to  note  the  time  required  to 
acquire  and  average  signals  using  the  solid  state  CJHF  pulse 
generator.  For  ten  averages  of  a  signal,  42  seconds  of 
elapsed  time  is  needed,  or  about  4.2  seconds  per  acquistion. 
Of  this,  about  1  second  is  needed  to  write  information  to 
the  DPO  screen  per  acquisition.    For  21  averages,   about  65 
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seconds  of  acquisition  time  is  required,  or  about  3.1 
seconds  pec  acquisition.  This  is  about  a  35%  reduction  in 
acquisition  time  overall  per  signal.  The  improved  time  per 
acquisition  is  due  to  the  reduced  frequency  with  which  data 
in  the  form  of  the  wavecount  is  written  to  the  DPO  screen. 
For  a  relatively  small  increased  investment  of  time,  a 
tremendous  advantage  is  realized  in  averaging  signals  and 
reducing  polluting  noise.  A  further  decrease  in  acquisition 
time  per  signal  is  realized  whan  10  1  or  more  averages  are 
specified.  However,  as  will  be  discussed  in  Chapter  7,  lack 
of  total  signal  stability  may  produce  errorneous  results  for 
such  a  large  number  of  acquisitions  due  to  the  length  of 
time  required. 

In  lines  2570-2820,  the  waveform  scale  factors 
are  determined  and  transferred  to  the  processor.  The  values 
are  returned  in  M$  (vertical  magnitude  scale  factor)  and  T^ 
(horizontal  time  scale  factor) .  The  discrete,  digitized 
verticle  values  are  appropriately  scaled  to  the  correct  real 
time  values  in  line  2700.  The  time  between  a  single 
discrete  sample  on  the  horizontal  axis  (sampling  rate)  is 
returned  in  Z1  at  line  2820. 

In  summary,  ACA7G  does  the  physical  acquisition, 
averaging   and  scaling   of   each   signal.    It   returns   the 
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current  waveform  count  to  the  DPO  CRT  for  display,  and  the 
appropriately  scaled  waveform  in  XO  to  the  driving 
subroutine  for  storage  in  the  appropriate  file.  It  further 
returns  the  appropriate  horizontal  and  vertical  axis  labling 
scale  factors  in  T  and  M  respectively.  The  flow  diagram 
for  ACAVG  is  provided  in  Figure  4.5 

The  basis  for  this  subroutine  was  derived  by 
Capt.  Hammond.  [Ref.  19]  His  program  in  its  entirety  has 
been  extensively  modified  to  fit  the  current  requirements  of 
the  TDRL  OS. 

e.   Subroutine  NEWPAR/ANTIN/TGTIN 

NEWPAR,  ANTIN  and  TGTIN  are  interactive 
subroutines  which  allow  an  operator  to  input  new  parameters 
for  defined  antennas  and  targets  into  the  mass  storage 
library  on  Data  Storage  Tape  2. 

NEWPAR  is  the  driving  routine  for  ANTIN  and 
TGTIN.  It  receives  its  input  from  the  DRIVER  and  selects 
either  ANTIN  or  TGTIN.  It  returns  the  new  parameters  to  the 
appropriate  file  on  the  storage  tape.  NEWPAR  then  determines 
if  further  inputs  are  required.  If  not,  program  control  is 
returned  to  DRIVER. 

ANTIN  defines  the  values  of  the  parameters  of 
the  new   antenna.    It   requests  the   operator  to   input  the 
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antenna  number.  An  array  of  63  elements  is  then  generated, 
each  third  element  being  fillsd  with  an  initial  antenna  run 
number  equal  to  1.  All  other  elements  are  set  to  zero.  If 
the  operator  specifies  an  antenna  value  greater  than  20,  the 
following  message  is  returned. 

"STORAGE  AREA  LIMITED  TO  20  ANTENNAS.  WANT  TO  CONTINUE?" 
This  is   necessary  as   current  space   limitations  limit   the 
number  of  antennas  to  be  defined. 

TGTIN  defines  the  values  of  the  parameters  of 
new  targets.  It  requests  the  operator  to  input  the  new 
target  number  as  well  as  the  length  and  radius  of  the 
target.  The  carraige  return  must  be  pressed  after  each 
entry. 

NEflPAR,   ANTIN   and  TGTIN   are  Located  in  lines 
4200-4800.   They  comprise  11%  of  the  program, 
f.   Subroutine  SCATWA7 

SCAT»AV  determines  the  target  backscattered 
waveform  based  on  the  differance  of  the  augmented  and 
incident  waveforms.  It  returns  the  backscattered  waveform 
to  storage  at  file  53  on  Tape  2.  It  also  sets  flag  P7  to  the 
appropriate  value  required  for  either  Prony  processing 
(P7  =  1)  of  Time  Domain  processing  (p7=0) .  Program  control 
is  then  returned  to  INPUT  DRI7ER. 
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2.   Program  MATH 

MATH  performs  most  of  the  significant  mathmatical 
computations  for  tha  OS.  Its  correct  performance  is  central 
in  the  production  of  accurately  processed  signals.  MATH  is 
stored  on  Tape  1,  File  3.  it  requires  14544-bytes  of  CPU 
memory  space.  With  all  arrays  and  constants  defined,  a 
maximum  total  of  49108-bytes  of  memory  are  required. 

This  sub-program  is  initiated  under  program  control 
only  from  INPOT  and  GRAPH.  There  is  no  interaction  between 
operator  and  system  during  MATH'S  execution. 

Inputs  to  MATH  are  from  INPUT,  GRAPH  or  Tape  2. 
Values  are  returned  to  mass  strorage  or  GRAPH. 

MATH  performs  the  following  functions  for  the  OS: 

1.  Fast  Fourier  transforms  the   Direct  and  Backscartered 
waveforms. 

2.  Determines  the  frequency  domain  impulse  response. 

3.  Performs  inverse  Fourier  transform  to  find  the  time 
domain  in pulse  response. 

4.  Determines  the  step  and  ramp  response  for  the  target. 

5.  Determines  the  physical  optics  shape  of  the  target. 

6.  Optimally   filters   the   frequency   domain   impulse 

response,   yielding   a  smoothed  time  domain   impulse 

response,  ramp  response  and  physical  optics  shape   of 

the  target  under  noisy  conditions. 
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7.  Compensates  the  ramp  response  and  optics  shape  for  EM 
propagation  attenuation. 

8.  Updates  optimal  antenna  parameters  for  the  antenna 
being  used,  based  on  the  targets  ramp  response, 
automatically. 

MATH  is  composed  of  four  subroutines  and  DRIVER.    A 
description  of  each  follows.  A  general  flow  diagram  for  MATH 
is  given  in  Figure  U.6. 
a.   MATH  DRIVER 

The  driver  routine  for  MATH  is  located  between 
lines  100  to  1210  and  comprises  59%  of  the  total  program. 
In  addition  to  initializing  the  program  and  driving  the 
subroutines,  DRIVER  performs  most  of  the  significant 
mathematical  computations  in  straight  line  data  flow 
fashion. 

MATH  DRIVER  is  initiated  at  line  100  from  INPUT 
or  line  130  from  GRAPH.  Parametric  data  is  input  directly 
from  INPUT.  It  is  recovered  from  mass  storage  in  line  160 
if  GRAPH  initiates  MATH.  Lines  100-221  serve  to  initialize 
the  local  program  parameters. 

Line  250  inputs  into  X0  from  mass  storage  file 
53  the  signal  (backscattered  waveform)  which  represents  the 
convolved  transient   response  of   the  transmitting  antenna, 
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the  receiving  antenna,  and  the  target.  This  data  is  then 
used  to  drive  subroutine  FFT.  XO  is  returned  from  FFT  in  a 
and  P  as  the  frequency  domain  magnitude  and  phase  and  saved 
in  arrays  E  and  S3  respectively. 

A  similar  procedure  is  performed  in  lines 
290-340  for  the  convolved  transient  response  of  the 
transmitting  and  receiving  antennas  without  target  on  the 
ground  plane.  The  magnitude  and  phase  are  returned  in 
arrays  F  and  F3  respectively. 

In  lines  360-420,  the  frequency  domain  impulse 
response  of  the  target  and  noise  optimization  are  performed, 
aagnitude  of  the  impulse  response  is  returned  in  SO,  with 
phase  returned  in  F0.  The  algorithm  for  determining  the 
impulse  response  and  performing  optimization  were  developed 
in  chapter  III  and  will  not  be  repeated  here. 

Lines  440-540  perform  the  inverse  Fourier 
transform  of  the  frequency  domain  impulse  response  of  the 
target  to  find  its  time  domain  equivalent.  The  frequency 
domain  data  entering  at  line  440  is  in  polar  form.  This 
form  is  most  suitable  for  the  previously  described 
mathematical  manipulations.  It  is  not  a  suitable  form  for 
the  inverse  Fourier  transform  method  employed  by  the  4052 
system.    In  order  to  perform  this  operation,   the  magnitude 
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and  phase  arrays  must  first  be  converted  into  two  source 
arrays  containing  real  and  imaginary  data.  Next,  the  first 
element  of  the  real  array  must  contain  the  signal  dc  term, 
while  the  last  real  element  must  be  the  Nyguist  freguency 
value.  As  both  these  elements  represent  pure  real  numbers, 
the  first  and  last  elements  of  the  imaginary  source  array 
muse  be  correspondingly  zero. 

Lines  440-470  establish  a  FOR/NEXT  loop  that 
returns  real  terms  to  M  and  imaginary  terms  to  P.  Real 
terms  are  established  by  multiplying  the  polar  magnitude, 
EO ,  by  the  cosine  of  the  corresponding  phase  FO: 

M(j)  =  E0(j)  x  (cos  (F0(j))  (4.3) 

Imaginary  terms  are  established  by  multiplying  the  polar 
magnitude  by  the  sine  of  the  corresponding  phase  term: 

P(j)  =  E0(j)  x  sin(F0(j))  t4-4) 

INPUT  had  previously  removed  all  dc  levels  prior 
to  MATH,  so  S(1)  and  M  (257)  ,  the  dc  and  Nyguist  terms,  are 
set  equal  to  zero  in  lines  480-490.  Lines  500-510  set  the 
corresponding  imaginary  terms  to  zero.  Note  that  both  of 
these  arrays  must  be  257  data  points  long,  as  the  output 
array  will  be  512  data  points. 
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Data  is  now  correctly  formatted  in  the  two 
source  arrays  so  that  the  command  "INL2AV"  can  interleave 
the  real  and  imaginary  data  into  a  single  destination  array, 
10,  whose  format  is  acceptable  to  the  GS   "IFT"  command. 

Line  540  calls  the  system  inverse  Fourier 
transform  command  "IFT"  which  transforms  the  spectral  data 
into  waveform  (time-domain)  data.  The  spectral  data 
originally  in  10  is  overwritten  by  the  new  time  domain  data 
and  returned  in  10.  10  is  now  the  time  domain  impulse 
response. 

The  inverse  Fourier  transform  can  be  expressed 
mathematically  by  the  following  summation: 


X(n)  =  (1/N)   V  xd(k)  eJ2™k/N   ,  for  n-0.1, . . . ,N-1   (4.5) 

K  =  U 


In  the  above  equation,  N  refers  to  length  of  the  array 
argument,  and  n  is  an  index  used  in  generating  the  various 
elements  of  the  output  array.  X  (k)  is  the  k-th  complex 
Fourier  coefficient,  and  X (n)  is  the  n  +  1-th  element  of  the 
real  data  output.  [ Hef .  20]  The  4052  uses  an  extremely 
fast  method  for  IFT  computations,  the  Sande-Tukey 
decimalization- in- frequency   algorithm. 
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Line  620-720  save  the  time  domain  impulse 
response  (file  54) ,  and  frequency  domain  direct  waveform 
(file  55)  ,  the  backscatt ered  waveform  (file  56)  and  the 
impulse   response    (file    57)     in   mass   storage. 

Next,  the  time  domain  step  response  of  the 
impulse  is  returned  to  X0  in  line  740  by  using  the  4052 
"INT"  routine.  This  is  an  integration  calculated  using  the 
trapezoidal  rule  for  approximating  the  definite  integral  as 
follows: 

B    =    0 
B     (t)     =    B(t-1)     ♦    .5    *     (A(t-1)     ♦    Mt))  (^  ^ 

for  t  =  2,3, ...,N 
where:    A  -    the  source  array 

3  =  the  destination  array  (result) 
N  =  the  number  of  elements  in  the  array 
Any  dc  value  in  the  result  is  removed  through  the  actions  of 
lines  750-760. 

Next,  the  time  domain  ramp  response  is 
determined  by  integrating  the  step  response  in  X0.  The  ramp 
is  returned  in  RO.  The  MATH  DRIVER  outputs  the  time  domain 
ramp  and  its  minimum  value  in  R2  to  Subroutine  OPTICS  which 
returns  the  physical  optics  shape  (shape)  and  ramp  response 
compensated   for   electromagnetic   propagation   attenuation. 
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Note  the   shape  does  not  as   yet  have  the   correct  amplitude 
values  in  its  array  elements. 

If  the  run  is  the  first  for  the  target  (prior  to 
optimization)  and  the  target  is  one  defined  in  mass  storage, 
MATH  DRIVER  inputs  the  last  computed  antenna  parameters  for 
the  compensator  for  shape  length  (S5)  and  radius  (S6) ,  as 
well  as  the  actual  defined  lenghth  (C6)  and  radius  (C5)  to 
subroutine  ANTUP.  S5  and  S6  are  then  automatically  updated. 
The  physical  optics  array  is  then  multiplied  by  the 
compensator  value  S5  so  that  its  elements  reflect  the 
estimated  dimensions  of  the  targets  in  inches.  This  is 
based  on  the  experimentally  observed  fact  that  a  particular 
size  and  configuration  of  stationary  target  will  return  some 
baclcscattered  field  relative  to  its  size  and  shape.  That 
this  value  is  not  exactly  the  same  from  one  run  to  the  next, 
but  only  approximately  so,  is  due  to  certain  stastically 
determinable  effects  produced  by  noise,  small  differences  in 
orientation,  etc.  A  large  number  of  runs  should  produce  an 
averaging  value  which  takes  into  account  these  factors.  If 
the  target  is  unknown  (the  general  case) ,  or  the  current  run 
is  for  noise  reduction,  previously  computed  values  of  the 
antenna  parameters  are  recalled  from  mass  storage.  For  the 
unknown  target  these   compensators  are  the  averages   for  all 
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targets  for  all  runs  for  a  particular  antenna.  In  fact, 
OPTICS  is  not  part  of  the  data  processing  flow,  the 
compensators  having  been  previously  determined  in  INPUT  at 
lines  1024-1025. 

Following  a  known  target  run,  the  minimum  value 
of  BO  is  returned  to  B8.  If  tha  run  is  not  for  optimization 
(L1=1),  then  36,  an  amplitude  compensation  optimization 
factor  is  set  equal  to  38  and  the  actual  length  of  the  shape 
is  retained  in  Q9.  For  graphing  purposes,  it  is  desirable 
that  B6  and  Q9  not  vary  from  optimization  run  to 
optimization  run.  MATH  DRIVER  will  input  BO,  Q9,  and  B6  to 
subroutine  OPTIM,  which  compensates  the  amplitude  and  length 
of  the  physical  optics  shape  for  the  low  pass  filtering  done 
to  reduce  noise.  An  arbitrary  target  is  treated  in  exactly 
the  same  manner.  (Note,  this  effects  the  shape  ONLY  and  not 
the  time  domain  ramp  response  during  optimization.) 

The  above  system  has  worked  remarkably  well  in 
deducing  the  actual  physical  parameters  of  both  known  and 
unknown  targets.  ?or  known  targets,  accuracy  is  invariably 
10%  or  less.  For  unknown  targets,  this  accuracy 
deteriorates  to  about  15%  or  less.  If  the  data  is  treated 
in  the  raw  form,  with  no  compensation  applied,  accuracy  is 
reduced  to  25%  or  worse,   up   to  50%  for  the  unknown  target. 
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This  can  be  noted  on  the  graphic  output  of  the  shape  where 
the  compensated,  uncompensated  and  known  (when  given)  values 
of  length  are  compared  run  to  run. 

Lines  1000-1140  save  the  updated  antenna 
parameters  in  the  appropriate  files,  the  shape  array  in  file 
58,  the  time  domain  ramp  response  in  file  59,  and 
significant  system  parameters  in  file  60,  all  on  mass  data 
storage  Tape  2. 

Lines  1160-1190  is  MATH  DRIVER  ending  routine, 
deleting  variables  and  passing  control  of  the  program  to 
GRAPH  (file  4)  on  Tape  1. 

In  summary,  MATH  DRIVER  is  the  principal  vehicle 
for  OS  processing  for  inverse  scattering.  A  flow  diagram  is 
given  in  Figure  4.7. 

b.   Subroutine  FFT 

FFT  is  the  fast  Fourier  transform  subroutine 
used  tc  convert  the  time  domain  Direct  and  Backscattered 
waveforms  to  spectral  data.  This  subroutine  is  located  in 
lines  2000-2140  and  requires  1%   of  the  total  program  space. 

X0  contains  the  real-valued  signal  data  on  which 
the  Fourier  transform  is  to  be  performed.  It  is  input  from 
MATH  DRIVER.  Line  2060  calls  the  GS  command  "FFT"  and 
performs  a   fast  calculation  (about   4  seconds)   of   the  DFT 
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(discrete  Fourier  Transform) ,   expressed   mathematically  by 
the  summation:   [Ref.  21] 


N-l 

Xd(k)    =      Z        X(n)ej27Tnk/N      ,    for  k=0,l N/2  (4.7) 

n=0 


where:      N  ~   array  length 

k  =  index  used  in  generating  various 
Fourier  coefficients 

X. (k)  =  k-th  Fourier  coefficient 

X(n)  =  n  ♦  1-th  elemant  of  the  real  data 
input 

After  the  command  has  executed,  the  original 
real  data  will  be  overwritten  by  the  results  of  the  FFT 
computation.  The  output  data  is  in  complex  rectangular  form 
with  the  first  element  of  the  array  XO  containing  a  real 
number  representing  the  dc  value  of  the  signal,  the  second 
element  containing  a  real  number  representing  the  value  of 
the  discrete  Fourier  transform  at  the  Nyquist  frequency,  and 
the  remaining  array  elements  containing  alternately  real  and 
imaginary  Fourier  coefficients. 

This  arrangement  is  not  suitable  for  the 
division  that  later  takes  place.  In  line  2070,  the 
rectangular  array  XO  is  converted  to  polar  form  by  calling 
the  GS  command  "POLAR".     Magnitude  components  are  returned 
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in  array  H  and  phase  components  are  returned  in  array  P.   At 
this  point,   subroutine   FFT  returns  control  to   JiATH  DRIVER 
with  the  arrays  XO,  a  and  P  for  further  processing, 
c.   Subroutine  OPTIC 

Subroutine  OPTIC  manipulates  the  ramp  response 
of  the  target  to  recover  the  physical  optics  shape,  but  it 
does  so  without  effecting  in  anyway  the  raw  data.  It  later 
compensates  both  the  ramp  response  and  the  shape  for  EM 
propagation  attenuation  experienced  along  the  body  of  the 
target.  Subroutine  OPTIC  resides  in  lines  2200-2700  and 
comprises  25%  of  Program  MATH. 

The  ramp  response  R0  and  the  location  of  its 
minimum  value  within  the  array  are  input  into  OPTIC  from 
HATH  DRIVES.  The  mean  value  of  the  response  is  calculated 
in  lines  2250-2260  and  returned  in  R3.  This  value  is  used 
later  in  determining  the  endpoints  of  the  shape  array. 

Lines  229  0-25  00  is  a  FOR/NEXT  loop  that 
seperates  out  the  targets  shape  from  the  ramp  array  data. 
The  search  for  valid  data  begins  at  the  point  in  the  array 
equal  to  the  location  cf  the  minimum  value.  It  is  assumed 
that  the  target  will  always  return  a  response  larger  than 
any  that  might  be  due  to  range  noise  or  other  excitation. 
For  the  size  of  the   targets   being  usei,   and   the  the 
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configuration  of  the  present  source,  this  has  always  proven 
to    be      an   accurate   assumption.  The   FOR/NEXT      loop   searchs 

first      for      a      trailing-edge   endpoint.  This      is     done      in 

ascending  time  steps  in  search  of  the  point  where  the  array 
elements  become  zero  or  positive,  representing  a  zero-axis 
crossing.  For      every     value      of    RO      less      than      zero,        a 

corresponding  value  is  input  in  BO,  the  shape  array.  If 
30  (k)  is  a  positive  value,  representing  a  zero-crossing,  the 
FOR/NEXT  loop  is  terminated.  If  BO  (k)  is  negative,  OPTIC 
then  checks  the  location  of  the  point.  It  is  important  to 
do  this  as  noise  and  instability  of  input  signals  from  an 
acquisition-to-acquisition  run  may  erroneously  shift-zero 
crossings  or  preclude  zero-crossings  altogether.  In  these 
cases,  determination  of  the  shape  must  be  somewhat  modified 
as    described   below. 

In  line  2350,  two  situations  are  checked. 
First,  is  the  preceding  array  element  of  RO  more  negative 
than  the  present  element  and  is  the  present  element  greater 
that      90%        of      the      ramp        response      average        value.  If 

affirmative,  the  next  array  element  is  fetched.  A  positive 
result  of  this  check  indicates  that  datum  is  on  the  skirt  of 
the  ramp  response,  and  that  the  data  magnitudes  are 
progressing    towards   a    zero-axis      crossing   normally.        if    the 
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check  returns  as  negative,  that  is,  the  present  array  value 
is  more  negative  than  the  preceding  element,  than  an  anomaly 
has  occurred  either  due  to  noise,  or  time-axis  shifting.  In 
this  case  a  special  endpoint  routine  must  be  employed. 

The  second  check  is  to  determine  if  the  array 
element  is  on  the  main  signal  portion  of  the  response  and 
not  the  skirts.  If  the  element  amplitude  is  greater  that 
90%  of  the  mean  value,  than  it  is  assumed  that  the  element 
is  on  on  the  skirt  ir regardless  of  the  relationship  between 
the  amplitude  of  the  present  and  preceding  elements.  The 
next  datum  is  fetched.  It  can  be  seen  that  this  arrangement 
is  subject  to  some  error  if  a  noisy  element  is  encountered 
just  at  the  point  where  the  ramp  response  passes  the  90-th 
percentile  of  the  signals  mean  value.  This  error  is 
considered  relatively  minor. 

The  90-th  percentile  was  chosen  following 
numerous  observations  that  noise  or  time-shift  in  the  main 
ramp  response  normally  do  not  produce  fluctuations  that 
cross  this  boundary.  Zrror  in  the  shape  representation  will 
be  introduced  if  it  does,  in  fact,  cross  this  levei.  The 
result  will  be  a  fore- shortened  drawing  of  the  target. 

If  the  need  for  an  endpoint  routine  is 
indicated,   processing   goes  to   lines  2360-2^70.    In  line 
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2360-2380,  the  distance  the  last  valid  data  element 
(R0(k-1))  is  from  the  zero  axis  is  determined.  In  line 
2360r  the  distance  between  the  previous  two  valid  data 
elements  is  found  and  returned  in  W(3): 

W(3)  =  R0(k-2)  -  R0(k-1)  (*H  •  «) 

The  above  is  for  the  ascending  time  series  run  of  the 
FOR/NEXT  loop  (trailing  edge) .  For  the  descending  time 
series,  the  values  in  the  parenthesis  have  opposite  sign. 
This  difference  represents  an  approximation  of  the  slope  of 
the   skirt.  A    problem  with      this   approximation   is      that    it 

often  occurs  where  the  slope  is  not  as  great  as  the  average. 
The     result    is      a      somewhat   elongated      figure.  It   may      be 

possible  to  defeat  this  problem  by  talcing  more  than  two 
points    for   the   average. 

In  lines  2370-2380,  a  check  is  made  to  see  if 
the  last  valid  data  point  is  within  10-times  or  less  of  the 
difference      determined      above.  If        it      is,        the      actual 

difference  is  returned  in  Q1  and  used  to  drive  the  values 
placed  in  BO  to  a  zero-crossing.  If  the  difference  is  less 
that  10%  of  the  amplitude  needed  to  effect  a  zero-crossing, 
then  one-tenth  of  the  actual  distance  is  returned  in  Q1.  In 
this      manner,    the  shape  rapidly    goes    to   zero    and  an   endpoint 
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is  determined  close  to  the  estimated  real  endpoint  had  a 
natural  zero-crossing  occured.  The  endpoint  is  located  with 
a  FOR/NEXT  loop  between  lines  2420-2470.  This  FOR/NEXT  loop 
simply  iterates  the  values  of  BO  by  Q1  until  a  zero-crossing 
is  effected.  Whenever  BO  becomes  positive,  all  loop 
processing  is  terminated  and  the  endpoint  is  returned  in 
W(1)  for  the  trailing  skirt  and  »(2)  for  the  leading  skirt. 

If  the  above  process  was  for  the  leading  edge  of 
the  ramp  response,  the  loop  is  reentered  and  the  trailing 
edge  check  of  the  data  is  performed.  The  method  employed  is 
identical  with  the  above  process  for  the  leading  edge.  Note 
that  up  to  this  point,  the  ramp  response  in  R0  has  been 
totally  unaffected  by  any  processing  in  OPTICS.  It  has  been 
used  simply  as  source  of  data. 

Once  the  target  response  has  been  located  in  R0 
and  entered  in  BO,  then  both  arrays  can  be  compensated  for 
propagation  losses  experienced.  This  is  done  in  lines 
2570-2640.  The  value  of  the  nominal  distance  the  target  is 
from  the  transmitting  antenna  is  returned  in  Z6.  A  FOR/NEXT 
loop  in  lines  2590-2640  corrects  the  arrays  for  the 
estimated  propagation  losses.  Line  2500  computes  the 
distance  along  the  body  the  wave  propagates  for  each  sample 
point  in  R0: 


135 


S    =    (Q    -    5?(1)     *    Z1    *     (10*8)/2 


Cm.<?> 


where:     (Q   -    W(1))    *    Z1    =    twice   the    time     of   travel 

of   the   propagated   wave 

In   line   2610,        the   propagation   loss   compensation      factor    is 
returned   using: 


U   -    CZ  +  RH2Z+R)  (4.10) 

22Z 


where:  Z  =  range  from  transmitting  antenna  to 
target  face 
It  is  also  assumed  that  the  distance  between  transmitting 
and  receiving  antenna  is  Z.  Z  is  in  fact  the  measured  value 
of  distance  read  from  the  "TIME/DISTANCE"  dial  on  the  DPO. 
Each  value  of  R0  and  BO  is  then  multiplied  by  this 
compensation  factor  in  turn.  Processing  continues  once  all 
values  for  the  computed  length  have  been  appropriately 
amended. 

Line  2260  is  a  FOR/NEXT  loop  that  transfers  the 
values  in  the  array  30  to  the  beginning  of  the  array  to 
facilitate  graphing.  Once  this  is  accomplished,  program 
control  is  returned  to  MATH  DRIVER  with  the  shape  unsealed 
for  actual   dimensions  of   the  target   (values  are  equal  in 
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amplitude  to  the  compensated  voltage  magnitudes  of  the  ramp 
response) .  The  ramp  resonss  is  returned  unef f ected  by  any 
of  the  above  processing  except  for  propagation  attenuation 
compensations.  Figure  4.8  provides  a  flow  diagram  for 
Subroutine  OPTICS. 

d.   Subroutine  AiliaP 

This  surboutine  updates  the  specific  antenna 
length  and  radius  compensation  values  if  the  target  is  one 
defined  from  storage  on  Tape  2.  If  the  target  is  arbitrary 
in  the  general  case,  no  updating  occurs.  ANTUP  is  located 
between  statements  3000-3080  and  comprises  H%  of  the 
program. 

The  unsealed  shape  array  B0  is  input  to  ANTUP 
from  HATH  DRIVER.  In  line  3060,  the  amplitude  compensating 
and  scaling  factor  is  computed  using; 

36=  (A  +  S*B)/(C*B) 

where:   S6  =  scaling/compensating  amplitude  factor 

C  ■  radius  of  target 

B  =  minimum  magnitude  of   physical  optics 
array 

A  =  the   number   of  previous  times   this 
antenna/target  combination  was  used 

S  =  sum  of  previous  compensation  values 
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The  result  is  simply  the  averaged  value  of 
compensation/scaling  for  ail  runs  using  this  particular 
antenna  and  target  combination. 

Line  3070  performs  a  similar  computation  for  S5, 
updating  the  compensation/scaling  factor  for  the  target 
length.  The  formula  is  identical  except  that  the  defined 
target  length  replaces  the  target  radius  and  the  observed  or 
measured  length  from  the  plane  replaces  the  minimum 
magnitude  of  the  shape  array.  The  concept  is  that  slight 
differences  in  target  returns  due  to  variations  in  target 
location  and  orientation,  noise,  etc.,  will  be  averaged  out 
over  the  long  term  for  any  particular  antenna.  When  the 
target  is  arbitrary  and  thare  is  no  a  priori  information 
then  ANTOP  is  bypassed  and  the  average  of  all  the 
compensation/scaling  factors  for  all  targets  for  the 
particular  antenna  is  used.  This  requires  that  any  new 
antenna  be  "seeded"  with  compensation/scaling  factors  from 
at  least  one  run  using  known  targets  prior  to  attempting  to 
run  an  arbitrary  target  on  the  antenna.  What  results  is  a 
"poor  man's"  compensator  for  the  currents  on  the  body  of  the 
target  based  on  empiracal  observations  of  many  target 
responses.  It  is  not  as  accurate  as  the  full  solution  of 
the  time   domain  integral   equations  would   be,   nor  is  the 
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exact  target  shape  returned.  What  is  obtained  is  an 
excellent  estimate  of  the  target,  whether  known  of  unknown, 
based  on  an  historical  compilation  of  data  for  the  TDRL 
system.  This  concept  can  be  generalized  to  any  system  using 
components  which  are  nominally  fixed,  such  as  large  scale 
bistatic  radars.  The  shape  itself,  at  least  for  simple 
ax  isy  in  metric  targets,  can  probably  be  returned  with  the 
development  of  the  proper  software  routines.  Perhaps  the 
rate  at  which  the  ramp  response  slopes  vary  may  be  analyzed 
for  this  purpose.  What  is  obtained  is  a  method  which 
employes  a  small  scale  processor  requiring  relatively  simple 
processing  to  obtain  accurate  time  domain  estimates  based  on 
the  ramp  response  of  the  target  almost  in  real  time,  and 
capable  of  being  generalized  to  almost  any  system  employing 
the  proper  source. 

The  compensation/scaling  factors  are  returned  to 
MATH  DRIVER  along  with  program  control.  30  is  unaffected  as 
yet  by  S5  or  S6. 

e.   Subroutine  OPTIM 

OPTIM  compensates  the  shape  array  for  the 
effects  of  low-pass  filtering  of  the  target  responses  using 
Riad's  Optimization  technique.  Optimization  is  done  to 
remove  high-frequency  noise  terms  and,   for  the   shape,   to 


140 


improve  the  overall  shape  of  the  drawing.  In  reality, 
without  compensation,  the  drawing  is  less  representative  of 
the  true  targets  shape  when  optimized.  Low  pass  filtering 
will  lengthen  the  overall  shape  and  reduce  the  radius,  the 
amount  being  relative  to  the  value  of  lambda,  the 
optimization  variable.  OPTIM  allows  for  this  and  introduces 
factors  which  will  pin  the  optimized  shape  to  the  original 
form  while  allowing  the  benefits  (reduced  noise  ripple  on 
the  body)  of  low-pass  filter  noise  reduction  to  be  shown. 
OPTIM  is  located  in  lines  3100-3180  and  represents  H%  of  the 
program. 

30,  B6,  (original  minimum  of  the  unoptimized 
shape) ,  38  (the  present  minimum  of  the  shape)  and  Q9  (the 
original  unoptimized  target  length)  are  input  from  MATH 
DRIVER.  The  amplitude  OPTIM  compensation  is  found  in  line 
3140  by  dividing  36  by  38.  It  is  returned  in  B7.  The  new 
compensated  target  length  factor  is  returned  in  Q7  and  is 
found  by  dividing  the  present  value  of  the  length  by  the 
original  value.  Finally  the  whole  physical  optics  shape 
array  is  scaled  by  the  new  amplitude  factor.  The  new  values 
are  returned  to  MATH  DRIVER,  along  with  the  original  factors 
in  36  and  Q9 .  In  fact,  after  the  initial  unoptimized  run, 
B6  and  Q9  are  constants  until  a  new  target  is  acquired. 
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3*      Program   GRAPH 

GRAPH   provides      a    central   means  to      visually   display 

the   results    of    processing   performed      in   MATH    .        The  correct 

evaluation   and   analysis   of   the      acquired  processed   data   is   a 

direct      function      of      the        effective      reproduction     of      the 

resultant      arrays     by    GRAPH      into     a      form    most      useful      and 

understandable   to  the   operator.      GRAPH   is   located  on   Tape    1, 

File   <i ,    and   requires    44712-bytes   of    memory   space. 

Program   GRAPH      is    normally      initiated    under      program 

control   from    Program    MATH,    in   line    100.      It    may    be   initiated 

directly   by      the   operator   if   the      results   of    a      previous   run 

have   been   retained  in    mass      storage.      The   commands   issued    by 

the  operator   are: 

FIND    4 

CALL    BOLD 

ROM 

Outputs   are    to    the  CRT    or    a     return   of   control    to  one   of   the 

other    programs      that    make-up  the      TDRL   OS.         No      results   are 

returned   to    mass   storage. 

GRAPH      subroutines        may      be         divided      into        three 

categories: 

1.      DRIV2R--provides      for   initialization      of   the     program 

and  drives   subroutines.    It  also  passes  control  to 

the   next    program  in   the   OS. 
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2.  AXIS  and  Labling  Subroutines  —  establishes  the  axis  of 
all  graphs  and  lables  as  appropriate.  Provides  for 
formatting  of  descriptive  parameters,  as  appropriate 
builds  grids  for  graphics.  Six  subroutines  comprise 
this  group.  They  are  AXIS,  TICKS,  AXLAB,  TITLE, 
GRPAR,  GRID. 

3.  Inverse  Scattering—Subroutines  in  this  section  graph 
the  results  of  processing  in  INPUT  and  MATH.  They 
reproduce  the  four  basic  input  waveforms,  the  Fourier 
transforms  of  the  Direct,  3ac)cscatter  and  Impulse 
Response,  the  Physical  Optics  Shape  of  the  target, 
and  the  time  domain  ramp  response.  Four  subroutines 
comprise  this  group.  They  are  GRIN,  GRFOUR,  GROP, 
and  GRAMP. 

A  general  flow  diagram  for  GRAPH  is  provided  in 
Figure  4.9.  A  disscussion  of  all  Subroutines  in  GRAPH 
follows. 

a.   Graph  DRIVER 

As  for  INPUT   and  MATH  DRIVER,   GRAPH   DRIVER  is 

the  principal   initializer  and   controiing  routine  for  this 

program,   coordinating  the  flow  of  processing.    Unlike  MATH 

DRIVER,   GRAPH   DRIVER  does  not   of  itself  perform   any  data 

processing.   That   function  is  retained  for   the  subroutines 

exclusively. 
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GRAPH  DRIVER  is  located  in  lines  100-1190  and 
comprises  19%  of  the  total  program.  GRAPH  is  completely 
independant  of  data  input  from  the  other  programs  in  the  OS. 
All  inputs  are  from  mass  storaga  Tape  2  or  the  GS  keyboard. 
All  results  are  returned  to  CRT  display  of  the  U052.  This 
routine  is  highly  operator-processor  interactive.  Appendix 
?  describes  a  typical  run  with  messages  as  produced. 

GRAPH  DRIVER  is  always  initiated  in  line  100 
when  being  entered  external  to  the  GRAPH  Program.  It  is 
initiated  from  line  290  by  returns  from  internal 
subroutines.  Lines  100-300  initialize  the  program  and  input 
parameter  data  from  mass  storage.  Lines  320-750  develope 
the  correct  Menu,  with  the  function/subroutine  selected 
input  to  P0  in  line  76  0. 

The  Menu  that  is  used  in  GRAPH  lis  shown  in 
Apendix  F.  Most  of  the  fiunctions  on  the  menu  are 
self-explanatory.  A  few  remarks  concerning  items  10,  1 1, 
12,  and  13  should  suffice  to  rsmove  any  ambiguity. 

Item  10,  "CONTINUE  OPTIMIZATION'1,  initiates  the 
Riad  low-pass  filtering  optimization  of  the  impulse  response 
in  lines  950-1030.  Program  control  with  the  operator 
designated  optimization  factor  is  input  to  Program  MATH. 
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Item  11,  "NEXT  TARSET",  passes  program  control 
to  INPUT  for  the  acquisition  of  an  entirely  new  target. 
Parameters  are  retained  and  passed  along  with  program 
control. 

Item  12,  "FINISHED  GRAPHING",  terminates 
processing.  The  program  can  be  reinitiated  by  typing  "RUN". 

Lines  820-900  select  the  appropriate  subroutine 
and  functions.  Line  940  returns  to  the  menu  for  further 
selection  for  all  items  except  10,11,  12,  13. 

The  next  series  of  subroutines  to  be  discussed 
perform  the  graphics  axis  formatting,  labling  and  titling. 
These  subroutines  sat  the  environment  that  makes  the 
information  of  use  to  an  observer.  AXIS, PICKS,  AXLAB,  and 
TITLE  are  subroutines  based  on  programs  previously  developed 
by  Tektronix,  Inc.  [Ref.  22]  All  have  been  modified  to  a 
greater  or  lesser  degree  dependant  on  the  needs  of  the  TDRL. 
b.   Subroutines  AXIS/TICKS/AXLAB 

Subroutine  AXIS  drives  Subroutines  TICKS  and 
AXLA3.  They  perform  the  general  axis  definition  function 
for  the  output  graphics.  Subroutine  AXIS  serves  as  the 
imput  portal  for  data  from  the  function  subroutine  that  was 
chosen  by  the  operator.  The  results  returned  to  the  CRT  of 
the  GS   are  "neat"  ticks,   labled  and  scaled  to   the  proper 
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values  of  the  input  function.  Subroutine  AXIS  occupies 
statements  2000-2250,  TICKS  is  in  lines  2300-26  10  and  AXLAB 
is  in  statements  2700-3190.  These  three  subroutines  occupy 
18s?   of    program   space. 

Data  input  from  the  function  to  be  graphed  are 
the  minimum  and  maximum  horizontal  values  of  the  viewport  in 
V1  and  72,  and  tha  minimum  and  maximum  vertical  values  of 
the  viewport    in    V3   and    74.  These   constants   are   in   Graphic 

Display  Units  (GDO) .  GDfl's  are  a  measure  of  the  resolution 
capability  of  the  3S  screen.  Tha  four  arguments  defined 
above  refer  to  the  actual  limits  for  the  graphic  information 
of  the  display.  The  GS  will  draw  no  lines  outside  these 
limits. 

Also  input  to  Subroutine  AXIS  from  the  driving 
function  subroutine  are  the  minimum  and  maximum  horizontal 
values  of  tha  function  itself,  511  and  #2,  and  the  minimum 
and  maximum  vertical  values,  H3  and  514.  These  values  are  in 
User  Data  Dnits  (UDO's).  That  is,  in  units  that  the  user 
defines,    such   as    millivolts,    seconds,    meters,    etc. 

Finally,  iXIS  receives  the  valua  of  the  flag  SO. 
This  specifies  whether  the  vertical  axis  (S0=0)  or 
horizontal  axis  (S0=1)  is  being  labled.  The  magnitude  of  the 
horizontal    scaling   function   for   the    input   wave,    S1     (received 
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from  INPUT),   and  the  definitions  in  UDU's  for  the  vertical 
axis  and  horizontal  axis  are  required. 

The  performance   of  AXIS,   TICKS,   and   AXLAB  is 
adequately  described   in  the  Tektronix  source   referances. 
[Ref.  23]  This   report  will  describe  only   those  statements 
which  have  been   changed  or  modified  to   produce  the  desired 
output  for  TDRL  purposes. 

Principally,  the  changes  effect  the  manner  in 
which  the  axis  is  lable  with  UDO's  in  Subroutine  AXLAB. 
These  changes  lie  in  statements  2860-3190.  Three  conditions 
must  be  checked: 

1.  Horizontal  (S0=0)  or  vertical  axis  (S0=1)  lables? 

2.  Graphics  of  a  Fourier  transform?  (4<P0<8) 

3.  Check  amplitude   of  maximum   vertical  value   (P6)   in 
UDO's. 

TABLE  17  defines  the  relationship  between 
statements  in  Subroutine  AXLA3  and  the  functions  to  be 
graphed. 

c.   Subroutine  TITLE 

This  subroutine  positions  the  title  of  the 
graphic  above  the  plotted  values.  It  is  located  in  lines 
3200-3310  and  comprises  2%  of  the  program.  The  value  of  the 
maximum  horizontal  coordinate  92      and  the   graphic  computed 
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TABLE  IV 


Graphic  Labling  Relationships — GRAPH  AXLAB 


GRAPH 


STATEMENTS 


REMARKS 


GRIN 
GRAMP 


HORIZONTAL 
AXIS 


HORIZONTAL 
AXIS 


2880-2890 


S1*10>1000.      Lable    Axis 


29  10 


S1*10<1000.      Lable    Axis 


3090 


Position  Dimension  Lable 


3110 


Dimension  Lable 


3140-3160 


Print   Dimension   Lable 


3180-3190 


Return 


VERTICAL 
AXIS 


VERTICAL 
AXIS 


2860 


Start 


2940-2970 


P(6) <1.      Lable    axis 


2980-2990 


1<P(6) <10.      Lable   axis 


3000-3010 


P(6)>10.       Lable    axis 


3130-3150 


Print    Dimension    Lable 


3180-3190 


Return 


. — 


GRFOOR 


HORIZONTAL 
AXIS 


HORIZONTAL 
AXIS 


2861-2862 


2863-2864 


"Magnitude"    Axis    Label 


2865-2866 


"Phase"    Axis    Lable 


3171 


Print   Dimension    Lable 


3180-3190 


Return 


VERTICLE 
AXIS 


VERTICLE 
AXIS 


2360-3170 


Same   as   for   Horizontal   Axis 


3173-3190 


Return 


NOTE; 


All    graphic  routines   use   lines    3020-3040 
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vertical  maximum  P(6)  are  input  by  the  graphic  subroutine. 
The  graph  title  centered  over  the  graph  is  returned  to  the 
CRT. 

d.   Subroutine  GRPAR 

GRPAR  graphs  the  parameter  input  by  INPUT  and 
developed  by  MATH.  This  subroutine  is  located  in  lines 
3400-3870  and  comprises  8%  of  the  program.  All  graphs  have 
the  run  number,  date,  target  distance  from  transmitter, 
antenna  number,  target  number,  averages  taken,  optimization 
and  any  specific  remarks  printed.  The  specialized 
parameters  printed  for  particular  graphs  alone  with  the 
statement  number  applicable  are  given  in  TABLE  V. 

TABLE  V 
Specialized  Parameter  Prints — GRAPH  GRPAR 


GRAPHIC 


PARAMETERS 


LINES 


GRIN 
GRAMP 


G20P 


MAXIMUM  PEAK  VALUE 


3510-3600 


MINIMUM  PEAK  VALUE 


RMS  VALUE 


MEAN  VALUE 


DIAMETER  (INCHES) 


WIDTH  (INCHES) 


=t 


|  3710-3750 
H 


1 
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e.   Subroutine  GRID 

This  simple  subroutine  draws  referencing  grid 
lines  on  the  appropriate  graphic.  The  computed  graph 
maximum  and  minimum  are  input  from  the  graphic  subroutine 
along  with  the  axis  tick  intervals.  Grid  lines  are  output 
to  the  CRT  from  GRID.  Lines  3900-4010  contain  the 
subroutine.  GRID  requires  2%   of  program  space. 

This  concludes  the  discussion  of  general 
graphics  formating  subroutines.  The  driving  graphic 
subroutines  will  now  be  discussed. 

Subroutines  GRIN,  GRPOUR,  and  GRAMP  are  all 
similar  in  formatting  structure.  With  some  judicious  choice 
of  programming  flags  and  routine  definition,  it  should  not 
be  too  difficult  to  combine  these  three  subroutines  into  a 
single  driving  graphic  subroutine.  GROP  is  formatted 
differently,  but  still  retains  a  few  features  common  to  all 
the  graphic  driving  subroutines.  The  similarities  between 
the  subroutines  will  be  discussed  first.  The  particular 
differences  applicable  will  be  discussed  secondly. 

Each  graphic  subroutine  is  composed  of  four 
basic  blocks.   They  are: 

1.  Input  waveform  data  from  mass  btorage 

2.  Define  waveform  and  graphic  boundaries 
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3.   Draw  the  graph  of  the  function  using  special  defining 
functions  when  necessary 

i*.   Labling  the  graph  with  the  correct  title  (s) 

GRAPH  DRIVER  inputs  the  function  value  PO  which 
is  used  to  identify  the  location  of  the  data  storage  file  in 
which  the  appropriate  data  rests.  This  single  row  array  is 
usually  returned  in  X.  The  waveform  window  and  viewport 
parameters,  bounding  the  output  graphic  display,  are  next 
determined.  The  maximum  horizontal  axis  values  are  found  by 
mutiplying  the  distance  between  a  single  data  point  (in 
UDU's)  by  the  total  number  of  data  points  in  the  array  along 
the  X-axis.  For  those  waveforms  listing  the  mean  graph 
values  as  part  of  the  descriptive  parameters,  the  sum  of  all 
the  array  element  values  is  obtained  and  then  divided  by  the 
discrete  number  of  array  points: 


N 

Mean  =  W5  =   E   x(i)/N  r4  in 

i-1 


To  ensure  that  the  function  is  graphed  along  the  zero  axis 
of  the  formatted  graphic,  the  means  are  subtracted  from  the 
array  elements.  In  similar  fashion,  the  root  mean  square  is 
found: 
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N 


1N      2      1/2 
RMS  =  W6  =  [  Z   x(i)^/512]1/^  C4.12) 


The  maximum  and  minimum  values  for  the  input 
data  are  found  by  implementing  the  4052  call  "MAX"  and  "MIN" 
routines.  These  values  are  returned  in  W4  and  W3 
respectively.  The  values  are  then  output  to  the  axis 
defining  subroutines  for  establishment  of  appropriately 
scaled  and  labled  axis. 

Next,  the  3S  cursor  is  positioned  to  the 
beginning  of  the  array  with  a  "MOVE"  command.  For  the 
number  of  discrete  data  in  the  array,  the  function  is 
graphed  on  the  G5  CRT.  Finally,  the  graphic  is  titled  and 
identified  by  the  appropriate  parameters.  Figure  4.  10 
provides  a  general  flow  diagram  of  the  graphic  driver 
routines. 

The  particular  aspects  of  each  graph  will  now  be 
briefly  discussed. 

f.   Subroutine  GRIN 

This  subroutine  is  the  driving  graphic  for  the 

display  of  the   inputted   Direct,    Incident  and  Augmented 

Waveforms,  and  derived  Backscatter  waveform.    It  resides  in 

statements  5000-5410  and  makes-up  7%   of  the  program. 
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The  general  description  of  graphic  routines 
provided  above  is  most  typified  by  GRIN. 

Two  special  notes  need  to  be  made.  First, 
graphing  begins  with  the  ninth  element  of  the  array.  This 
is  to  prevent  the  display  of  an  erroneous  spike  input  by  the 
DPO.  (These  erroneous  values  have  a  minimal  effect  on  the 
processing).  Secondly,  the  function  is  graphed  in  discrete 
steps  equal  to  1/N  of  the  total  array  value  along  the  axis: 

Graphed  Point  =  (Q  *  Z1)  ,  X  (Q)         Cm.  il) 

where:     Q  =  array  number 
X (Q)  ■  array  element 
The  title   of   the  graph   is   selected  by   the 
routine  located  in  lines  5290-5330. 

The  horizontal  dimenion  for  the  graphic  is  time 
while  the  vertical  axis  is  scaled  in  volts, 
g.   Subroutine  GRFOUR 

This  subroutine  is  also  fairly  typical  of  the 
general  graphic  method.  It  displays  the  magnitude  and  phase 
components  of  the  spectral  data  in  polar  form  for  the  Direct 
and  3acXscattered  waveforms  and  for  the  impulse  response  of 
the  target. 
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Only  the  first  50  elements  of  the  magnitude 
array  are  graphed.  Values  beyond  the  first  60  are  of  a 
random  nature  due  to  the  processing  performed  to  derive  the 
transforms.  Spectral  content  much  beyond  10  GHZ  is, 
therefore,  negligible.  All  values  of  the  phase  array  are 
plotted,  although  the  note  concerning  the  randomness  of  the 
data  is  also  applicable  here. 

GRFOUR  is  actually  performed  in  two  distinct 
operations.  The  first  graphs  the  magnitude  of  the 
transform.  The  horizontal  axis  dimensions  are  converted  to 
frequency  from  time  by  finding  the  frequency  of  each  data 
point: 

Graphed  Point  =  Q/(Z1*257),  X  (Q)         t**.\H) 

The  second  operation  graphs  the  transform  phase. 
The  horizontal  axis  dimensions  are  converted  to  radians: 

Graphed  Point  =  (Q*2  7Y) /257,  X  (Q)         Ch.i^J 

GRFOUR  is  programmed  in  lines  5500-6080.     It  is  10%  of  the 
total  program. 

h.   GROP 

GROP  is  the  graphic  that  reproduces  the  shape 
representing  the  target.   That  is,  the  target  shape  based  on 
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physical   optics    principals   and   uncorrected    for    field   effects 
induced        by        residual        currents.  It        represents        a 

two-dimensional      cross-sectional      area   cut      of      the     target. 
Therefore,    only   target    length   and   height  are    shown. 

This  subroutine  is  constructed  generally  in  the 
same  format  as  the  other  graphic  subroutines.  One  important 
difference  is  that  a  true  one-to-one  relationship  must  be 
esxablished  for  the  coordinate  axis.  However,  before  this 
is  done,  the  metric  dimensions  of  the  target  length  must  be 
bounded.  This  is  done  in  lines  6210-6240.  Note  that  the 
array  BO,  source  for  the  shapes  storage  in  MATH,  has  already 
been  scaled  for  the  metric  dimension  in  the  y-direction  or 
amplitude.  The  x-dimension,  or  length,  is  scaled  for  metric 
measures  in  line  6220.  This  is  accomplished  by  multiplying 
the  time  between  samples  by  the  total  number  of  samples  in 
the  array  representing  the  target.  This  value  is  then 
compensated  by  the  averaging  antenna  parameter  established 
in    MATH,   S5: 


W2   -    (Z)(L)(1.2xl010US5)  (4>16) 


where:   Z  -  sampling  time  between  single  samples 
L  =  measured  length  in  number  of  samples 
S5  =  antenna  length  compensator 
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Q   =    low-pass    filter   compensator 
Q      is    always      one     for   zero     optimization,        or    no     low-pass 
filtering.  It      has      effect        only      when      optimization      is 

performed,  and  is  designed  to  compensate  for  the  spreading 
of  the  target  dimensions  in  the  x-direction  during 
filtering.  Note  also,  that  all  defined  targets  are  measured 
in  inches.  All  outputs  are  scaled  inches.  It  should  be  a 
relatively  simple  procedure  to  generalize  the  program  to 
handle   either   metric   or  english    units. 

The  next  major  difference  between  GEOP  and  other 
graphic  routines  is  that  all  graphing  is  done  in  GROP  itself 
without  calling  AXIS  or  any  other  external  formatting 
subroutine.  As  mentioned  above,  the  important  concept  is  to 
properly  relate  the  co-ordinant  axis  on  a  true  one-to-one 
plot.  This  is  done  is  statements  6290-6360  and  is  called 
"axis  scaling".  Axis  scaling  is  done  in  two  steps.  First 
the  general  x-axis  scaling  factor  A0  is  found.  This  is  done 
in  6290  by  dividing  the  maximum  horizontal  bounding  UDU  by 
the   viewport    established   by    the    display: 

A0    =     (W2    -    W1)     /    (V2-V1-5)  Ol.tl) 

This      yields      0DUfs/GD0      in   the      x-direction.  The      y-axis 

scaling     factor     A1      is     then   fixed      to      this      value.  The 
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dimensions  in  the  y-direction  are  then  checked  by  the 
program  in  lines  6320  to  ensure  that  "clipping"  or  windowing 
out  of  the  data  does  not  occur.  This  is  necessary  as  the 
vertical  viewing  axis  is  only  62£  that  of  the  x-axis.  If  the 
y-dimension  were  more  than  52%  of  the  x-dimension,  data 
would  be  lost  during  graphing.  To  prevent  this,  the  x-axis 
values  are  fixed  so  that  the  maximum  x-dimension  will  always 
be  plotted.  This  is  an  effective  reduction  in  the  size  of 
the  total  displayed  shape.  To  retain  the  true  one-to-one 
axis  relationship,  the  y-direction  dimensions  must  be 
reduced     by      a   corresponding      percentage.  This      reduction 

factor  is  computed  in  line  6330  and  applied  to  the  array  in 
6UU0.  Note  that  in  6330  twice  the  maximum  y-dimension  is 
used.  This  is  done  because  the  graphic  output  produces  the 
mirror  image  of  a  axisymmetric  object  while  the  array  itself 
is    single-sided. 

Once   scaling      is   established   in      statement   6380, 
the      array   is      plotted.  As  data      is     graphed      it   must      be 

dimensioned  in  english  units,  compensated  for  optimization 
effects,  and  compensated  for  any  viewport  clipping.  This  is 
done  on  a  point-by-point  basis  in  line  6410.  Note  that  the 
array  has  already  been  totally  dimensioned  and  compensated. 
It    is      in   fact    the      plotting   point      which   now    must     be    taken 
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into  consideration.  This  compensated  plotting  point  is 
returned  in  01  and  is  used  in  drawing  the  graph.  Finally, 
the  mirror  image  of  each  point  is  graphed  to  produce  a 
"true"  image.    A  flow  diagram  for  GROP  is   given  in  Figure 

4.  11. 

i.   Subroutine  GRASP 

GRAMP  graphs  the  time  domain  ramp  response  in 
lines   6600-6910.    It   comprises   5/5    of   the   program. 

GRAMP  is  identical  to  GRIN  in  practically  every 
detail.  Note  that  the  array  X  is  multiplied  by  1000  in  line 
6670.  This     converts     array        amplitudes      to      millivolts. 

GRIN   also   employs   this   conversion. 


160 


Ui 

Ss 

CD  * 

■4- a: 

C    X    UI 

o 

^    0.    3 

4^      ^      H-l 

iy  oc  oc 

CJ 

Ui 

X 

C. 

x§ 

Ui 

X 
•— •  Ui 

a: 

LOO 

a. 

o 

3 

k— < 

*r 

5p! 

-/7^ 

1 

Qi 

,J5 

Q. 

^1     •     1 

I 

1i 

a 

O 

O 

viz 

JT 

31 

03 

ID 

u 

3 

0 

CO 

0") 

7) 

/o  \ 

/    X 

/     *- 1 

^ 

r— T      <!-  ■> 

3 

o 

\     Q- 

\    *** 

/  o 

u. 

\  — '   i 

r    c 

\°  / 

• 

J- 

3 

tt 

-r4 

u. 

/ 

.      Ui 

( 
X      L 
CO     • 

o  c/> 

li     <-5 

(?) 

H-     2 
-     u 

!    <r 

•    a 

_J      < 

m    •; 

r    _ 

Zi 

\?J 

<x    : 

>w^ 

i— < 

H- 

H- 

n  i— 

ui  >- 

>■ 

■      f/) 

'/)      < 

D    Ui 

a: 

X      Ui 

c  a.    2- 

Ui      i 

a  <>} 

O    <C      «-t 

(.  a    a: 

v»-    i  j     <2, 

151 


Z-   EXPERIMENTAL  OBSERVATIONS 

In  this   chapter,    the   results   of   several  selected 

measurements  will  be  presented.    The  observations  made  were 

chosen  to   demonstrate  the  flexibility  of   the  Operating 

System   as  well   as   verify  experimental  results  with   the 
theoretical. 

A.   MEASUREHENT  1 — 7ARY  BORESISHTED  TARGET  DISTANCE 

1 .  Description  of  Se  tu£ 

The  purpose  of  this  measurement  is  to  observe  what 
effect  locating  the  target  a  distance  different  from  the 
transmitting  antenna  then  that  specified  for  the  receiving 
antenna,  but  maintaning  the  boresighted  relationship  of  the 
antennas  has  on  measurements  taken  by  the  OS.  The 
receiving  antenna  is  placed  2.54ft  from  the  transmitting 
antenna.  The  target  is  defined  from  the  target  library  on 
Tape  2. 

2 .  Results 

The  results  are  given  in  Figures  5.1  to  5.3  with  the 
backscattered  waveform,  ramp  response  and  physical  optics 
shapes  being  presented. 
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The  ramp  response  is  seen  to  be  typical  for  a 
cylinder  although  somewhat  reduced  in  size.  This  is  to  be 
expected  since  the  propagation  attenuation  losses  are 
greater  for  the  distance  at  which  the  target  is  located  in 
this  experiment  than  it  is  for  the  normal  control  case. 
With  the  distance  taken  into  consideration  by  the  Operating 
System  software,  it  can  be  seen  that  the  resultant  physical 
optics  shape  is  a  close  approximation.  The  target  width 
(length)  is  within  10$  of  the  true  length  and  the  diameter 
is  15,  US  of  the  true  diameter.  The  uncorrected  width  is  U9% 
greater  than  the  true  length.  The  backscattered  waveform 
shows  the  effect  of  some  spreading  in  comparison  to  the 
control  case  and  also  displays  a  reduced  S/N. 
3 •   Conclusions 

The  distance  between  the  receiving  antenna  and 
target  can  be  arbitrary  and  the  Operating  System  will  be 
able  to  provide  a  reasonably  accurate  interpretation  of  the 
target.  As  the  distance  increases,  this  accuracy  can  be 
expected  to  decrease  directly  to  the  limit  of  the  Systems 
ability  to  seperate  the  target  and  noise. 
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target:  run  date    dist   ant  tgt  remarks 

2  8-7-81   2.54m  i   1   location  difference 


MAXIMUM  PEAK  UALUE. 

MINIMUM  PEAK  UALUE 

RMS  UALUE 

MEAN  UALUE 

NUMBER  OF  UAUEFORMS  AUERAGED  ■  58 


♦1.68  nU 
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OPTIMIZATION  UALUE  -  8 


BACK  SCATTER  HAUE  FORM 


8. 38  1.88  1.38  2.88  2.58  3.88  3.58  4.88  4.58  5.88   nS 


Figure  5.1.   MEASUREMENT  1  —  Backscattered  Wavefora. 
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target:  run  date    dist 
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Figure    5.2.       MEASUREMENT    1—  Ramp   Response 
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TARGET:  RUN  DATE 


OIST   ANT  TGT  REMARKS 


2   8-7-81    2.54N   1 
DIAMETER  CINCHES-CORRECTED) 
WIDTH  < INCHES-CORRECTED) 
WIDTH  (IHCHES-UHCORRECTED)t* 
NUMBER  OF  WAUEFORMS  AUERAGED  -  58 


1    LOCATION  DIFFERENCE 
6.9625  (  6.9808) 
13.2847  (12.8888) 
17.3711  (12.0888) 

OPTIMIZATION  UALUE  - 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure  5.3.   MEASUREMENT  1—  Physical  Optics  Shape 
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B.   MEASUREMENT  2 — TARGET  AND  RECEIVING  ANTENNA  NOT 
30RESIGHTED 

1 .  Description  of  Setug 

An  attempt  is  made  to  observe  the  effect  of  placing 
a  target  at  an  arbitrary  location  on  the  TDRL  ground  plane. 
The  transmitted  field  is  not  axially  incident  on  the  target. 

For  this  experiment,  the  target  (cylinder  1)  was 
placed  2.11  meters  from  the  transmitting  antenna  and  at  a 
90°  angle  (approximately)  to  the  receiving/transmirting 
antenna  axis.  The  receiving  antenna  is  located  in  the 
control  area  of  the  TDRL.  The  target  is  at  about  a  60° 
angle  to  the  receiving  antenna  and  located  2.44  meters  from 
it.    The  experimental  setup  is  diagramed  in  Figure  5.4. 

2.  Results 

The  results  are  given  in  Figures  5.5  to  5.7 
representing  the  ramp  response,  physical  optics  shape  and 
backscattered  waveform.  Note  that  the  distance  input  to  the 
system  is  incorrect.  The  correct  value  as  read  from  the  DPO 
should  be  1.57  meters.  The  diameter  should  be  multiplied  by 
the  factor  1.53.  Doing  so  yields  a  diameter  of  6.3341 
inches  as  the  true  corrected  value.  The  resulting 
compensated  values  are  thus  very  close  to  the  targets  actual 
dimensions.   These  figures  may  be  suspect  as  subsequent  runs 
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RECEIVING 


Figure  5.4.   Measurement  2 — Physical  Setup 

never  yielded  the  same  degree   of  accuracy,    although  they 
approached  it.     The  important   item  is  that  the  physical 
optics  shape  is  a  close  approximation  of  that  expected  to  be 
representative  of  a  cylinder. 
3 .   Conclusions 

It  would  appear  that  the  targets  locatior  can  be 
generalized  if  the  shape  only  and  not  the  dimensions  is  the 
significant  item  of  interest.  However,  time  precluded 
further  verification  of  this  observation.  It  cannot  be  said 
with  absolute  certainty  that  this  observation  is  completely 
general.    Other   conf igurations  with  the  target   located  at 
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different  distances  from  the  boresighted  transmitting  and 
receiving  antennas  should  be  observed.  Note  that  for  proper 
compensation  for  prapagatior.  attenuation,  the  distance  the 
target  is  from  the  receiving  antenna  is  determined  by 
reading  the  DPO  Distance  Scale.  This  should  be  the  OS 
distance  input  value  vice  the  straight  line  measured 
distance  from  transmitting  antenna  to  target. 
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Figure   5.5.      Measurement  2~R»»P  Response 
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TfterFf  RUM  DATE     DtST   ANT  TGT  REMARKS 
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Figure  5.6.   Measurement  2 — Physical  Optics  Shape 
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Figure  5.7.   Measurement  2 — Backscattered  Waveform 
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C.   3EAS0REMENT  3  — ANTENNA  EFFECTS  ON  TARGET  RESPONSES 

1 .  2§scrip_t ion  of  Setu2 

Th-e  purpose  of  this  measurement  is  to  verify  that 
the  effects  of  the  antenna  are  in  fact  deconvol/ed  out 
during  processing.  Sost  previous  runs  have  been  performed 
using  antenna  1  (TSM  Horn)  .  This  measurement  is  performed 
using  three  different  length  monopole  antennas.  The 
antennas  used  were  antenna  2  (see  TABLE  II)  (85  mm 
monopole) ,  antenna  8  (205  mm  monopole)  and  antenna  12  (285 
mm  monopole) .  The  target  used  was  target  1  (see  TA3LE  I) . 
Both  target  and  antenna  were  placed  in  the  control 
locations,  1.27  meters  from  the  transmitting  antenna. 
21-acguisitions  were  averaged  for  each  run.  The  results 
obtained  should  be  compared  with  those  obtained  previously 
for  the  TEM  horn. 

2.  Results 

The  results  are  presented  in  Figures  5.7  to  5.24 
Each  antenna  has  the  results  of  its  run  shown  in  six  graphs. 
The  first  graph  depicts  the  direct  waveform  as  seen  by  the 
individual  antenna,  the  second  is  the  ramp  response 
resulting,  the  third  figure  is  the  physical  optics 
representation  of  the  target  shape,  and  the  final  three 
figures  are  the  optimized   frequency  domain   impulse,   time 
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domain  impulse  and  ramp  response  for  the  target.  Figure  5.7 
to  5.12  are  for  monopole  antenna  2.  Figures  5.13  to  5.18 
are  for  monopole  antenna  8  and  figures  5.19  to  5.24  are  for 
monopole  antenna  12. 

The  first  obvious  result  is  that  in  comparison  of 
the  ramp  responses,  the  smaller  ihe  antenna  is  in  relation 
to  the  target,  the  less  typical  is  the  ramp  and, 
correspondingly,  the  physical  optics  shape.  However,  for 
antennas  considerably  taller  than  the  height  of  the  target, 
a  significant  amount  of  noise  pollutes  the  ramp  response. 
In  comparing  the  target  impulse  responses,  it  can  be  seen 
that  the  exact  opposite  appears  to  be  true.  The  most 
representative  impulse  is  for  the  the  shortest  monopole. 
The  frequency  domain  impulse  graphs  show  that  the  longer 
monopoles  accept  a  higer  frequency  content  than  the  shortest 
monopole. 

3 .   Conclusions 

First,  the  closer  the  antenna  is  to  natural 
resonance  with  the  incident  pulse,  the  more  faithful  appears 
the  resulting  impulse  response  to  theory.  But  it  does  not 
naturally  follow  that  the  corresponding  ramp  response  will 
be  as  true  a  repsentat ion.  The  opposite  appears  to  be  true. 
It  may  be  that  the  short  antenna  is  shadowed  by  a  portion  of 
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the  target.  More  likely,  the  short  antenna  is  just  not 
large  enough  to  accept  the  full  radiated  backscatter  of  the 
target.  Conversely,  the  235-mm  monopole  is  sensitive  to 
environmental  noise.  As  can  be  seen  in  figures  5.20  and 
5.14,  the  ramp  responses  display  the  relatively  greater 
sensitivity  of  the  285-mm  monopole  to  environmental  noise 
than  the  205-mm  monopole.  Although  the  observation  is 
limited  to  just  these  runs,  it  may  be  possibly  concluded 
that  there  is  a  "best"  monopole  antenna  length  for  a 
particular  target  size. 

In  conclusion,  it  can  be  said  that  the  effects  of 
the  antenna  are  in  general  deconvolved  out.  However,  other 
effects  indirectly  related  to  the  target/antenna 
configuration  may  introduce  factors  which  materially  effect 
the  reproduction  of  the  target  shape  and  its  overall 
responses. 
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Figure   5.8.      Measurement   3 — Direct   Waveform      (85   mm) 
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TARGET:  RUN  DATE     OIST   ANT  TCT  REMARKS 
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Figure  5,9.   Measurement  3 — Ramp  Response  (85  mn) 
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TARGET:  RUN  DATE     DIST 
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Figure  5.10.   Measurement  3  — Physical  Optics  Shape  (85  ma) 
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TARGET:  RUN  DATE 

13  8-11-91 


DIST   ANT  TGT   REMARKS 

1.27M   2    1    DIPOLE  ANTENNA  #2<85MM> 


MAXIMUM  PEAK  UALUE 8.9099  hU 

MINIMUM  PEAK  UALUE -9.999?  nU 

RMS  UALUE 9.  9992  hU 

MEAN  UALUE 9.9989  nU 

NUMBER  OF  WAVEFORMS  AUERAGED  *  21      OPTIMIZATION  UALUE  ■  8.3 


flU 


9.8985 


8.8888 


-8.8885 


IMPULSE  RESPONSE 


/ 

/\ 

V 

y 

\, 

^ 

\ .  y 

'Y 

J 

\^ 

^* 

**  II  II 

■~s* 

\ 

/ 

V 

H 

V 

8. 58  1.88  1.38  2.88  2.58  3.88  3.58  4.08  4.58  5.88   nS 


Figure  5.  II, 


Measurement  3 — Time-Domain  Impulse  Response 
(85  mm) 
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C<*'' 


TARGET:  RUN  DATE     DIST   ANT  TGT   REMARKS 

13  8-11-81   1.27*1  2    1    DIPOLE  ANTENNA  #2<83MM> 
NUMBER  OF  WAUEFORMS  AUERAGED  »  21      OPTIMIZATION  UALUE  ■  3.5 


IMPULSE  FFT— MAGNITUDE 


1   A 

3.340 

i 

■iM 

8.920 

i 

/  ^ 

9.9*8 

) 

0.89    2.89    4.89    6.89    3.89   19.99   12.99    GHZ 


IMPULSE  FFT  — PHASE 


8.3   8.3    1.9    1.5    2.8    2.5    3.8   3.5    4.9    4.5    5.8   5.5    6.8        RAO 


Figure  5.12, 


Measurement   3 — Frequency-Domain   Impulse   Res- 
ponse   (8  5   mm) 
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TARGET:  RUN  DATE 

13  3-11-81 


DIST   ANT  TGT  REMARKS 

1.27M  2   1    OIPOLE  ANTENNA  »2<85MM> 


MAXIMUM  PEAK  UALUE +613.71  nU-M 

MINIMUM  PEAK  UALUE -345.20  nU-M 

RMS  UALUE +284.65  iiU-M 

MEAN  UALUE +141.79  nU-M 

NUMBER  OF  WAUEFORMS  AUERAGED  =  21      OPTIMIZATION  UALUE  ■  8.5 


nU-M 


RAMP  RESPONSE  —  TIME 


690 

488 

• 

288 

9 

-  ^ 

"N    | 

?00 

\ 

/ 

/ 

V 

J 

8.58  1.08  1.59  2.89  2.59  3.99  3.59  4.99  4.59  5.99   nS 


Figure  5.13.   Measurement  3 — Ramp  Optimized  (85  mm) 
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TARGET:  RUN  DATE     DIST   ANT  TGT 
14  8-11-81   1.27H  8   1 


REMARKS 

DIPOLE  ANTENNA  18  <283NH> 


MAXIMUM  PEAK  UALUE +184.92  mU 

MINIMUM  PEAK  UALUE -116.61  nU 

RMS  UALUE +59.38  mU 

MEAN  UALUE 8.  09  mU 

NUHBER  OF  WAUEFORMS  AUERAGED  »  21      OPTIMIZATION  UALUE 


nU 


OIRECT  WAUEFORM 


188 

58 

■/ 

~\ 

8 

•/ 

\    , 

-58 

188 

9.58  1.88  1.58  2.88  2.58  3.88  3.58  4.88  4.58  5.89   nS 


Figure  5.  14.      Measurement  3  —  Direct    Waveform    (205   am) 
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target:  run  pate    dist   ant  tct  remarks 

14  8-11-91  1.27m  8  1   dipole  antenna  »8  (285mm) 

maximum  peak  ualue +5.33  nu-m 

minimum  peak  ualue -689.89  nu-m 

RMS  UALUE *357.89  nU-M 

MEAN  UALUE -299.  54  nU-M 

NUMBER  OF  WAUEFORMS  AVERAGED  =  21       OPTIMIZATION  UALUE  *  8 


nU-M 


RAMP  RESPONSE  —  TIME 


-288 


-488 


-688 


- 

\y 

\ 

• 

\ 

• 

V 

m 

8. 58  1.88  1.58  2.88  2.58  3.88  3.58  4.88  4.58  5.88   nS 


Fiqure  5.15.   Measure.ent  3-Ha-P  Hesponse  (205  ..) 
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9 


TARGET:  RUN  OATE     DIST 

14  3-11-81   1.27M 
DIAMETER  (INCHES-CORRECTED) 
WIDTH  (INCHES-CORRECTED) 
WIDTH  (INCHES-UNCORRECTED)** 
NUMBER  OF  WAVEFORMS  AVERAGED  ■  21 


1 


ANT  TGT   REMARKS 

DIPOLE  ANTENNA  f9  <293MM> 
3.3341  (  6.8088) 
12.8734  (12.8098) 
19.3945  (12.8888) 

OPTIMIZATION  VALUE 


8 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure  5.16.   Measurement  3 — Physical  Optics  Shape  (205  am) 
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TARGET:  RUN  DATE     DIST   ANT  TGT  REMARKS 

14  8-11-31   1.27M  8   1    DIPOLE  ANTENNA  «8  <285Mf1> 

MAXIMUM  PEAK  UALUE 8.0086  nU 

MINIMUM  PEAK  UALUE -8.0886  nU 

RHS  UALUE 8.8082  nU 

MEAN  UALUE 8.8000  hU 

NUMBER  OF  WAUEFORMS  AUERAGED  =    21      OPTIMIZATION  UALUE  -  8.5 


«U 


8.8885 


8.8888 


-8.8885 


IMPULSE  RESPONSE 


■ 

rv 

A 

r 

V 

/T\ 

./ 

<"> 

\ 

\ 

/'  J 

\7 

\ 

i 

8.58    1.88    1.58    2.88    2.58    3.88    3.58    4.88    4.58    5.88       nS 


Pigure  5.17. 


Measurement   3 — Time-Domain  Impulse   Response 
(205  mm) 
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TARGET:  RUN  DATE     OIST   ANT  TCT   REMARKS 

14  8-11-81   1.27M   8    1     DIPOLE  ANTENNA  #8  <285MM> 
NUNBER  OF  UAUEFORMS  AUERAGED  •  21      OPTIMISATION  UALUE  -  8.3 


IMPULSE   FFT 

—MAGNITUDE 

k 

8.848 

) 

■  fki 

L.    -  — .     ... , 

8.828 

) 

yiVV 

F1 

8.888 

r~ 

8.88    2.88    4.88    6.88    8.88   18.88   12.88    GHZ 


IMPULSE  FFT— PHASE 


8.8    8.S    1.8    1.3    2.8    2.3    3.8    3.3    4.8   4.5    3.8    5.3    6.8        RAD 


Figure   5.18. 


Heasurement    3 — Frequency-Domain   Impulse   Res- 
ponse   (2  05   mm) 


186 


TARGET:  RUN  DATE     DIST   ANT  TGT  REMARKS 

14  8-11-81   1.27M  8   1    DIPOLE  ANTENNA  18  <283MM> 

MAXIMUM  PEAK  UALUE *8.  19  nU-M 

MINIMUM  PEAK  UALUE -537.10  nU-M 

RMS  UALUE +318.08  «U-M 

MEAN  UALUE -255.31  mU-M 

NUMBER  OF  WAUEFORMS  AUERAGED  =  21  OPTIMIZATION  UALUE  =  8.5 

RAMP  RESPONSE  —  TIME 
nU-M 


8 

'N 

280 

400 

- 

0. 53  1.80  1.50  2.00  2.50  3.00  3.30  4.00  4.50  5.00   nS 


Figure  5. 19.   Measurement  3 — Ramp  Optimized   (205  mm) 
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TARGET:  RUN  DATE 

12  8-U-91 


OIST   ANT  TGT  REMARKS 

1.27M   12   1    DIPOLE  ANTENNA  il2<285MH> 


MAXIMUM  PEAK  UALUE +101.31  nU 

MINIMUM  PEAK  UALUE -122.78  «U 

RMS  UALUE +67.  32  fiU 

MEAN  UALUE 9.38  mU 

NUMBER  OF  WAUEFORNS  AUERAGED  =    21      OPTIMIZATION  UALUE 


nU 

188 
58 

8 

-58 

-188 


DIRECT  WAUEFORM 


- 

•/ 

•/ 

- 

8.59  1.89  1.59  2.99  2.59  3.98  3.59  4.99  4.59  5.99   nS 


Figure  5.20.   Measurement  3  — Direct  Waveform  (285  mm) 
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target:  run  date    dist   ant  tgt  remarks 

12  8-11-31  1.27n  12  1   dipole  antenna  #12<283hm> 

maximum  peak  ualue +139.35  nu-m 

minimum  peak  ualue -563.13  nv-m 

RMS  UALUE *328.  92  nU-M 

MEAN  UALUE -252.  72  nU-M 

NUMBER  OF  WAVEFORMS  AUERAGED  =  21      OPTIMIZATION  UALUE  *  8 


nU-H 

J\  A 

RAMP 

RESPOI 

HSE   — 

•    TIME 

8 
-208 

£\ 

i 

S\ 

y 

J^V 

f 

t 

V 

\r 

-488 

\„A 

Y 

8.58  1.88  1.58  2.88  2.58  3.88  3.58  4.88  4.58  5.88   nS 


Figure  5.21.   Heasurement  3— Samp  Response  (285  im) 
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target:  run  date 


DIST   ANT  TGT   REMARKS 


12  8-11-31   1.27M   12   1 
DIAMETER  (INCHES-CORRECTED) 

WIDTH  (INCHES-CORRECTED' 
UIDTH  (INCHES-UNCORRECTED)t* 
NUMBER  OF  UAUEFORMS  AUERAGED  »  21 


DIPOLE  ANTENNA  #12<283HH> 
6.8809  C  6.9889) 
12.8989  (12.9889) 
16.2395  (12.9998) 

OPTIMIZATION  UALUE  -  9 


RAMP  RESPOHSE  —  PHYSICAL  OPTICS 


Pigure  5.22.   Measurement  3 — Physical  Optics  Shape  (285  mm) 
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(r* 


TARGET:  RUN  PATE 

12  8-11-81 


DIST   ANT  TGT  REMARKS 

1.27M   12   1    DIPOLE  ANTENNA  #12<285HH> 


MAXIMUM  PEAK  UALUE 9.0906  «U 

MINIMUN  PEAK  UALUE -9.9893  nU 

RMS  UALUE 9.9992  r«U 


MEAN  UALUE 

NUMBER  OF  WAVEFORMS  AUERAGED 


21 


9.9999  mU 

OPTIMIZATION  UALUE  «  9.5 


nV 


9.9995 


0.9099 


-0.0995 


IMPULSE  RESPONSE 


- 

A 

\  . 

A 

^r 

\ ,/ 

\l 

v 

0.50    1.00    1.58    2.00    2.50    3.00    3.50    4.09    4.50    5.08       nS 


Figure   5.23. 


Measurement    3— Time-Domain   Impulse   Response 
(285  mm) 
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TARGET:  RUN  PATE     DIST    ANT  TGT   REMARKS 

12  8-11-81   1.2711   12   1     DIPOLE  ANTENNA  il2(285MM> 
NUNBER  OF  UAUEFORMS  AUERAGED  »  21      OPTIMIZATION  UALUE  •  8.3 


IMPULSE  FFT— MAGNITUDE 


8.828  3 
8.888  3 


■M 

3 

rl.\ 

Vw  . 

3 

8.88    2.88    4.88    6.88    8.88   18.88   12.88    GHZ 


IMPULSE  FFT— PHASE 


5.88 

8.88 

-5.88 


;fc 


^ 


^vfe 


8.8    8.3    1.8    1.3    2.8    2.5   3.8   3.3    4.8    4.5    5.8    5.5   6.8        RAD 


Figure   5. 24. 


Measurement    3 — Frequency-Domain    Impulse   Res- 
ponse   (2  85   mm) 
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TARGET:  RUN  DATE 

12  8-11-81 


DIST   ANT  TGT   REMARKS 

1.27M   12   1     DIPOLE  ANTENNA  il2<285MM> 


MAXIMUM  PEAK  UALUE +28.80  nU-N 

M1HIMUM  PEAK  UALUE -596.67  *U-M 

RMS  UALUE +368.85  *U-M 

MEAN  UALUE -314.  14  nU-M 

NUMBER  OF  WAUEFORMS  AUERAGED  «  21      OPTIMIZATION  UALUE  -  8.3 


nU-M 


RAMP  RESPONSE  —  TIME 


-200. 


-408 


0.50  1.80  1.30  2.80  2.50  3.00  3.50  4.00  4.50  5.80   nS 


Figure  5.25.   Seasurement  3 — Ramp  Optimized  (285  am) 
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D.   MEASUREMENT  4  —  ARBITRARY  COMPOSITE  TARGET 
1 •   Description  of  Setup 

This  measurement  was  made  to  observe  the  response  of 
the  OS  to  a  completely  arbitrary,  but  axisymmetric  target 
not  encountered  by  the  system  previously.  The  composite 
target  was  a  combination  of  target  11,  a  cone,  and  target  3, 
a  cylinder.  A  drawing  of  the  composite  is  given  in  Figure 
5.25. 


Pigure  5.26.   Measurement  4 — Composite  Target 

2.   Results 

The  results  are  given  in  Figures  5.26  and  5.27,  the 

physical  optics  shape  and  ramp  responses  respectively.  The 

corrected   diameter  is   11%   larger   than  actual,    and  the 

corrected  length   is   35%   shorter   than    actual.  The 
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uncorrected  length  is  13%  greater  than  actual.  The  shape  is 
a  close  approximation  of  the  targets  actual  shape,  although 
somewhat  shortened.  The  cone  portion  appears  to  be  about 
one-half  of  the  the  total  figure,  as  it  was  in  the  actual 
target. 

3 .   Conclusions 

Based  on  this  observation  and  the  many  others  taken 
during  the  implementation  of  the  OS,  it  can  be  said  that  the 
OS  can  be  relyed  upon  to  produce  recognizable  representative 
shapes  for  known  or  unknown  axisymmetric  metallic  targets. 
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TARGET:  RUN  DATE 


DIST    ANT  TGT   REMARKS 


2?  9-11-81   1.27N   1    ???   COMPOSITE  ARBITRARY  TARGET 


DIAMETER  (INCHES-CORRECTED) 
UIDTH  (INCHES-CORRECTED) 
WIDTH  (INCHES-UNCORRECTED)t* 
UMBER  OF  WAUEFORMS  AUERAGED  »  21 


6.5659  (  9.9908) 

8.9994  (  0.9998) 

13.5352  (  9.9988) 

OPTIMIZATION  UALUE 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure  5.27.   Measurement  *— Physical  Optics  Shape 
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TARGET:  RUN  DATE     DIST   ANT  TCT  REMARKS 

27  3-11-91   1.27M  1  ???  COMPOSITE  ARBITRARY  TARGET 

MAXIMUM  PEAK  UALUE +38.86  nU-M         j 

MINIMUM  PEAK  UALUE -33S.77  mU-H      ^ 

RMS  UALUE +171.95  mU-M 

MEAN  UALUE -115.28  nU-M 

NUMBER  OF  WAUEFORMS  AVERAGED  ■  21      OPTIMIZATION  UALUE  *  3 


mU-M 

e 

-53 
-130 
-153 
-293 
-253 
-389 


RAMP  RESPONSE  —  TIME 


\^        ,\ , , , , , , . 

: ^"^zj- 


3. 59    1.33    1.S8    2.38    2.58    3.88    3.58    4.88    4.50    5.88       nS 


Figure   5.28.       Keasur anient   <* —  Ramp  Response 


197 


E.   MEASUREMENT  5  — DIELECTRIC  CYLINDERS 

1 .  Description  of  Setup 

Two  runs  were  made  using  different  dielectric 
cylinders  in  place  of  metallic  objects  as  targets.  The 
dielectric  cylinders  used  are  similar  in  configuration  to 
cylinder  targets  2  and  3  machined  from  aluminum  stock,  but 
are  made  of  plastic.  For  this  measurement,  the  dielectric 
cylinders  will  refered  to  as  target  2  and  3,  respectively. 

The  purpose  of  the  measurement  was  to  observe  the 
response  of  the  dielectric  cylinders  to  axially  incident 
electromagnetic  fields  and  to  also  see  how  the  OS  responds 
to  a  target  which  was  non-metallic.  The  targets  were  placed 
in  the  standard  control  location,  1.27  meters  from  the 
transmitting  antenna.  The  runs  assumed  the  targets  were  of 
arbitrary,  unknown  shape  for  3S  observation. 

2.  Results 

Target  2  was  run  first.  The  results  are  given  in 
Figures  5.28  to  5.31.  The  physical  optics  shape  in  Figure 
5.23  is  a  very  poor  representation  of  the  true  targets 
shape.  The  reason  for  this  can  be  observed  in  the  ramp 
response  of  Figure  5.29.  The  S/N  is  extremely  low.  A  noise 
spike  at  about  1.6  nsecs  has  caused  the  entire  resultant 
shape  to  be  foreshortened.   The  amplitude  discrepancy  can  be 
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attributed  to  the  fact  that  the  amplitude  compensator  was 
developed  for  the  relatively  large  valued  returns  of  the 
metallic  cylinders.  nevertheless,  when  the  ramp  response  is 
optimized,  the  general  shape  of  the  target  can  be  seen  as  in 
Figure  5.30  and  5.31. 

The  results  using  Target  3  are  shown  in  Figures  5.32 
to  5.35.  Because  this  target  is  physically  larger  than 
target  2,  it  has  returned  a  substantially  larger  signal, 
large  enough  so  that  system  noise  is  less  of  a  significant 
factor  (although  still  seen  on  the  waveform)  .  The  amplitude 
is  again  much  below  that  to  be  expected  for  a  metallic 
target  of  similiar  size. 
3 .   Conclusions 

Although  the  information  provided  by  this 
measurement  is  limited,  it  appears  that  dielectric  targets 
might  provide  returns  which  can  be  identified  by  time  domain 
technigues.  This  could  be  an  important  area  for  more 
intensive  research  as  the  use  of  high  density  polymers  and 
resins  in  real  targets  of  interest  is  of  growing  importance 
to  surveillance  systems  attempting  to  locate  and  identify 
such  targets. 
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TARGET:  RUN  DATE 


DIST   ANT  TCT   REMARKS 


1 


24  9-11-91   1.27H 
DIAMETER  (INCHES-CORRECTED) 
WIDTH  (INCHES-CORRECTED) 
WIDTH  (INCHES-UNCORRECTED^t* 
NUMBER  OF  WAUEFORMS  AUERAGED  -  21 


???  DIELECTRIC  CYLINDER  <3X12  IM) 
0.7662  (  8. 8000) 
2.1596  (  8.8000) 
3.2813  (  8.0000) 

OPTIMIZATION  UALUE  -  8 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


\  \l\^NSNs|vf*i«^#-**- 


Figure  5.29.      Measurement   5 — Physical  Optics    (Target    2) 
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TARGET:  RUN  DATE 

24  8-11-31 


DIST   ANT  TCT   REMARKS 

1.27H  1    ???  DIELECTRIC  CYLINDER  <3X12  IN) 


MAXIMUM  PEAK  VALUE 

MINIMUM  PEAK  UALUE 

RMS  UALUE , 

MEAN  UALUE 

NUMBER  OF  UAUEFORMS  AUEPAGED 


21 


♦11.15  nU-M 
-38.59  nU-M 
♦17.97  nU-M 
-13.94  mU-M 

OPTIMIZATION  UALUE 


RAMP  RESPONSE  —  TIME 


3.30  1.80  1.58  2.88  2.58  3.88  3.58  4.88  4.59  5.88   nS 


Figure  5.30.   Measarenent  5  —  Ramp  Response  (Target  2) 
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TARGET:  RUN  DATE 


24  9-11-81   1.27M   1 
DIAMETER  (INCHES-CORRECTED) 
UIDTH  CINCHES-CORRECTED) 
WIDTH  <INCHES-UNCORRECTED)t* 
NUMBER  OF  WAUEFORMS  AUERAGED  -  21 


DIST   ANT  TGT  REMARKS 


'??  DIELECTRIC  CYLINDER  (3X12  IN) 
9.7662  <  8.0890) 
2.1596  (  8.0988) 
9.8828  C  8.9898) 

OPTIMIZATION  UALUE  -  1 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


Figure   5.31. 


Measurement    5 — Optimized   Physical   Optics    (Tar- 
get  2) 
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TARGET:  RUN  DATE 


24  8-11-81   1.27M   1 

MAXIMUM  PEAK  UALUE 

MINIMUM  PEAK  UALUE 

RMS  UALUE 

MEAN  UALUE 

NUMBER  OF  WAUEFORMS  AUERAGED  =  21 


OIST   ANT  TGT  REMARKS 

???   DIELECTRIC  CYLINDER  <3X12  IH> 


+8.29  nU-M 
-27.71  nU-M 
♦  13.28  riU-M 
-10.39   nU-M 

OPTIMIZATION  UALUE 


RAMP  RESPONSE  —  TIME 


8.38  1.88  1.38  2.88  2.58  3.88  3.S8  4.88  4.38  5.88   nS 


£ 


Figure  5.32.      Measurement   5--Optimized   Ramp    (Target    2) 


203 


TARGET:  RUN  DATE 


DIST   ANT  TGT  REMARKS 


23  8-11-31   1.27M   1 
DIAMETER  (INCHES-CORRECTED) 
WIDTH  (INCHES-CORRECTED) 
WIDTH  vINCHES-UNCORRECTED>t* 
NUMBER  OF  WAUEFORMS  AUERAGED  "  21 


???  DIELECTRIC  CYLINDER  <6X6  IH> 
1.4154  (  0.0808) 
5.7075  <  8.0888) 
3.571?  <  8.0880) 

OPTIMIZATION  UALUE  ■  8 


RAMP  RESPONSE  —  PHYSICAL  OPTICS 


NiMgilgipil*"**** 


5 


Figure  5.33.   Measurement  5--Physical  Optics  (Target  3) 
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TARGET:  RUM  DATE     01  ST   ANT  TGT   REMARKS 

23  8-11-81   1.27M   1    ???   DIELECTRIC  CYLINDER  (.6X6   IN> 


MAXIMUM  PEAK  UALUE. 
MINIMUM  PEAK  VALUE. 

RMS  UALUE 

MEAN  UALUE 


NUMBER  OF  WAUEFORMS  AUERAGED  »  21 


+94.83  tiU-M 

-71.29  nU-M 

♦33.88  mU-M 

-3.59  nU-M 

OPTIMIZATION  UALUE  •  9 


nU-H 


RAMP    RESPONSE   —    TIME 


<£w 


vVyM 


«V^ 


9.S9    1.08    1.59    2.00    2.58    3.89    3.38    4.08    4.50    3.99       nS 


•  % 


Figure  5.34.   Measurement  5 — Ramp  Response  (Target  3) 
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target:  run  date    dist  ant  tgt  remarks 

23  8-11-81  1.27m  1  ???  dielectric  cylinder  <6x6  in> 

maximum  peak  ualue 8.0903  nu 

minimum  peak  ualue -9.0884  nu 

RMS  UALUE 0.0881  nU 

MEAN  UALUE 0.  0000  mU 

NUMBER  OF  WAUEFORHS  AUERAGED  «  21      OPTIMIZATION  UALUE  «  8.3 


fiU 


8.8882 


8.8088 


-0.8002 


IMPULSE  RESPONSE 


8.30    1.80    1.30    2.88    2.38    3.80    3.39    4.80    4.30    3.80       nS 


Figure   5.35. 


Measureient    5 — Optimized   Physical    Optics    (Tar- 
get   3) 
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TARGET:  RUN  DATE     OIST   ANT  TCT   REMARKS 

25  8-11-81   1.2711   1    ???   DIELECTRIC  CYLINDER  (6X6  IN) 
NUMBER  OF  UAUEFORMS  AUERAGED  ■  21       OPTIMIZATION  UALUE  ■  0.3 


IMPULSE  FFT— MAGNITUDE 


8.828  3 
8.8153 
8.8183 
8.885  3" 
8.880  3 


8.88 


:.88    4.88    6.88    3.88   18.88   12.88    GHZ 


3.88 

IMPULSE 

FFT--PHASE 

0.80 

1'tiiU 

K 

\iW 

mM 

A  ihkf^AH^  r 

iiuii^iL 

3.88 

•Pr     It" 

Upi 

VT 

fWv  \Ji    yy 

fylWj^'   r 

8.9  9.5  1.8  1.5  2.8  2.5  3.8  3.5  4.8  4.5  5.8  5.5  6.0    RAD 


*    •   % 


Cif 


Figure  5.36.   Measurement  5 — Optimized  Ramp  (Target  3) 
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VI.   SUMMARY 

A  software  Operating  System  has  beer,  developed  for  the 
Time  Domain  Radar  Laboratory.  The  System  has  been  designed 
to  be  relatively  easy  for  inexperienced  operators  to  use 
after  an  initial  introduction  -co  the  OS  code  and  -he 
constituent  hardware  controlled.  The  OS  extends  the 
capaoility  of  the  TDRL  to  the  point  where  accurately 
processed  transient  time-domain  data  is  presented  in  a  form 
suitable  for  operator  evaluation.  The  OS  should  serve  as  a 
primay  research  tool  for  the  study  of  electromagnetic 
transients  at  the  Naval  Postgraduate  School. 

The  OS  is  highly  operator-processor  interactive. 
Although  the  code  that  makes-up  the  05  is  relatively 
complicated,  its  actual  implementation  and  use  is  not 
difficult.  There  are  numerous  directives  issued  by  the 
System  and  displayed  on  the  Graphic  Systems  Crt  to  the 
operator  requesting  interactively  the  input  of  information 
through  the  GS  keyboard  at  vital  points  in  the  programs. 

The  OS  is  composed  of  three  subprograms.   They  are: 
1.   INPUT  for   the  inputting   of  signal   information  from 
the  TDRL  ground  plane 
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2.  MATH  for  performing  the  significant  mathematical 
computations  for  the  05 

3.  GRAPH  for  the  presentation  of  the  processed 
information  in  a  form  suitable  for  user  evaluation 
and  analysis 

INPUT  receives  tha  desired  signal  in.  discrete  digitized 
form  from  the  systems  Digital  Processing  Oscilloscope,  the 
principal  analog  to  digital  converter  use!  in  the  TDRL.  The 
source  signal  is  a  rapid  rise  time  gaussian  pulse  simulating 
an  impulse.  Tha  acquired  information  is  averaged  a  number 
of  times  as  interactively  requested  by  tha  operator.  This 
is  done  to  reduce  the  white  gaussian  noise  content  of  the 
acquired  analog  signal. 

3ATH  performs  the  Fouriar  conversion  and  signal 
decon volution  to  ietermine  the  impulse  response  of  the 
target.  It  also  implements  a  noise  reduction  algorithm  -co 
optimize  tha  processed  data.  The  amount  of  optimization  is 
under  the  control  of  the  operator  and  is  dependant  on  the 
degree  to  which  loss  of  small  surface  features  can  be 
tolerated  and  yet  obtain  a  good  representation  of  the 
targets  true  shape  and  dimensions.  The  step  and  ramp 
responses  are  also  obtained  and  stored  for  later  evaluation. 
HATH   determines  the   Physical  Optics   Shape   of  the   target 
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based  on  the  ramp  response  uncorrected  for  the  effects  of 
induced   currents      on    the      boiy    of      the   target.  It    further 

computes  and  applies  compensation  factors  for  length  and 
width  based  on  historically  observei  results  from 
pre-defined  targets.  The  compensators  are  applied  directly 
to  unknown  targets  resulting  in  an  accuracy  of  10%  or  better 
in    determination   of   the  true   target    dimensions. 

GRAPH  outputs  the  processed  data  to  the  GS  CRT  of  the 
4052  Graphic  System.  Information  graphed  includes  the  real 
time        direct,         incident,  augmented      and        backscattered 

waveforms,  the  time-dcmain  ramp  response,  physical  optics 
shape  of  the  target  and  the  frequency  domain  spectral  data 
for  the  direct  and  backscattered  waveforms  as  well  as  the 
frequency  domain  impulse  response.  A  means  for  inputting 
the  optimization  variable  and  for  the  return  of  program 
control   to    INPOT   to    acquire   a    new    signal   is    also    provided. 

The  total  result  is  a  highly  versatile  and  powerful 
transient      analysis   tool.  The    OS      has     flexiblity    in      its 

internal  organization  to  allow  for  future  expansion  and/or 
modification  so  as  to  allow  the  control  of  constituent 
hardware  and  software  in  applications  involving  other 
aspects     of      time-domain   target      identification.  This      is 

demonstrated   in    the   complete    integration   of    a    Prony*3    Method 
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Hybrid  Solution  foe  the  location  of  a  targets  poles  and 
residues  with  the  requirements  for  pure  transient 
time-domain  analysis  procedures. 

The  code  itself  is  complete  as  it  stands.  However,  there 
are  areas  in  which  it  could  oe  made  mere  efficient. 
Particularly,  a  number  of  graphic  subroutines  could  be 
combined  into  one  routina  so  as  to  reduce  memory 
requirements.  Also,  a  method  for  the  display  and  selection 
of  menu  items  using  placement,  of  the  scraen  cursor  to  select 
the  item  might  be  developed.  This  would  allow  for  the 
display  of  more  information  graphing  processes  than  is  new 
allowed. 

One  of  the  specialized  algorithms  developed  for  the  OS 
is  a  method  for  expanding  certain  displayed  graphs  while 
retaining  a  true  plot  of  relevant  information.  The  result 
is  that  the  GS  has  essentially  infinite  resolution  and  would 
allow  for  a  close  look  at  particular  areas  of  the  display 
without  disturbing  the  true  graph  relationships.  Although 
presently  applied  to  only  ona  subroutine  in  GRAPH,  it  could 
easily  be  adapted  to  any  other  graphic  subroutine. 

A  second  specialized  routine  is  one  whereby  a  true 
one-to-one  relationship  between  coordinate  axis  is 
maintained.    The   unique  feature  of   this  algorithm   is  its 
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generality  in  being  applicable  -co  any  size  lata  field  on  any 
size  viewport  without  ever  losing  information  due  to 
"viewport  clipping". 

The  results  of  field  observations  using  the  OS  have 
proven  remarkably  faithful  to  theoretical  rasults.  However, 
the  present  use  of  the  OS  has  been  principally  directed 
towards  verifying  that  the  system  was  up  and  working.  As 
such,  data  was  collected  for  axisymmetrio  targets  almost 
exclusively  as  there  was  a  large  amount  of  literary  support 
in  this  research  area.  There  still  is  a  need  to  explore 
system  operation  with  this  configuration  of  targets.  In 
particular,  the  OS  needs  to  be  extended  so  as  to  be  able  to 
solve  the  Magnetic  Field  Integral  Equation  for  EM  transients 
on  the  surface  of  a  target.  This  calculation  would  then  be 
compared  with  field  results  to  obtain  a  corrected  and 
accurate  representation  of  the  target  both  in  shape  and 
dimension. 

Hardware  has  proven  very  well  suited  to  the  needs  of  the 
TDRL.  However,  a  monostatic  antenna  more  directional  than 
the  present  omnidirectional  radiator  presently  used,  may 
make  the  system  more  generally  applicable.  Improvements  in 
reception  by  better  match  of  the  systems  antenna  components 
should  allow  for   the   determination   of  finer   feature   on 
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smaller  targets.  A  more  rapid  rise  time  generator  having 
greater  power  outputs  than  presently  supplied  by  the  UHF 
impulse  generator  would  be  useful.  Finally,  a  fully 
computerized  analog  digital  processor  with  a  microprocessor 
controlled  sweep  rate  and  distance  windowing  function  would 
add  much  to  the  automation  and,  consequentially,  consistancy 
of  results  in  the  TDRL. 

This  project  is  dynamic  and  holds  promise  for 
potentially  practical  applications  of  applied  research. 
There  remain  many  areas  to  explore  by  using  the  OS 
imaginatively.  The  system  is  now  ready.  The  true  need  is 
for  the  refinement  of  procedures  in  using  the  OS  and 
applying  it  to  research  areas  having  merit  both  instructive 
and  practical. 
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APPENDIX  A 
DESCRIPTION  OF  PERIPHERAL  GRAPHIC  3Y5EM  ROM 

As  stated  in  Chapter  IV,  three  peripheral  ROM  pacics  have 
been  added  to  the  Graphic  System  to  tailor  it  to  the 
specialized  needs  of  the  TDRL.  Thsy  are  the  Signal 
Processing  ROM  Pack  No.  1,  the  Signal  Processing  ROM  Pack 
No.  2  (FFT)  and  the  editor. 

Ihe  Signal  Processing  ROM  Packs  provide  fifteen 
mathmatical  operations  that  are  essential  to  advanced  signal 
processing  and  data  analysis.  Signal  processing  is 
performed  on  waveforms  or  lata  stored  in  one-dimensionsal 
arrays.  These  capabilities  do  no  alter  the  operation  of  the 
3S  or  occupy  any  of  its  available  RAM  space.  A  summary  of 
each  ROM  command  follows.  [Raf.  2H  25] 
A.   Signal  Processing  ROM  Paok  No.  1 

1.  MAX  --  finds  the  maximum  of  an  array. 

2.  MIN  —  finds  the  minimum  of  an  array. 

3.  CROSS  --  finds   the  location  of  a   specified  crossing 
level  within  an  array 


214 


4.  DIF2  —  performs  a  two-point  differentiation  of  an 
array . 

5.  DIF3  —  performs  a  three-point  dif f arentiation  of  an 
array . 

6.  INT  --  integrates  an  array. 

7.  DISP  —  displays  a  graph  of  an  array  in  raw  form 
(without  axis) . 

B.  Signal  Processing  ROM  Pack  No.  2 

1.  FFT  —  computes  the  fast  Fourier  transform  of  an 
array,  placing  the  results  in  the  same  array. 

2.  IFT  --  computes  the  inverse  Fourier  transform  of  an 
array,  placing  the  results  in  the  same  array. 

3.  CONVL  --  convolves  cwo  input  arrays,  placing  the 
results  in  a  third  array. 

4.  CORR  --  correlates  two  input  arrays,  placing  the 
results  in  a  third  array. 

5.  POLAR  —  converts  an  array  of  FFT  data  from 
rectangular  form  (reals  and  imaginaries)  to  polar 
form  (magnitude  and  phase) . 

6.  TAPER  —  Multiplies  an  array  by  a  cosine  window  of 
program-selectable  tapering  weights.  when  tapering 
is  selected  to  be  50%,  it  provides  a  Hanning  window. 
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7.  UNLEAV  --  sorts  an  array  of  interleaved  FFT  data  into 
two  arrays,  one  containing  real  and  one  containing 
imaginary  components. 

8.  INLEAV  —  interleaves  the  real  and  imaginary  data 
from  two  input  arrays  into  a  third  array  whose  format 
is  acceptable  to  the  IF!  command. 

The  EDITOR  is  a  ROM  ievice  that  contains  firmware 
routines  that  allow  ASCII  magnetic  files  to  be  altered  or 
edited  off-line.   [Ref.  26] 
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APPENDIX  B 

GENERAL  IMQSMATION  --  IDRL  OPERATING  SYSTEM 

A.   GRAPHING  VARIABLES 

The  following  named  variables,  numeric  and  constant,  are 
used  in  more  than  one  subprogram.  A  brief  description  of 
their  purpose  and  the  file  in  which  they  are  initiated  is 
given  below.  Note  specific  program/subroutine  for  local 
numeric  variables,  constants  and  arrays. 


r 1 

VARIABLE 
NAME 

PURPOSE 

1 

r      1 

FILE 

INITIATED 

i- 

STRING  CONSTANTS 

1 

SS 

DATE  OF  RUN 

2 

7$ 

TARGET  DISTANCE 

2 

G$ 

ANTENNA  NUMBER 

2 

H$ 

TARGET  NUMBER 

2 

m$ 

VERTICAL  SCALE  FACTOR  - 

mVOLTS 

2 

1: 

?$ 

i  _  „,.. . 

RUN  NUMBER 
_  ■  - 

1 

t 

2 
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1 —     — 

VARIABLE 
NAME 

PURPOSE 

.... 

FILE 
INITIATED 

R$ 

"RUN" 

2 

HORIZONTAL  SCALE  FACTOR  -  SECONDS 

-> 

is 

■  ■  ■          ■  ■  ■ 

SPECIAL  REMARKS 

2 

NUMERIC  CONSTANTS  AND  VARIABLES 

■"   i 
36 

OPTIMIZATION  COMPENSATION 

3 

r         i 
37 

OPTIMIZATION  COMPENSATION 

3 

B8 

OPTIMIZATION  COMPENSATION 

3 

C5 

TARGET  LENGTH 

2 

C6 

TARGET  RADIUS 

2 

—   —   —  —  ■ 

C7 

ANTENNA  PARAMETER  AVERAGING  VAL 

i 

2 

L 

OPTIMIZATION  FACTOR 

3,4 

L1 

FLAG 

3 

S2 

SIGNALS  AVERAGED/WAVEFORM 

"  ■ 

2 

P7 

FLAG 

2 

....... 

P8                    RUN  NUMBER 

I 

2 
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r       —  i 

VARIABLE 

NAME 

j 

r           —                      —  -  1 

PURPOSE 

1 

FILE 
INITIATED 

27 

OPTIMIZATION  COMPENSATION 

3 

1     ...-.—  . 

ENDPOINT  DISTANCES  OPTICS 

3 



51 

r  '  ' '   '  '  ' 

MAG.  TIMS  SCALE  FACTOR 

2 

S4 

FLAG 

2 

.  .... 
55            ANTENNA  PABASiSTSB- LENGTH              3 

1 
S6 

ANTENNA  PARAMETER- WIDTH 

3 

n 

DPO  SWEEP 

2 

W 

MISC  CONSTANTS  ARRAY 

2,3 

Z1 

j 

SAMPLING  RATE  DPO 

2     1 

i 

,.  _.  , j 

APPENDIX  C 
LIST  OF  INPUT  PROGRAM 

A.   GENERAL  DESCRIPTION 

1.  Location:  File  2 

2.  Memory  space  required  {less  mass  storage): 
24552-bytes 

3.  Maximum  no.  of  steps:  272  (212  less  remarks) 

4.  Description:  INPUT  provides  the  means  to  initialize 
and  drive  the  TDRL  Operating  System.  It  allows  input 
of  all  waveforms,  averages  acquisitions  and  stores 
data.  It  allows  for  input  of  antenna  target 
parameters  and  initializes  most  graphing  parameters. 
Inputs  are  from  the  keyboard  of  the  OS,  mass  storage, 
and  File  4.  Output  is  to  mass  storage  and  File  3. 
This  program  is  highly  operator-processor 
interactive.  Numerous  messages  and  directives  guide 
the  user. 
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B.   LOCAL  PROGRAM  ?ARI ABLES  AND  CONSTANTS 


1 •   Significant 


TYPE 


ARRAY 


1 — 


NAME 


-" 


USE 


DIWAV 


REMARKS 


(512) — TEMPORARY  ARRAY  FOR 
DIRECT  WAVEFORM 


ARRAY 


ARRAY 


AUGWA V 


(512) — TEMPORARY  ARRAY  FOR 
AUGMENTED  WAVEFORM 


-t- 


XO 


ACAV3 
INWAV 
DIWAV 
AUGWA V 
SCAT  WAV 


(512) — TEMPORARY    ARRAY    FOR 
SIGNAL    ACQUISITION 
AND    AVERAGING 


CONSTANTS 


?0 


ALL 


MENU    ITEM 


2.      Miscellaneous 

a.  ARRAY:   BO  (512) 

b.  SIRING  CONSTANTS:  AS ,  33,  C$,  D$,  J3,  L$,  K$ 

3  ir  YS 
C.   NUMERIC  VARIABLES  AND  CONSTANTS:  C1,  C2,  C3 ,  C4, 
C5,  I,  MO,  N,  N1  ,  P9,  50,  T,  X 


221 


C.   MEMORY  REQUIREMENTS 

1  •   10.12.  Internal  Data  Storage 


DESCRIPTION 


EACH  ACQUIRED  WAVEFORM  TO  BE  AVERAGED 


AUGMENTED  WAVEFORM  ARRAY 


DIRECT  WAVEFORM  ARRAY 


i 1 

BYTES 


41  14 


1 


4114 


41  14 


VARIOUS  NUMERIC/STRING  VARIABLES  AND  CONSTANTS 


5575 


TOTAL 


17917 


■ 


MAXIMUM  MEMORY  REQUIREMENTS  (INTERNAL)  FOR  INPUT 


42469 


2.   External  Mass  Storage  Requirements  (Data  Ta£e) 


DESCRIPTION 


T 


BYTES 


1 


FILE    # 


IL 


DIRECT    WAVEFORM    ARRAY 


INCIDENT    WAVEFORM    ARRAY 


4114 


50 


4114 


51 


AUGMENTED  WAVEFORM  ARRAY 


BACXSCATTERED  WAVEFORM  ARRAY 


4114 


52 


4114  J   53 

I 

-f 


TOTAL 


16456J 
I 
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D.       SUBROUTINES 


r  ■■      "i 

CLASS 

.,  ., 

NAME 

LINE 

■— n 

STEPS 

-i 

r  — i 

DESCRIPTION 

DRIVER 

DRIVER 

100- 
1380 

129 

40 

-  ■ 

INITIALIZE  AND  DRIVE 
PROGRAM.  PASS  CONTROL 
TO  NEXT  FILE 

INPUT 

ACAVG 

2000- 
2730 

74 

23 

ACQUIRES,  STORES  AND 
AVERAGES  DIWAV,  INWAV, 
AUGWAV  INPUTS 

INPUT 

INWAV 

3000- 
3170 

18 
19 

6 
6 

INITIALIZES  INPUTS. 
STORES  INCIDENT  WAVE 

INPUT 

DIWAV 

3200- 
3380 

19 

6 

INITIALIZES,  INPUTS, 
STORES  DIRECT  WAVE 

INPUT 

AUG  WAV 

3400- 
3590 

20 

6 

INITIALIZES,  INPUTS, 
STORES  AUGMENTED  WAVE 

INPUT 



SCATWAV 



3600 
3630 

9 

3 

INPUTS  TARGET  SCATTER 
FROM  DIFF.  OF  INWAV  - 
AUGWAV.  STORES  DATA 
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. ,   ,-,.-,■■ 
CLASS 

NAME 

LINE 

STEPS 

1 

% 

DESCRIPTION 

SPECIAL 

ANTLOC 

4000- 
4060 

7 

2 

1 

INPUTS  ANTENNA  PARA.  | 
FROM  MASS  STORAGE 

SPECIAL 

IGTLOC 

4100- 
4150 

7 

2      TARGET  DIMENSIONS 

(IF  KNOWN)  FROM  MASS 
STORAGE 

1 
SPECIAL 

NEWPAR 

4200- 
4330 

r  ■—  ■"■"  ■  —  ■  1 
4 

■  ■ _  —  " 

1 

.  ....  _  . . i 

INITIALIZES  ANTIN  AND 
TGIIN  FOR  STORAGE  OF 
NEW  ANTENNA  AND 
TARGET  PARAMETERS 

SPECIAL 

ANIIN 

4400- 
4600 

21 

7 

INPUTS  NEW  ANTENNA 
PARA.  TO  MASS  STORAGE 

SPECIAL 

TGTIN 
_. 

4700- 
4800 

L    _  _ 

1  1 

p 
3 

INPUTS  NEW  TARGET 
PARA.  TO  MASS  STORAGE 
_              j 
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APPENDIX  D 
LIST  OF  HATH  PROGRAM 

A.   GENERAL  DESCRIPTION 

1.  Location:  File  3 

2.  Memory  space  required  (less  mass  storage) : 
14544-bytes 

3.  Maxinum  no.  of  steps:  206  (169  less  remarks) 

4.  Description:  MATH  performs  the  significant 
mathematical  computations  for  the  operating  system. 
These  are  fast  Fourier  and  inverse  Fourier 
transformations,  determination  of  impulse,  step,  and 
ramp  responses,  noise  optimization  and  compensation, 
and  antenna  parameter  determination.  Inputs  are  from 
File  2,  File  4  and  mass  storage.  Output  is  to  File  4 
and  mass  storage.  This  program  is  competely 
automated  with  no  operator-processor  interactions. 
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B.   LOCAL  PROGRAM  VARIABLES  AND  CONSTANTS 


1 .   Sign  if icant 


r   —  " i 

TYPE 

—  ■■  l 
NAME 

■  —  --   — r   — i 

USE        REMARKS                             1 

ARRAY 

30 

OPTIC 
DRIVER 

(512) --METRIC  PHYSICAL  OPTICS 
FROM  RAMP  RESPONSE 

ARRAY 

E 

DRIVER 

(257) --FREQUENCY  DOMAIN 
BACKSCATTER  MAG. 

L  .                      .  .  ,. . 

ARRAY 

EO 

'"    ■    i 
DRIVER 

(257) --FREQUENCY  DOMAIN 

IMPULSE  RESPONSE  PHASE 

1 
ARRAY 

E3 

DRIVER    (257) --FREQUENCY  DOMAIN 

BACKSCATTER  PHASE 

ARRAY 

F 

DRIVER 

(257) --FREQUENCY  DOMAIN 
DIRECT  WAVE  MAG. 

ARRAY 

FO 

DRIVER 

(257) — FREQUENCY  DOMAIN 

IMPULSE  RESPONSE  MAG. 

ARRAY 

F3 

DRIVER 

(257) --FREQUENCY  DOMAIN 
DIRECT  WAVE  PHASE 

ARRAY 

10 



DRIVER 

(512) --TIMS  DOMAIN  TARGET 

IMPULSE  RESPONSE 
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h 1 

TYPE       NAME  j USE 

1 
1 

1 

REMARKS 

ARRAY 

RO    (DRIVER    (512) --TIME  DOMAIN  RAMP 

RESPONSE 
i 

1     "  "           ' 

CONSTANTS  MISC  | ALL 

1 
1 

1 

... 
MISC  CONSTANTS  AND 

VARIABLES 

L             _  j 

2«   Miscellaneous 

a.  ARRAY:   XO  (512)  ,H  (257)  ,P  (257) 

b.  NUMERIC  VARIABLES  AND  CONSTANTS:  I,  13,  K,  M1,  H2, 
N,  N1f  P1,  P2,  Q,  31,  Q2,  Q4,  32,  S2,  U 

C.   SEaORY  REQUIREMENTS 

1  •      {£0.52    Internal    Memory   Requirements 


1 —        -     -         -   -                   1 

DESCRIPTION 

i -  i 

BYTES 

EACH  ACQUIRED  WAVEFORM  TO  BE  AVERAGED 

41  14 

FREQUENCY  DOMAIN  BACK  SCATTER  PHASE 

2074 

FREQUENCY  DOMAIN  BACK  SCATTER  MAG 

2074 

FREQUENCY  DOMAIN  DIRECT  WAVE  PHASE 

2074 

FREQUENCY  DOMAIN  DIRECT  WAVE  MAG 

2074 

FREQUENCY  DOMAIN  IMPULSE  PHASE 

207U 
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DESCRIPTION 


FREQUENCY  DOMAIN  IMPULSE  MAG 


BYTES 


:07u 


EACH    ACQUIRED    WAVEFORM    TO    BE    AVERAGED 


FREQUENCY    DOMAIN    TEMP    STORAGE    PHASE 


in  14 


2074 


FREQUENCY  DOMAIN  TEMP  STORAGE  MAG 


TIME  DOMAIN  TARGET  IMPULSE  RESPONSE 


TIME  DOMAIN  RAMP  RESPONSE 


PHYSICAL  OPTICS  FROM  RAMP  RESPONSE 


S CATCH  PAD  TEMPORARY  STORAGE 


VARIOUS  NUMERIC/STRING  VARIA3LES  AND  CONSTANTS 


:074 


41  14 


4114 


4114 


4114 


1516 


TOTAL 


34564 


MAXIMUM  MEMORY  REQUIREMENTS  (INTERNAL)  FOR  INPUT 


49108 


:3  3 


D.       EXTERNAL    MASS    STORAGE    REQUIREMENTS     (DATA    TAPE) 


DESCRIPTION 

r            ■  i 

BYTES  FILE  #| 

I 

FREQUENCY  DOMAIN  IMPULSE  RESPONSE 

4148 

56   | 

■  "    "  '■■'■"  i 
FREQUENCY  DOMAIN  DIRECT  WAVEFORM                4148 

55 

FREQUENCY  DOMAIN  BACKSCATTER 

4148     57 

TIME  DOMAIN  IMPULSE  RESPONSE 

4114     54 

PHYSICAL  OPTICS  FROM  RAMP  RESPONSE 

r 
4114 

58 

TIME  DOMAIN  RAMP  RESPONSE 

4114 

59 

GRAPHING  PARAMETERS 

8  36  4 

60 

TOTAL 

29036 

. 
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E.   SOBROOTINES 


T 


CLASS 


DRIVER 


NAME 


■t 


D  RIVER 


FFT 


T 


OPTIC 


ANTUP 


oprm 


LINE 


100- 
1210 


STEPS 


122 


2000-    15 
2140    I 

(_ 


2200-   |51 

2700    | 
I 


3000- 
3080 

J__. 


3100-   |9 
3130    I 


L 


5  9 


25 


-r 


DS3CRIPTION 


INITIALIZE  AND  DRIVE 
PROGRAM.  PASS  CONTROL 
TO  NEXT  FILE 


PERFORMS  FAST 
FOURIER  TRANSFORM 


SEPERATES  TARGET 
FROM  RANGE  NOISE 


-f 


UPDATES  ANTENNA 
PARAMETERS 


COMPENSATES  PHYSICAL 
OPTICS  OPTIMIZATION 


240 


■  • 
•  • 


•  »  ••  ••  •• 

•  •         •• 

•  •  •-»  *  «   t_ 

•  •00  ••  h) 
••!••> 

•  •  «^  ••  •** 

•  •  o  ••  &. 

••  I    ••'^ 

•  •09  •• 

•  •  ••  £ 
......  ^ 

e 

......  e» 

......  o 

£- 

......  a 

••     » • 

••x  ••  ff 

.  .  H"  .  .  •»* 

« » <X  •«  4 

•  •X  ••  £ 

•  •  •  • 

......  iy 

....  ..  N 

«.  «•  •«  ••* 

.•••••  d 

.«   ..   *•  ••« 

..  .«  •••f 

..        ..  •*< 

••ro ••  c 

•.     ••  •.< 

..  ui .. 

•  •  __i  *• .. 
..  •— «  •• .. 
..  ti. » . .. 

•  •      .... 


•• 

<«. 

•• 

rw 

•• 

in 

•• 

CM 

•• 

«-< 

w 

•• 

1- 

m 

•• 

•k 

u. 

•• 

rw 

m 

•• 

u 

/■N 

•• 

■k 

r*. 

iO 

m 

c_ 

u 

CM 

V 

•k 

^ 

+* 

in 

ro 

++ 

o 

Ui 

o 

■. 

•> 

u 

CM 

y% 

>A 

X" 

rw 

JC 

■k 

in 

»J 

«H 

IM 

d 

M 

\* 

J2 

■k 

a. 

rw 

m 

-M 

o. 

/N 

V 

» 

i^ 

9» 

09 

in 

(. 

a. 

CM 

O 

•• .. ..  SI  CD  ii 
uJuJUJOi-l-J 

a:  oca: 

OCDO 

r-  co  o*i  ©  ~*  cm 


CD 
CM 


OJC9 
X«^ 
H-  •• 

II  <2* 
T 
'COO 


4* 


X_l 

4*-« 

X_J 

OOO 

4*  to 

U.C0 

uico 


~u. 


•  •  •  • 
T       CM  CM 

3 

xcd<z<z 

•-<   .1  LULu 


a 

o 

00 

coco 
o  * 
(\jm> 

CO 

z  - 

UJfw 

xao 

H-  •• 

CM 

It  &> 
0.0 
U.U1 


fw 

m 
u. 


cut 

>•  *.• 
ox 
u.  - 

rw-« 

mm 

CM~ 


4. 
O 

n 


i. 

O 


•  • 

•  • 

•  • 

•  • 

•  • 

•  • 


•  • 

•  • 


•  • 

•  • 
•• 

•  • 

V 

> 
c* 

3 


w 


£ 

o 

^*- 
W 

C 

d 
i. 


i. 
o 


/•».  /•>     •  •  ro  o  . .  o  cd 

r*-<M     ..inxo         Mtfjxo 

U1^       ......'O  ......vo 

ojin     ..oic\ico         ••cicvico 

^•^  ..iMOICVJ  ••OJCVJCVJ 

QO  • «  o>  'i*  » »  liHiJ 

LUX                         GO  ..  CO 

<>J  ..Cic2»3  CL..OOD 

—isxoxxoio^x  n  xx<rtnx 
II  hh»-i  u  Liit-«LUO  ii  roui«-«uiO  M 

—icoxuiU.ixoujuiuiU.acou. 


•• 
•• 
•• 

G 
O 


<3 
H 


a 
o 

« 

C 

o 

a 

i. 

* 

a 
c 


0 

s 
o 

-a 

3i 
u 

c 
« 

« 
U 
*♦. 

.. 

•• 
•• 
•• 
•• 
•• 

«"««U.CU 

a.-**  k.uiu. 
H  xx  n  M 
m  n  ujoco 

U.O.0TUJU. 


CO  >S  O  ••""»  CD  CD  CD  CD  CD  CD  CD  ^^  CD  CD  CD  CD  CD  CD  CD  CD  CD  CD  CD  CD  O  O  CD  CD 

ro^mmvofwCoc^cD«M(\jcviro^rm^>r>.cocriCD^<v4rn^rmvorwoo 


24  1 


m 

♦    N.  +  aa 

O  O       P*l  •• 

U.  UlUiUl  ••       *"•* 

n  n  u  «  zoc-? 

U_  UJ  Ql  U.  QC  U-  Z 


P 

o 


0 

c 
I 
o 

¥1 

s. 


c 

V 
O 


> 


O 


-a 
* 

N 

•wi 

a—         <«w 

O  /"». 

-«  -> 

■**  %• 

•w*  is 

C      U. 

•a*  Vrf" 

♦OCT) 

••too 

...MO 
•  •       ■*■ 

••o  /*• 
..►-"-> 

■  ■       -w 
•>  H  LU 


■  • 

■  • 


•  • 
•a 

•  a 


Of 
in 

c 
o 
a 
m 
v 

a. 
« 

VI 

a 
s 


-> 

u. 
z 
(ft 


o 

uj-) 

il 

»•  txi 
O.Z 


a 
c 
o 

CD"D 

-  w 

a.  s 

3  ..Min 
<x  ••  «.\ 

Ui  »•  =      '"N 

_J •  •!— CD 

CD     q  z  •  •  u.  *■*  ro 

U         II  «-»•  .twv  •-• 
O  ^  ©  '"■-  *     ••  ■     21    I 
«  N   II  N        -a       ^)«3 

~«CM~«<M-JZ_i  n    II 

vwv<iuj<iroo 

£ra.a.uaoMM 


•  a 

aa 
a« 


aa 


C 

o 

Wl 

1= 


•a* 

s. 
o 


V 

c 


<M  C 

-«  o 

o 

QQ  aa 

a.  aa 
/•S  a  a 

o_cm  .. 
zm  •• 

v  *• 

UJ  o  »• 
H-QC  ••© 
Ui        »•  <— • 

JSZ  ft 

UJ'—UJO 


•  a 

•  a 


•  • 
aa 

•• 
a* 
aa 
aa 
aa 
aa 
aa 
aa 


a* 
aa 


aa 
aa 


aa 
aa 


•  a 
aa 


n.     po 

u. 

<s>  z      - 

t- i-iuitnu. 

ID  ••XU")  •• 

•  «Oi>-  ..cm 

CM  Oi  CM  CM 
CVi  (»;  ^  oi  G; 
(3)        n  <±> 

Lurw     uj 

OH-O.OH- 

Z  **       Z  •"< 

U.Z»->U.Z 


ui 

'£>Ui 

in  •• 

•  •CM 

CM  CM 
CMC2J 
C*> 
Ui 

ai- 

z  ►•* 


o 
u. 

CO 

r»ui 
in  •• 

..CM 
CM  CM 
M(2) 
Of 

Ui 
OH- 
Z«-« 

i-.cc: 


M 
CD  a 
U.       O 

»  a 
O  w 
ui      «j 

•     t. 
CD 
-•       Q, 

-CD  V 
rOiTl+'CD 
U.CD  ttX 

r*>     «%co 
uiz—dc 
..UI  ■•  » 
U.X  ••« 

UI      *«z 

ui  a  ••« 
*-rw  %a 

uitv.  ••-! 

-i  S-J 
UiU.Ui<C 
o*-«ot:<_> 


CD  O  CD  *S  CD  O  CD  CD  CD  *33  CD  CD  CD  CD  CD  (S  t9  O  CD  CD  <S  CS  CSJ  CD  CD  <S  CD  O  CD  CD  CD  CD  CD  CD  CD 

aso^cMnn,uD^^oD<T»o^cMr^^«nvDoo<y(o^cvinTru^u3rw03cr«cD»-«cMf,r)'^ 
ro^^^^^v^^^tru^u^toinu7inioininvijvDv^ciDVD^^vDvovor^rvrwrwrw 


24 


•• 

•• 

■  • 

V 

•  • 

•  • 

•• 

•• 

%• 

u 

•  • 

•  • 

•  • 

•• 

•  • 

cs 

•  • 

•• 

■  • 

•• 

•• 

9 

•• 

•• 

•• 

•• 

•  • 

4* 

•  • 

•  • 

•• 

•• 

•• 

* 

•• 

•  • 

•• 

•  • 

—4 

•  • 

M 

•• 

•  • 

•  • 

"d 

•  • 

u 

•• 

•• 

•  • 

•  • 

•9+ 

•• 

•  • 

•  • 

V 

•  • 

** 

•• 

•• 

•  • 

4* 

•  • 

a 

•• 

•• 

•  • 

0 

•  • 

o 

•• 

•• 

•  • 

M 

•  ■ 

•• 

•  • 

•  • 

p 

•  • 

c 

•  • 

•« 

•  • 

w 

•  • 

o 

■• 

•• 

•  • 

a 

•  • 

•  • 

•  • 

•  • 

c 

•  • 

VI 

•  • 

•  • 

o 

•  • 

4* 

W 

•  • 

•  • 

u 

•  • 

u 

4* 

•  • 

•  • 

\ 

•  ■ 

01 

d 

•  ■ 

■  • 

\ 

•  • 

*♦» 

91 

•• 

•  • 

4* 

•  • 

>*. 

&. 

•  • 

•  • 

d 

•  • 

d 

0 

•  • 

•  • 

91 

•  • 

A+ 

•  • 

c 

•  • 

e 

V 

•  • 

9 

o 

I 

« 

•  • 

4* 

IA 

•»« 

o 

e 

•  • 

e 

4* 

H. 

0 

•  • 

*v 

d 

o 

ft 

•• 

a 

4* 

N 

VI 

M 

•  • 

d 

•94 

&. 

«f 

•  « 

M 

c 

s 

V 

c 

•  • 

u 

o 

••4 

4-> 

•  • 

•m* 

c 

4* 

d 

a 

•  • 

4* 

o 

a 

c 

c 

•  • 

a 

a 

o 

•3 

o 

•  • 

0 

o 

9 

£_ 

LCD 

CM 

c. 

V 

»— 

on 

e 

o 

a© 

<^ 

"O 

vt 

9 

|S 

**i 

•■4 

is. 

(S 

c 

u 

z 

V 

9 

04 

d 

a 

••4 

Ul 

4* 

V 

cz 

■ 

a. 

M 

X 

c 

+» 

PU1 

v4 

u 

w 

31 

»- 

d 

« 

vz 

cu 

\ 

•^ 

V 

JZ 

in 

■^H 

<- 

VO 

4» 

L 

a 

s 

V 

e 

(= 

<*N 

N 

a 

<«■• 

4* 

V 

o» 

^ 

rw 

>• 

*■ 

o 

a 

HH 

s 

e- 

0 

1—1 

a 

**• 

+ 

Ol 

\ 

♦ 

E 

«h 

(S 

e- 

"»3 

-co 

c 

WN 

/^» 

• 

ID'O 

>^ 

V 

OCDCM 

•«4 

s 

a 

cor* 

o 

>e- 

«¥ 

•^ 

C0C0-»h 

v^ 

> 

S.C£Q£ 

vw 

II 

•a 

0QCN 

u 

OS 

UJ 

\ 

♦  +  +L3 

o 

CM 

■1 

1      »s 

4* 

■k 

m  it 

\^ 

«^ 

vo 

/>./>.  /"N  >• 

M 

<m* 

*  •  o  o 

■  • 

z 

•  • 

CDZ 

•  • 

«*PO 

-J 

4- 

NCJP0^< 

_im 

9 

to 

■•XX 

•  • 

•  • 

CQUi 

^s 

■  • 

•  •Z-* 

<xinr* 

m^- 

^)  id  M)<x  <X 

••09  CD 

\ 

•  •     •■ 

■<• 

>•• 

QC 

•  • 

-z 

«-» 

•• 

•  ■       /~. 

/-•. 

ZXLUV 

*^  \s  ^ 

3  •• 

•  •If)  •• 

<«s 

■■  s 

s 

•>(9C 

•  •  CD 

s    H- 

s^ 

•  •  CD 

••oc^ro 

•  %  ■• 

II 

«  kommio 

•  •CM 

•  ■  .«  CVI 

CO 

••►- 

22 

■  ■CD 

•  •  (S) 

z 

z 

■  ■  CD 

..oiyjiMW 

1^ 

^-.  co  <r<c  <x  cm  (\j 

•  •CM  CM 

xro 

••z 

1— < 

••eg 

•  •O  '0 

t— «  ^4 

t 

•  •  w-l 

•  •       »v 

v^- 

<M  OJ  OJ 

^  II 

(i  «  a 

CMC=) 

•  •CM«» 

-•oc 

•  •  •— 1 

s 

•  •CM  XV 

••roco 

S/s. 

/-«. 

•  •1*0 

■  •CD_l«n  <2jC1> 

1 

I     />>  /^.  -'V  *^. 

a> 

•  •  Cfi» 

rr  i 

•  •  s 

■ 

•  • 

•*• 

•  •      ** 

s    V 

CM 

■  • 

•  •  II  <x 

<r 

r^  is.  is.  ca  ro  -< 

Ui 

••      U| 

30 

•  • 

•  •00  ^r 

•  •  Q0  CD 

•-•00 

s* 

•  •  00 

•  •T  Z> 

(00<1 

..Ci}— 

COX 

..-J 

-J 

••3(0 

•  •300— 1-J0Q 

T 

•  •10  CD 

••en  ii 

•*" 

z<x 

*>■  >^   Vrf" 

v^  Vrf-  »• 

2h 

zz»~« 

h   n 

—  _J_I 

ZCO 

Z  CO  11 

-J         M 

il 

ZCO   II 

z     r»- 

t—t 

^uimu^inioioio 

i-<  atiii  t—oc 

roouj<i 

cc 

LUOU.UiOCD^rU.USOMJJO'H 

uiu_<x£>u.Q£<x<x<x<xa:<xu.za:u.:z 

XXKU4 

—* 

xuaiUMQ        tuja« 

CD 

»^ 

<\jroTfto*£r».aooM3 

-«  CO  CO 

O  O  O  O  O  O  CO  O  CD  CD  O  CO  CD  CD  CO  CD  CD  CD  CD  ON  CD  CD  CD  CD  QOOOOOWflfl'f  IT) 

to  <r>  rv  oo  Ci  cd  »* 

00  ON  CD  *h 

OirO^ltfVOr^OOONONCDOl^CDCDCDGDCDCDCOCDCDCDCDCD 

r^r>-rwr^rs-COOOCOCOO*«CNCkiO>iO*iCr»i7>0*iCX»CT»<X»<-* 

243 


u 

■> 

o 

w 
o 

». 

CM 

Z 


M 

«. 

a. 

CO 

a. 


qc 

~3 


© 

GO 

U.OCQC 

w  ••  •» 

CViCMOO 

LUH- 

Ol— LU 

Z»-«_J 
*-»0£Lu 


Z_J 

♦*to 
uco 

»   «kVO 

<»•<?•.-« 

LUCO— « 

**-«z 

OO-COLU 

*o  •  ■  •  •  X 

•  -CVJCVif- 
AJCMCVJ 

o  II 

©<Z»>-H-a_ 

*r  m  q:  a:  u. 

COU-ZZh-i 


M> 


</) 


o 
vow 

OCT) 

iniC 

CM  If) 

z<x 

—  lO 
MM 

N2 
Q.    - 

00    •» 
CL_i 


acc/> 


a  cog; 

*»  o    •»   ** 
CO  t-^CO 

at     a_m 

*o  covo 
a.  a*  coco 


r- 
cq 

»  • 
CM 
CM 
CD 

LU 


•  ON 

•  ocas 

•  XLU 

•  LULUT 
■  »— h- 

•  LULUO 


N 
\ 

N. 
N. 

V 

\ 
CO  V 
CO  \ 

»\ 

s  N 
^  V. 
Z  V 


Ct.  LU  UJ  LU  •—• 


<XUJ 


\  S   \  N.  N.  S 

N.  \   \   S   \  \ 

N  V  \  S,  K  \  •••••• 

N  V.  >.  N  X  N.  •••••• 

V.   V.   X.   N.  \   \  •■  .... 

S  \  \  \  V.  \  • 

N.  \  \  \   \.  ..  ..  »• 

«\\\\\ 

a:  \  \  \  \.  \  ..  ••  •• 

L3  \  \  \  \  S  ..  ••  •• 

0  \  \  v.  \  \  ..  ..  •• 

QC  \   \  V.  \  \  ••  ••  •• 

Q.\\\\\  •«  «• 
\  \  \  \  \  »•  •■  V 

1  \  V  \  N.  \  —  ..  C 
O.  \.  \  N.   S.  \  ••  mm—t 

1VWN.  \ii.iV 
Q.  \  \  \  \  \  .»  ..  3 
U\\\NN.....O 
v.  v.  \  V  \  ••  ••  i. 
OWWVmm 
J—  \  \  \  V  \  ..  ..  S 

O  \   \  \  \  N.  ....  O 
UN  \.  N.  N.  N.  ••••»*■ 

i  \  s.  \      n.  ..  ••  in 

i  \  \  \  w  \  ••  ..  e 

<r  v  %  vzM-Ht, 

Q£\  \  \  »-i  V  Ll.  ..-^ 

O  V  V  N  Z>  \  I  •«  C 
QC  V  \  O  \  I  ..  « 
O.  \  QC  \  QC  S  LU  •••-« 
\  Ul  \  to  \  Z  .«  J. 
X  \  Z-  v.  Z)  \  •-•  ••  3 
»—  %>.»-«\.0^\H-..  O 

<i  n  on  "v.      \zj..'*-     a. 

S  *».  iZ)  V  z  v  o  ..  * 

Z\2\UJ\3'«d       CO 

Hi  VLU  \  00  \  V)  «^-       >C 

\        \        \       ..       CO    « 

v.  s  v  \  \  \  ..  ••  ..x- 

V.  \  V  V  \  \  ..  ..  —    ~Q£ 

V.  S,  V  \  \  \  ..  ..  ..s    <X 

V  N  X   \  \  \  ......  H»  _J 

S  \  S.  \  \  \  ..  ..  ..U.O 

\  \  X.  V  \  \  ».  •«  ■•U.OL 
\  \  V   \\  \mmm8    -    Z 
\  N  N  \  \  \  ..  ....  QC 

Lulu*" 

QCQC 


LULU 

acoc 


LULU 
QCQC 


lu  lu  uj  a:  •  z  lu 

QC  CC  £X  CJ  O  Gi 


••        •• 

••        •• 

«  " !!  v 

tt 

......  a 

......  o 

......  a 

......  m 

......  v 

j- 

......  9| 

......  c 

« 

......  c 

"""!•  c 
......  o 

e 

■  •  a.  a.^O 

..  .r  .*  of 

..  ....  VI 

......  p 

......  o 

......  a 

..      ..  w 
..ai  ••  v 
>*o  *.  c 
..  t— < .« 
.•I—  •••** 

». Cl.  ..  4f 

•»o ..  c* 

..I  ••  (. 

..   I    II  o 

•  •LU  ■••»* 

.  .  Z£  .. 
.  *  i— i  ..  Qf 

•  «h-  ••** 
..  Z>  *.  tf 

•  »0  «»  i- 

.•0£  ••  Q| 

..  ca  ••  & 

•  «C/>  *•  (A       CM 
....  *>* 

........     n 

••••••••     \ 

..   ..  ..  ..  '"N 

COCf. 

........        QC     • 

........     n-»cs> 

.. .. .. ».     z  ■♦♦ 

........     z>ro 

...... ..ocna 

zzzz  u  n  a  «* 

LU  LU  LU  LU  o  r-)  M  1 
tXQCQCQCaStaCOCZ 


CO  CO  CD  CD  CO  CD  CO  *S  CD  CD  CD  CD  CO  CO  CO  CO  CO  CO  CO  CO  CO  CD  CO  CO  CD  CD  CD  CD  CD  CD  CD  CD  CD  CD  CO 

^ois-cocr>o^cMr^^in'ONCo^o^oDcno^cv»roTruTvaNcvji3^cvjf^^woM)N 

OOOO^^^^^^^^^^CiCVJ<T»Cr\OCDCOCDCOOOO'^CVJCMCMCVJCVJCMv\lC\J 
^^^^^^^^^^^^^^•^^^^C4CMCMCMCVJCyCMCMCM^lCMCMCVICV»CMW<M 


2'4  4 


9 

© 

CM 

Z 
UJ 


•  • 
•• 
•• 
•• 
•• 
mm 

•  m 

irt 
u 


o 


o 


Z       <5> 
If) 

a.     ro 

UJ       CM 

c/>     z 
Ul 

<~       X 
Z      t- 

/•N 

H-  »^*  V 
©  /^ 

W  It    It  «-"  tf 
«■ <  i^  /•».  ©  /^. 

in     ^m^ 

«  g;^     >• 

-«OOU.© 
ZU.(jQ*-«C0 


oc 
a 
v 

>• 

o 

IX 

© 

ac 
•\ 

a 
a:/*. 

z 
a  i 
z^ 

<X^'© 

Q  ~** 

<~.  a:-* 

z  i  ~ 

i  /-ro 

~-*3 
©CM   I 

©as  i  ^ 
-*  •*  -^  z 
m^^  i 
cm:*:©:*: 
*.•  ct  »-• 
-^o©  it  © 

I  t-QL^QC 

:*:         ro  a 

II  OU.VW 


© 
© 

rr© 
CM«-« 

\ 

z~ 

UJZ 
X  I 

C©  ©  <"n 
CM  U  s*» 

ro  u 

M30 


••4 

O 

-2 


O 

a 

V 

(= 

•*■*  z     © 

£  is- 

«UUi       CM 


V 

■O 


•  « 


I 
t 

t—  s*- 

© 

It     II 

a 

00 

U.CQ 


UJ 


© 

in 
01 


v 

St 

(= 


w 

c 
o 

V 
91 

s. 
o 


13 

in 

<z 
w 
a. 


© 

CA 

10 

CM 


©        ©  © 

a  ©  cm  cm 

>-•  it      a  ^ 

©  /■».  o  o 

CQOI-H  »— h- 

U.QGUJ  It  OUi 
mQIUZ^UZ 


UJ 

X 


Z~« 

II 

u.~ 


o     z 
u     ui 

z 


o     ••«**>. 

CM        ••U.V 
CM        ••%•  /^. 

11  o  n  ««<rv 
1  cm     -*z  n 

It  vOvlUUJU. 

Z3U3o:n«-« 


CM 

© 

CM© 
♦w  ■ 
3M 

■*- 
O^ 

•*♦ 

3^ 

—  «  3 

n  a  1 

>■<     a 

Vr-O    H 

3U.0C 


245 


CVI 

*- 

M 

** 
CM 

QE 

+ 
M> 

rM3ZJ 
■*■•«• -W 

^ao 
^♦•^^ 
«-.  oo 

ctcc  CO 
+   il    u  o 

ti)  <~. /*w 

HOOHO 

a  oow^ 

=>££CQl=t3 


CI 

10 

o 


3-« 

II   + 

OC  II 


IS  O  CD  CD  CO  CD  O 

^  o  j  ro  t  to  '>d  r» 
u>  <o  ^o  vd  vo  ^x>  *o 

(M  CM  CM  CM  CM  CM  CM 


•  •        •• 

•  •       •• 

......  Of 

•  •  •• ••■f 

•• ....  o 

«T5 

......  a. 

......  3 

.....I  J. 

..  ....  V 

..  ..  ..-*-• 

•  01 

......  s 

......  (3 

•  e 

•• .. ..  o 
......  a 

•  •     .. 

..  fl_  *.  <3 

••3  *.  a 

«.h-  ••  c 
■•z  ■■  ftf 

..  <X  •  «  •*-' 

■  •  I  ••  (= 

..    I    .«  1^ 
.. LU  •  « 

•  •Z  •  •  0> 

•  ■•-*.*     C 

..  H"  «.  "»*J 
■«3  ■■  £ 

•  a  O  ••    i.  '— ' 

•»ac «.  a*   » 

..  ^  ..  iy  CD 

..  co  «.  TJ   » 

• «      ..     o 

/•k  ......aaCQ 

C3  ........    - 

>•  a.  ..  .a   aaS 

O  ......  ..2 

CD  ••  ••  ••  •  »•— « 

u         s: 

<^  OZ s 

-<r     or*.  ......     in 

at-z>  .. .. ..  .._jco 

^XH- SSI^ZI-J  II 
QLULULULUUJUiOCCM 

XlZiQCQLLX:0£.Q£U(S) 

CD  CD  CD  CD  CD  CD  CD  CD  CD 

co  (Tt  cd  o  »h  cm  ro  «*■  to 
CMCMCMPirororororo 


••     •• 

......  V) 

......  ft* 

..  a.  .«-** 
......   J. 

o 

"O 

......  a 

:* 

......  o 

-J» 

CM        ...... 

CO         ......r- 

«•  o 

is.  .. .. ..  c 

cj  ■••••••«« 

v  .. ....  w 

/~S  a  .....  -.4 

CM *. 

co      ..  -•  •>  o 
I       ...... 

/^.      . .      . .  o 

CD       ..JC  ••■*■' 

UJ  ..  cu  ..  ft 

CM  »«0  ••  O 

•  •»!••$- 

IjO  v*  *  ■   I    . .  9) 

CO  *  •  •  LU  •• 

•4-/^.  •  »  ~ZL  »  •  U 

N  A  .  «  >— I  ••  ■*-! 

CJ  ^*  aa|—  •••*> 

V    Nrf-  a.  ZD  ••    CV 

/-s  3       . .  o  ..  o  *-* 

ko  i      ..q:  ••        on      •» 

CO  «••»        ..uQ««t=       G>       CO 

ICM  ..ZJ  •••!*  N.  CO 

/-S    -w*  •  •  . J")     •  •     ii.  /"*.  ■» 

CMS       •»       ••  /^       CD 

03   •-'  aa  a.  ••  ••  -*  CD 

%•  ■»*  ••••••••  »w           ■. 

CO-^  ........  3  » 

COfM  ........  I  X 

<lv      ........  oo  ^  r^  <i 

\  \.     .. .... ..cqcmcos: 

yj'inz .w#»  z 

CJOCE  ••  ••  ••  ..UD3CD       or 

».,•  '^  Z3  .a  a.  ..  ..CO^'CQ-JZ) 

ii  h>-s:2:s:s:ii  n  n  _jh- 

lOUTUJUJULiJUjrvNOIUJ 
COCOCtOrCCOiCeCQOCOUQC 

CD  CD  CD  CD  CD  CO  CD  CD  CD  JS  CD  CD 
«JDP>-COCD^CMrO'rrLOM3rwCO 

CD  CD  O  »■*  ^^  •**  "■!  ^~*  •■*  *^  ^^  •■* 

rorororororoforoforopofo 


uin 

cct 

CO 


2'4  6 


00 

u> 

00 

CO 

GO 

00 

CO 

CO 

CO 

CO 

CO 

CO 

fw 

y> 

U3 

'■£> 

<o 

<>£> 

V£> 

^0 

<£ 

•iJ 

*o 

N. 

r^ 

IS. 

IS- 

P^ 

r>- 

p^ 

Pw 

Pw 

Pw 

COCOCOCOOOOOCOCOCOOO 

CO  «..0  *-D  VO  *£>  U3  * O  U3  ^D  "43  *■£> 
VJ3PwP>-PwP>-pwJS.PwPwJS.Pw 

Pw 


-*OJP«>'«rUf>U>PwOOO>G>-« 


H-l-f-h-   h>   f»   I-   h-   t-t—   I— 


<X 


<Z 

o 


<x 


<x 


<r 


<I 


> 
<r 


> 
oc 
<i 


<x 


> 
a; 

<x 


lal 

2   **   *— i   >— «   •— <   »-* 
00   00   00   <t^   i^ 


> 

<r 

2 


<r 


5- 


> 


00   0B   00 


2   2   2  Hi  ui  UJ  Ul  UJ IU  UJ  LU  111  <X 

*-*       •-•   *-"»22222222Z-i 
00   00   00 


OJ 


po-«rtr>v0p^cooMS>-«oj 

04 

PO 

•«• 

IP 

i<0 

rw 

GO 

CM 

<3> 

-* 

ojporopopoporopo^r^^r 

M 

OJ 

CM 

OJ 

OJ 

OJ 

OJ 

OJ 

PO 

PO 

po 

247 


CVJ 


a_ 

•x 
»- 

i 

i 

CO 

<x 
z 
z 

Ui 

»- 
z 
<x 

u_ 

o 

CO 

z 
o 


<x 

u 

o 


Ui 


00 

CO 

00 

00 

00 

00 

00 

00 

00 

00 

oo 

u> 

'4) 

'■£> 

*0 

vo 

to 

10 

10 

*0 

10 

<0 

rs- 

JN- 

Pw 

r- 

Pw 

Pw 

Pw 

Pw 

Pw 

p>- 

Ps- 

0000000000000000 
OOi0l0l0l0i0t0<0<0 

^0Pwf^ps.pwpwp<-p».rw 
p«- 


-WW 


m     k0 


z     z     z     z     z 
<x     <x     <r     <r     x 


00         CPl         CO  "•-•  CVJ 

<*•«  *>4  «-4 

I—     1^ —     t —  I —  •— 

z     z     z  z  z 

<x     <r     x  i  ox 


<x 


<r 

z 


<x 


0- 

QC 

<x 


<x 


QC 
X 


<x 


> 

QC 
<X 

z 


X 


> 

QC 

X 

z 


X 


QC 
X 
Z 


X 


<X 

z 


X 


<x     <r     <x 

H-  t—  h- 
XXX 
Ct       Ci        Ci 


>  >-  > 

QC  QC  QC 

<r  x  x 

z  z  z 


QC 
X 


XX 
<X>-H- 

(-<!<! 
XOC. 
Ck 

>0- 

>QCQC 

tt:<r<i 
xzz 


x<xx<x<xx 

t-hHHHH 

xxxxxa: 


QCQCOCQCQCOC 
XXXXX<X 

zzzzzz 


00       CO       00       00       ^^       ^ft       jjQ       00       <Tft       00       CQ       00 


igqcqoogqcqodcdoo 


u. 

■• 

•  •  CO 

z  •-• 

UJ  -J 

QC  H 


«      (U      M 


«n     <0 


oo 


CTi        —•        *-i 


rO"«Mn<0P>-OOO\® 


2U8 


APPENDIX    2 
LIST    OF    GRAPH    PROGRAM 

A.       GENERAL    DESCRIPTION 

1.  Location:    File  4 

2.  Memory  space  required  (less  mass  storage): 
44712-bytes 

3.  Maximum    no.    of  steps:    584     (473    less    remarks) 

4.  Description:  GRAPH  provides  the  means  to  visually 
display  the  results  of  signal  acquisition  and 
processing.  Subroutines  are  divided  into  three 
class  es: 

a)  GENERAL — those  subroutines  that  are  applicable  to 
most  of  the  other  routines  in  GRAPH.  They 
initialize  the  program,  format  graphics,  title, 
lable,    and   overlay   a   grid. 

b)  INVERSE  S3ATTERING--those  routines  applicable  to 
the  inverse  scattering  solution  of  time  domain 
electromagnetic  transients.  They  are  the 
secondary    drivers    for    the   development    of   graphics. 
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c)     FRONT'S      TECHNIQUE-- routines      applicable        to      the 
Prony    method    of    solution. 
Inputs   to   GRAPH    are   from   File   3,    File    7,    the  GS   keyboard 
and  aass   storage.        Outputs   are    to    the    GS    Z2T ,    mass   storage, 
File   2,    File    3,    File    6,    and   File    8.         This    program    is    highly 
operator-processor      interactive.  Numerous        messages      and 

directives    guide      the    user.         Signicant   graphs      are   provided 
for   operator   evaluation     of    processing   results. 

B.       LOCAL    PROGRAM    VARIABLES    AND    CONSTANTS 
1 •      Significant 


1     -      ■ 1 

TYPE 

NAME 

-i 

USE 

REMARKS 

ARRAY 

X 

VARIOUS 

(512) — TEMPORARY    STORAGE 
AND    PROCESSING 

ARRAY 

- 

Y 

VARIOUS 

(512) --TEMPORARY    STORAGE 
AND    PROCESSING 

ARRAY 

P 

—  ■ 

AXIS            |  (8)  —AXIS    LABLING 

TICKS          | 

1 

AXLAB          | 

1 1                                     1 
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IPS 


CONSTANTS  SO 


NAME 


USE 


ALL 


REMARKS 


FLAG 


CONSTANTS 


L$ 


ALL 


HORIZONTAL  AXIS  LABLE 


CONSTANTS 


1 

N$ 


ALL 


VERTICAL  AXIS  LABLE 


CONSTANTS 


PO 


ALL 


MENU  ITEM 


CONSTANTS 


-*- 


L2    (DRIVER 


OPTIMIZATION  FACTOR 


CONSTANTS 


AO 


GROP 


HORIZONTAL  SCALING 
FACTOR 


CONSTANTS 


A1 


GROP 


VERTICAL  SCALING 
FACTOR 


CONSTANTS 


I 


GROP     |  SCALING  VIEWPORT 

I 

j  COMPENSATOR 


2.   Miscellaneous 


a.   ARRAY:   none 


b.   STRING  CONSTANTS:  A$,  3$,  S$ 

C.   NUMERIC  VARIABLES  AND  CONSTANTS:  D5,  G1,  G2,  G3f 

G4,  C5,  I,  J,  Q,  P1,  P2,  93,    W,    71,  72,    V3,  VU , 

W1  ,  92,    W3,  WU,  H5,  »6 
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C.   MEMORY  REQUIREMENTS 

1  •   ^£0 5.2  I nt ernal  Data  Storage 


DESCRIPTION 


VARIOUS  TEMPORARY  STORAGE/PROCESSING 


MATRIX  LA3LING  ARRAY 


VARIOUS  NUMERIC/STRING  VARIA3LES  AND  CONSTANTS 


TOTAL 


| MAXIMUM    MEMORY    REQUIREMENTS     (INTERNAL)     FOR    INPUT 


BYTES 


3228 


82 


1566 


9866 


54573 


2.       External    Mass    Storage    Requirements    (Data   Tape) :    none 


D.       SUBROUTINES 


CLASS 

NAME 

LINE 

STEPS  \% 

DESCRIPTION 

DRIVER 

DRIVER 

100- 
1190 

110    |  19 

I 

I 

i 

I 
i 

INITIALIZE  AND  DRIVE 
PROGRAM.  PASS  CONTROL 
TO  NEXT  FILE 

GRAPH 

[_  -   —  ^ 
AXIS 

2000- 
2250 

l 

26     |  4 
I 

! 

,  i  . 

DRIVES  AXIS  GRAPHING, 

| GRAPH 

I- 

TICKS 

2300- 

i 

22     |  5 
1 

h 

AXIS  TICKS  INTERVALS 
h                — 1 
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CLASS 

NAME 

LINE 

STEPS  | % 
I 
i 

- 
DESCRIPTION 

GRAPH 

AXLAB 

2700- 
3190 

50     ! 

9 

LA3LE  VERTICAL  AND 
HORIZONTAL  AXIS 

1 
GRAPH 

TITLE 

3200- 
3310 

12 

2 

POSITION  GRAPH  TITLE 

GRAPH 

GRPAR 

3400- 
3870 

48 

8 

PRINTS  PARAMETER 
VALUES  AS  REQUIRED 

GRAPH 

GRID 

3900- 
4010 

"■"I 
12 

2 

DRAWS  GRIDS 

INVERSE 
SCATTER 

GRIN 

5000- 
5410 

42 

7 

GRAPHICS  FOR  DIRECT, 
INCIDENT,  AUGMENTED 
AND  BACKSCATTER 
TIME  DOMAIN  WAVES 

f- 

INVERSE 
SCATTER 

GRFOUR 

5500- 
6080 

■  -  ■■■ ■  ■"■" 

59 

10 

GRAPHICS  FOR  FREQ 
DOMAIN  PHASE 
AND  MAGNITUDE  OF 
DIRECT,  BACKSCATTER 
AND  IMPULSE  RESPONSE 

| , 

INVERSE 
SCATTER 

GROP 

6100- 
5410 

43 

8 
! 

GRAPHICS  FOR  PHYSICAL 
OPTICS  OF  TARGET 

- 

INVERSE 
SCATTER 

GRAMP 

6600- 
6910 

32 

I 

GRAPHICS  FOR  RAMP 
RESPONSE 
| 1 

253 


CLASS 


PRONY 


PHONY 


PHONY 


NAME 


PRIMP 


T 


PROTHAN 


-■ 


PROLES 


LINE 


7000- 
7340 


-f 


7400- 
7710 


7800- 
8300 
j 


STEPS 


35 


32 


f- 


51 
1 


+ 


DESCRIPTION 


GRAPHICS  FOR  IMPULSE 
RESPONSE 


j  GRAPHICS  FOR  PHONY 
jTRANSIENT  RESPONSE 


| GRAPHICS  FOR  PRONY 
I 


_L 


| POLES  OF  IMPULSE 
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APPENDIX  F 
XlilCAL  OPERATING  SYSTEM  PROCEDURE 

In  this  Appendix  is  presented  a  typical  run  in  which  a 
sphere  and  a  con?  are  illuminated  by  the  incident 
elect. omagnet  ic  field.  The  messages  presented  and  the 
results  displayed  are  those  as  would  be  encountered  during  a 
normal  run  of  the  system. 

The  OS  is  first  used  to  acquire  an  I  display  a  normal 
signal  (the  sphere) .  This  signal  is  then  optimized  to 
remove  noise.  The  results  of  the  effected  waveforms  are 
then  represented.  Finally,  the  system  is  returned  to  INPUT 
and  a  new  target,  a  cone  is  acquired.  The  results  for  the 
cone  are  presented.  Note  that  the  impulse  response  is 
optimized  by  a  factor  of  .5  while  all  other  graphs  for  the 
cone  are  shown  unoptimized. 
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